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Surtmary 

Ionospheric  modification  experiments  over  Piatteville,  Col.1,  and 

2  r 3 

Arecibo,  Puerto  Rico  are  operated  at  HP  intensities  sufficient  to  drive 
electrostatic  waves  unstable  in  the  F-layer  plasma,  through  the  process  of 
parametric  instability.^’^  We  have  studied  the  nonlinear  effects  ^>7, 8, A 
the  waves  excited  by  the  HP  modifier  on: 

(a)  the  modifier  propagation,  reflection,  and  deviative  absorption^, 

(b)  UHF  backscattering  models^’®;  and 

(c)  the  nonequilibrium  distribution  of  electrons  in  the  F-layerB. 

It  has  been  found  that  the  electron  plasm  waves  excited  by  the  Arecibo 
modifier  determine  the  general  features  and  magnitudes  of  the  obseived 
UHF  vertical  backscattering  cross-sections,  although  there  are  still  some 
unresolved  Issues.  The  most  widely  accepted  current  theory  is  based  on  ion- 
acoustic  frequency,  35  cm  density  fluctuations  beating  with  the  vertical 
component  of  the  modifier  electric  field  to  produce  density  fluctuations 

lib 

around  the  plasma  frequency.  The  scattering  cross-section  for  other 
frequencies  (VHF  to  Ghz)  and  other  scattering  angles  has  also  been  found6 *A. 
Mich  larger  cross-sections  are  predicted  at  scattering  angles  which  probe 
plasma  waves  propagating  nearly  parallel  to  the  geomagnetic  field6.  Below 
one  GHZ  the  scattering  model  is  based  on  the  saturation  of  plasm  waves  by 
induced  scattei a ng  off  ions  ’  (also  known  as  nonlinear  ion  Landau— damping^) , 
provided  photoelectrons  are  ignored.  Above  about  one  GHZ  the  scattering 
model  is  based  on  orbit-perturbation  saturation6  ,A.  The  latter  theory  is 
probably  also  relevant  to  scattering  below  1  GHZ,  when  photoelectrons  are 
taken  into  account,  although  this  calculation  is  still  in  progress. 

We  have  found  that  the  electron  plasma  waves  excited  by  the  Piatteville 
modifier  operating  at  a  peak  output  of  50  yui/m  vacuum  intensity  at  the 
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F- layer  cause  an  additional  10  to  H0%  nonlinear  deviative  absorption  of  the 

modifier,  for  a  Chapman  scale  height  of  75  km>  when  the  modifier  frequency  is 

below  f  F_.  At  1500°Kit  is  stronger  (20-402).  Photoelectrons  have  the 
o  2 

effect  of  reducing  the  size  of  the  nonlinear  altitude  range,  and  hence 
decreasing  the  anomalous  absorption.  We  therefore  expect  stronger  anomalous 
absorption  at  nigit,  compared  to  in  the  daytime.  The  anomalous  absorption 
which  occurs  at  higher  modifier  pump  powers  has  also  been  computed  and  may 
be  of  interest  if  self-focusing10  can  intensify  the  present  peak  power,  or 
if  future,  more  powerful  transmitters  are  developed. 

Hie  nonlinear  propagation  of  the  modifier  into  the  inhomogeneous 
ionosphere  has  been  studied0,  neglecting  the  geomagnetic  field.  The  modifier 
is  essentially  a  standing  wave  at  altitudes  within  several  kilometers  below 
the  reflection  point  even  when  there  is  strong  anomalous  absorption.  There 
are  steep  oscillations  found  in  the  Poynting  flux  pattern  and  hence  bands  of 
altitudes  at  which  local  anomalous  heating  should  be  strongest.  The  above 
theory  is  based  on  saturation  of  plasma  waves  due  to  induced  scattering  off 

ions. 

The  short  wavelength  electron  plasma  waves  excited  by  the  Platteville 
modifier  operating  at  peak  output  have  been  found  to  alter  the  electron 
velocity  distribution0.  We  have  found  the  rate  at  which  electrons  are 
accelerated  by  10-15  cm  wavelength  plasma  waves  to  energies  above  3eV.  The 
flux  of  such  energetic  electrons  along  the  geomagnetic  field  is  typically  of 
the  ordt  r  of  10^  electrons  cm-2  sec  1,  which  is  sufficient  to  explain  the 
6300  8  airglow  enhancements  observed  during  modification,  and  possibly  the 
5500  8  enhancements  as  well.  The  theory  is  based  on  perturbed- orb it 
saturation  of  electron  plasma  wavesA,D.  It  predicts  that  the  energetic 
electrons  originate  at  altitudes  about  122  of  the  Chapman  scale  height  below 
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the  reflection  point. 

Interesting  new  plasma  wave  effects  associated  with  the  geomagnetic 

field  have  been  found11' *  *  .  Platt eville  modifier  intensities  should  be 

sufficient  to  drive  a  new  parametric  instability  involving  both  very  low 

frequency  (“10”^Hz^  long  wavelength  (600  meter)  electrostatic  waves 

propagating  transverse  to  the  geomagnetic  field  and  backscattered  HP 
P 

radiation  .  A  crude  linear  stability  analysis  predicts  a  growth  rate  of 
several  minutes,  which  is  probably  too  long  to  account  for  the  wide  band 
attenuation  observed  in  ionosonde  experiments.  However,  given  enough  time 
to  develop,  these  wav.-s  would  constitute  field  aligned  irregularities.  A 
nonlinear  saturation  theory  has  been  developed  for  the  very  low  frequency 
field-aligned  waves.  Ihe  theory  predicts  a  fluctuation  level  for  these  waves 
which  is  about  1%  of  the  ambient  density. 


4 


I.  Introduction 


This  report  is  concerned  with  the  observable  consequences  of  various 

plasma  wave  instabilities  which  can  be  excited  in  the  F-layer  of  the  iono- 

1 

sphere  by  "punp"  radio  waves  in  the  5-10Mhz  frequency  range  having  vacuum 

p 

intensities  of  25-50  noi/ra  .  Three  parametric  instabilities  which  can  be 

7 

excited  are  examined  here  in  varying  degrees:  the  "decay  instability,"  ’  in  which 

an  electron  plasma  wave  and  an  ion-acoustic  wave  are  excited;  the  "purely-growing" 
11  B 

instability,  *  in  which  an  electron  plasma  wave  at  the  pump  frequency  and  a 
purely  growing  mode  are  excited;  and  "stimulated  diffusion  scatterlng"Fin  which  the 
excited  waves  consist  of  a  very  low  frequency  diffusion-like  mode  propagating 
across  the  magnetic  field,  together  with  a  backscattered  electromagnetic 
wave.  We  have  been  interested  principally  in  the  decay  instability,  as 
excited  by  an  ordinary-wave  punp  in  the  absence  of  a  magnetic  field,  for 
electrons  and  ions  assumed  to  be  the  same  temperature . ^  The  excited  waves 
produce  enhanced  electron  density  fluctuations  which  have  been  measured  in 
UHF  backscattering  experiments.  In  order  to  theoretically  conpute  scattering 
cross-sections,  study  the  effects  of  the  plasma  waves  on  the  electron  dis¬ 
tribution  function,  etc., it  was  necessary  to  know  the  nonlinear  plasma  wave 
Intensity.  The  determination  of  this  intensity  required  that  we  know  not 
only  the  linear  growth  rate  of  the  waves,  but  also  the  nonlinear  decay 
rate  which  is  ultimately  responsible  for  their  saturation.  A  conplete  theory 
further  required  a  knowledge  of  the  nonlinear  thermal  'Toise"  sources  of  the 
plasma  waves. Two  complementary  such  nonlinear  saturation  theories 
we»'e  developed  for  a  homogeneous  plasma.  The  first  involves  stimulated 
scattering  of  the  excited  plasma  waves  off  ions  into  less  unstable  regions 
of  phase  space  * ' ;  the  second  involves  absorption  of  the  excited  plasna 
waves  by  electrons  whose  orbits  have  been  nonlinearly  perturbed  by  ;he 

g 

waves.  We  shall  refer  to  the  first  of  these  as  "induced-scattering 
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saturation”  (sometimes  also  referred  to  as  nonlinear  ion  landau-damping) , 

and  to  the  second  as  orbit-perturbation  saturation.  The  induced-scattering 

saturation  mechanism  was  first  presented  by  DuBois  and  Goldman  at  the  Ivory 

Coral  Meeting  in  Albuquerque,  New  Mexico,  August,  1971.  A  similar  theory 

was  also  being  developed  independently  at  that  time  by  Valeo,  Oberman  and 
o 

Perkins."  The  orbit-perturbation  saturation  theory  was  first  presented  by 
Bezzerides  and  Weinstock  at  the  ARPA  meeting  in  Boulder,  Col.,  Noveirber,  1971. 
In  Section  II  below  we  describe  both  of  these  theories  and  further  refine¬ 
ments  by  the  present  authors  and  by  others.  Also  described  are  the  con¬ 
sequence.-  of  these  theories  for  the  plasma  wave  spectral  intensity  distrib¬ 
ution,  total  wave  energy,  and  for  the  associated  anomalous  resistivity. 

Section  III  below  is  devoted  to  effects  of  the  weak  spatial  density 
inhcmogene Ity  in  the  F-layer  on  the  plasma  wave  level,  and  the  observable 
consequences  for  absorption,  scattering  and  airglow.  For  the  ionosphere, 
a  sinple  model  is  found  to  be  adequate:  the  homogeneous  results  are  taken 
over,  but  the  density  and  hence  the  electron  plasma  frequency  and  principally- 
excited  plasma  wavelength  vary  with  altitude  for  a  fixed  pump  frequency. 

The  consequences  of  this  model  for  the  punp  propagation  towards  its 

reflection  point  are  explored  by  numerically  solving  the  nonlinear  punp 
r 

wave  equation  for  a  vertically  plane-layered  ionosphere,  neglecting  geo¬ 
magnetic  field  effects.  Sore  nonlinear  deviative  absorption  is  found,  and 
we  display  the  pump  electric  field  profile  with  altitude.  The  solution  gives 
us  the  electron  density  fluctuation  spectrum  as  a  function  of  altitude,  and 
hence  determines  the  UHF  scattering  cross-section  for  frequencies  in  the 
vicinity  of  the  electron-plasma  frequency  and  the  ion-acoustic  frequency. 

The  shorter  wavelength  plasma  waves  are  shown  to  be  capable  of  accelerating 
electrons  to  sufficient  energy  to  excite  6300  8  and  possibly  5500  8 
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colllsionally-induced  airglow.0 

Finally,  Section  IV  treats  an  interesting  Instability  associated  with 
the  geanagnetic  field;  a  process  in  which  the  pump  undergoes  stimulated 
scattering  from  a  very  low  frequency  (~10  ^  Hz)  diffusion- like  mode  propa¬ 
gating  across  the  magnetic  field  with  a  wavelength  given  by  one-half  the 

E  F  G 

local  pump  wavelength.  5  5  This  leads  to  very  slow  growth  of  both  the 
diffusion  mode  and  the  scattered  radiation.  Some  comparisons  are  made  with 
stimulated  Brillouin  scattering.  Preliminary  work  on  the  nonlinear  saturation 
of  the  "stimulated-diffusion  scattering"  instability  is  described.0 
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IX.  Homogeneous  nonlinear  saturation  of  the  parametric  decay  Instability 
ard  anomalous  resistivity 


a)  Definitions  ard  general  properties  of  nonlinear  saturation. 


The  Intensities  of  the  waves  excited  In  the  parametric  decay  Instability 
are  measured  by  the  correlation  functions  of  electrostatic  field  fluctuations. 


integrated  over  frequencies  near  the  electron  plasma  frequency  and  near  the 

6  7 
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ion  acoustic  frequency.  These  frequency-integrated  correlation  functions 


are  called,  respectively,  I^k.x)  and  I2(k,x),  the  former  referring  to 
electron  plasma  waves  (which  we  sometimes  call  Langnuir  waves ) ,  and  the 


latter  to  ion-acoustic  waves,  k  is  the  propagation  wave-number  and  x 


6  *k  is  the  cosine  of  the  angle  the  propagation  vector  makes  with  the 

—  rt  - 


pump  electric  field  (=  parallel  to  the  geomagnetic  field  for  an  ordinary 
wave  near  its  reflection  point).  For  equal  electron  and  ion  temperatures. 


the  Langjnulr  wave  intensity  spectrum  generally  obeys  a  nonlinear  rate 


equation  of  form 


£hm+2i 


NL 

Vyl'yl. 


jh  +  2S1* 


where  y_  *  YrPfx  . 
g  u 

Y  is  the  linear  growth  rate,  yl,  the  linear  damping  rate,  and  Y^the  nonlinear 


damping  rate.  The  fourth  term  represents  thermal  (noise)  emission  at 


Langnuir  frequencies.  Y^(k)  =  Yq  +  Y^(k)  consists  of  linear  collisional 
damping  and  (k-deperxlent)  Landau  damping.  For  k  below  about  0.2  kp, 
where  kp  is  the  Debye  wavenumber,  the  Landau  damping  is  negligible.  The 
quantity  P  is  a  measure  of  the  pump  field  strength  and  is  given  by 


O  -l 

l»unQ  lbYjVuip 
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where  n  is  the  electron  density  and  0,  the  corrmon  eleci  ron  and  ion 
temperature  in  energy  units.  P=1  defines  the  minimum  tl^eshold  field, 

^*0 1  min.  thresh. ,flDr  exci-fcation  parametric  instability,  The  P-value 
tells  us  how  much  the  pump  power  is  above  or  below  the  minimum  threshold. 

The  factor,  f  =  f(k,x)  is  called  the  frequency  mismatch  factor  and  is 
equal  to  one  for  k  equal  to  the  frequency-matched  wave  number,  k^,  defined  by 

wo‘  “t/V  +  uia(V*  ^ 

vhere  u)Q  is  the  pump  frequency,  oj^Ck)  is  the  Langnuir  frequency  (wp  +  3k2V^)lii 
and  <*>ia(k)  is  the  icn  acoustic  frequency,  ~  (nigA^)*5  kVg.  For  k<k^,  the 
mismatch  factor  f  is  <  1.  Note  is  fixed  by  the  density,  temperature, 

MT 

and  punp  frequency.  The  nonlinear  damping  rate  y‘  (k,x)  depends  on  I,  ,and 

has  been  the  subject  of  various  theories  and  approximations  to  be  described. 

Ihe  emission  term,  consists  of  linear  Cerenkdv  and/or  bremunstrahlung 

emission  at  the  langnuir  frequency  o^Ck)  plus  a  term,  first  derived  by  two  of 
13  4 

the  authors  *  ,  representing  the  beating  of  Cerenkov  emission  at  ion-acoustic 
frequencies  with  the  pump  field  to  produce  emission  around  the  Langmiir 
frequency.  It  dominates  the  usual  linear  emissions  even  for  P  <  1. 

Since  we  have  assumed  azimuthal  symmetry  about  the  punp  electric  field, 
the  solution  to  equation  1  for  some  appropriate  model  of  the  nonlinear 
damping  rate  gives  Ip(t,k,x)  for  a  homogeneous  plasma.  Asymptotic 
time-independent  solutions  L^CkjX)  obtained  by  setting  the  right  side  of 
equation  1  equal  to  zero  and  solving  the  resulting  nonlinear  equation  are 
called  saturated  steady-state  Intensity  spectra.  1^  and  are  necessary 
in  order  to  compute  scattering  cross-sections  or  anomalous  resistivities. 

In  all  of  the  saturation  theories  developed  for  the  equal  temperature 
decay  instability  the  ion  acoustic  spectrum  I ^  is  simply  proportional  to 
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an^  given  by  I2  =  ^Ygk  /kp]  1^.  In  the  ionosphere  the  intensity  of 
parametrically  excited  ion-acoustic  waves  turns  out  to  be  much  lower  than 
that  of  the  Langmir  waves .  Physically,  the  ion  acoustic  waves  are  heavily 
dairped  quasi-modes  which  are  "dragged”  along  by  the  strongly  enhanced 
Langnui?  modes. 

Another  quantity  of  interest  is  the  total  energy,  E,  in  Langrair  waves 
for  a  given  electron  density,  temperature  and  pump  frequency  (hence,  for 
a  given  k^) 

Cd\  L(k) 

E  fep  m 


Ihe  two  theories  for  y^  to  be  described  below  have  both  been  solved  in 
various  approximations  for  I^k),  E  and  for  the  high  frequency  electro¬ 
magnetic  anomalous  dissipative  conductivity,  o^.  An  expression  for  the 

x 

(transverse)  anomalous  dissipative  conductivity,  a^,  applicable  to 
either  theory  is  given  by. 


I-^k) 

w~ 


(5) 


Ihis  is  a  simple  expression  of  energy  conservation.  The  energy  loss  rate 

of  the  pump  (left  side)  is  balanced  by  the  total  energy  gain  (sunned  over 

k  on  the  right  side) of  Langnuir  waves,  due  to  their  linear  growth  at  the 
2 

rate  yg  =  x  fPyL,  indicated  in  the  rate  equation  1. 

b)  Saturation  of  Langnuir  waves  by  Induced  scattering  off  ions 
In  this  theory,  the  parametrically  unstable  Langruir  waves  are  scattered 
off  ions  out  of  the  unstable  region  of  phase  space  where  y  (k)  is  largest 
and  into  more  stable  .regions.  The  scattering  rate,  is  proportional 
to  the  intensity  of  Iangoiir  waves  and  the  process  is  therefore  said  to  be 
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"Induced”  or  "stimulated."  (The  scattered  waves  are  entirely  electrostatic 
In  this  process,  which  is  not  to  be  confused  with  electromagnetic  scattering.) 
The  earliest  theory  *  made  the  crude  approximation  that  1^  was  independent 
of  x  and  therefore  that  the  Langmuir  waves  were  isotropically  distributed  in 
angle  about  the  direction  of  the  pump  electric  field.  We  will  refer  to 

n 

this  as  the  isotnopized  theory.  In  the  published  form  of  this  work  1^00 
spread  continuously  from  to  lower  k  as  P  was  increased,  and  the  total 
energy,  E,  was  found  to  increase  as  P  .  An  estimate  of  the  radar  backscatter 
cross-section  was  also  made  at  this  time  and  appeared  to  give  too  large  a 


result  in  comparison  with  the  Arecibo  experiments.  Later,  more  refined 
computations^ * ^ ^ ,  still  based,  however,  on  resonant  approximations  for 
dielectric  functions  (and  for  f),  and  on  the  so-called  derivative  approximations 
for  Yj^  led  to  311  angular  distribution  of  Langmuir  waves,  I^(x)  which  was 
highly  peaked  along  the  pump  electric  field  (i.e.,  at  x  =  1).  For  the  iono¬ 


sphere,  this  means  peaking  along  the  geomagnetic  field  direction.  The  width 

f+1  to 

in  k  of  the  angle-averaged  anisotropic  spectrum,  L^k)  ^  I^kjX) 

was  found  to  be  0.67  (rr^/m^  )  t .  The  anomalous  dissipative  conductivity 
was  found  from  the  energy  balance  equation  (5)  by  further  noting  that 
the  induced  scattering  process  represented  in  this  saturation  theory  by 
merely  redistributes  the  energy  among  Langnuir  waves  of  different  k 


without  acting  as  an  energy  sink  for  the  entire  collection  of  waves  integrated 


over  k.  By  integrating  the  rate  equation  (1)  over  k  it  was  therefore 
possible  to  show  that  the  right  side  of  equation  (5)  was  still  correct  if 
the  growth  rate  yg  were  replaced  by  the  damping  rate  yl>  provided  the 
emission  term  S.^  was  negligible  (valid  for  P  >  1).  Physically,  this  says 
that  the  only  real  energy  sink  in  this  theory  is  the  linear  damping  rate  of 
plasma  waves .  The  energy  flows  from  the  pump  to  the  plasma  waves,  is  then 
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redistributed  In  k- space  by  the  nonlinear  saturation  mechanism,  and  is 

finally  delivered  to  the  plasma  particles  by  the  collisional  and  Landau 

damping  present  in  y For  in  the  collision-damped  region,  therefore, 

NL 

the  anomalous  dissipative  conductivity  is  given  by 

°tl”  (6) 

Ihe  results  of  Reference  6  and  of  other  and  more  refined  theories  of  this 

saturation  all  give  slightly  different  numerical  factors  for  E,  but  all 

yield  ao^in  terms  of  E,  through  equation  (6).  The  results  of  these 

theories  can  be  sunmarized  by  using  a  multiplicative  parameter  n,  such  that 

E  *=  iJn  P2n0  (7) 

r 

The  values  of  n  corresponding  to  each  theory  are  given  in  Thble  1.  Some 

of  these  theories  make  more  accurate  predictions  for  the  structure  of 

Lj^k.x).  The  time-dependent  solutions16*17  of  equation  1  provide  estimates 

of  the  time  scale  necessary  for  the  establisiiment  of  the  steady  state 

spectrum  in  the  ionosphere.  In  the  case  of  rapid  pump  turn-on  a  time  of 

15  m  sec  is  required  to  establish  the  steady  state  spectrum  if  P  =  5,  while 

5  m  sec  is  required  if  P  =  20.  The  spectrum  may  develop  adiabatically  with 

the  pump,  if  the  latter  is  turned  on  slowly  enough.  The  rise  time  must 

be  slower  than  0.75  seconds  for  P  =  5,  or  3  seconds  for  P  *  20.  The  most 

accurate  of  the  1-D  time-dependent17  solutions  yields  a  discrete  structure 

in  k,  consisting  of  peaks  at  values  of  k  corresponding  to  langnuir  waves 

shifted  in  frequency  by  multiples  of  the  ion  acoustic  frequency.  This 

18 

structure  was  also  found  independently  by  Kuo  and  Fejer.  The  most 
complete  refinement  of  our  theory  to  date  appears  to  be  tne  numerical 
integration  by  Kuo  and  Fejer  of  the  rate  equation  (1),  in  time,  wave 
number,  and  angle,  using  the  exact  dielectric  function  and  no  derivative 
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approximation  for  y^*  The  most  interesting  feature  of  their  solution  is 
angular  spreading  v/hich  increases  with  P  and  is  out  to  about  20°  at  P  =  2. 

Th;  earlier  k-structure  is  again  found,  leading  to  umbrella-shaped  discs 
ot  Langnuir-wave  intensity  in  k-space.  No  values  for  the  total  energy  E 
are  given  in  this  work,  but  our  earlier  isotropized  model  leads  us  to  expect 
an  n-value  of  between  1  and  2. 

Induced  scattering  of  Langmuir  waves  off  ions  appears  to  be  the  dominant 
saturation  mechanism  for  matched  wavenumbers ,  <  .23'^,  in  the  collision- 

damped  region,  and  the  orbit-perturbation  process  to  be  discussed  in  the 
next  section  appears  to  be  the  dominant  saturation  for  matched  wave  numbers 
in  the  landau-damped  region,  k^  >  .231^.  In  Appendix  C,  in  addition  to 
treating  effects  of  spatial  inhomogeneity  on  the  nonlinear  pump  propagation, 
we  have  extended  the  induced-scattering  saturation  theory  of  c^scmewhat , 
to  include  Landau-damping  (qualitatively),  four-mode  effects  (which  raise 
the  threshold  when  wia  <  yl),  and  broad  spectrum  effects  (at  high  pump 
powers).  Refer  to  Appendix  C  for  details. 

c)  Saturation  of  Langmuir  waves  by  orbit-perturbation 

Tills  saturation  mechanism  involves  absorption  of  the  parametrically 
excited  Langmuir  waves  by  electrons  whose  orbits  have  been  nonlinearly 
perturbed  by  the  waves.  Sometimes  it  is  called  resonance-broadened  Landau¬ 
damping.  The  theory  of  the  nonlinear  damping  rate,  and  the  spectrum 
I^k)  is  described  in  detail  in  Appendix  A  and  in  Reference  8.  Energy 

comparisons  have  shown  that  for  frequency  matched  wave-numbers  k  in  the 

m 

Landau-damped  region  k^  >  .23k^  orbit-perturbation  is  the  dominant  saturation 
mechanism.  For  smaller  k^,  landau-damping  is  exponentially  weak  since  there 
are  so  few  electrons  whose  velocities  match  the  wave  phase  velocity.  The 
resonance  broadening  effect  for  <  .23^  (at  least  for  the  pump  intensities 
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of  Interest  In  the  ionosphere)  can’t  reach  back  far  enough  into  the  velocity 
distribution  to  find  an  appreciable  number  of  electrons.  For  these  wave 
numbers,  the  induced-scattering  saturation  described  in  the  previous  sec¬ 
tion  dominates.  Physically,  the  distinction  between  the  two  theories  is 
made  by  asking  whether  or  not  a  wave  is  scattered  out  of  the  unstable  region 
(large  y  )  of  phase  space  by  the  induced-scattering  mechanism  before  it  is 
locally  absorbed  by  resonance  broadened  Landau-damping.  Thus,  a  comparison 
of  the  rates  for  the  two  processes  should  be  made.  A  careful  treatment  of 
the  region  k^  ~  0.2k^,  where  the  two  rates  are  apparently  comparable  has  not 
been  made  to  date.  As  we  shall  see  in  the  next  section,  plasma  waves  with 

km  7311165  111  both  nonlinear  saturation  regimes  are  always  excited,  although 
at  different  altitudes. 

The  validity  of  orbit  perturbation  theory  requires  that  the  corre¬ 
lation  time  of  the  fluctuating  electric  field  be  short  compared  to  an 
electron  trapping  bounce  time.  This  is  guaranteed  if  the  total  wave  energy 
density  is  less  than  the  particle  energy  density,  a  condition  which  is 
satisfied  for  ionosphere  parameters. 

Hie  spectrum  obtained  is  broad  in  angle,  provided  Px2  >  1,  and  narrow 
in  k.  In  the  Landau-damped  region  of  the  ionosphere  the  spectrum  is  con¬ 
fined  to  a  half  width,  Ak  =  3  x  10  about  the  frequency-matched  wave 
number,  k^.  The  angular  part  of  the  spectral  intensity  is  given  by 


Ux)  =  V^A  (Lx)  •  time  x  3.1) ^  to  (1  +  r)Px2  -  r) 

1  J(*rf  1  ^  (  (kp/k/  7  i).5  )  *  W 

where  r  is  the  ratio  of  the  collision  damping  rate  y  to  the  Landau- 

c 

damping  rate  evaluated  at  k^.  This  expression  is  valid  provided  Px2  >.  1. 

We  see  from  eq-  (8)  that  the  spectrum  becomes  increasingly  bread  in  angle. 
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as  the  pump  is  increased,  finally  reading  isotropy.  Figure  1  is  a  sketch 
of  I^(x)  as  a  flinction  of  0  =  arc  cosine  x.  The  maximum  spreading  angle  is 

2  _i 

given  by  x^^  =  P  ,  0^^  =  arc  cosine  x^^.  The  total  saturated  wave  energy 
density  is  given  by 


E  6.75 

**  [<W!  + 


[find  +  r)p]2  -  i»[ln(l  +  r)P  -  2]  - 


8 


/TTFp 


(9) 

For  km  in  the  Landau-damped  region,  so  that  r  «  1,  the  total  wave  energy 
density  E  is  seen  to  increase  rather  slowly  with  the  pump  ratio  P  for  P 
large.  In  Figure  2  we  illustrate  E  as  a  function  of  P  for  *  0. 
r  *  1.  Fran  the  curve  one  observes  that  E  increases  approximately  as  P. 
Equation  (5)  may  again  be  employed  to  give  the  nonlinear  high-frequency 
dissipative  conductivity,  A  direct  application  of  (5)  yields 


orbit  perturbation 


a  Y  P  E 
3.  L%2 


(10) 


a  is  a  numerical  factor  resulting  from  the  angular  dependence  of  I(x).  It 
varies  slowly  between  1  and  2  as  a  function  of  P.  E  is  given  by  equation 
(9),  and  by  Figure  2. 

Whereas  for  the  induced-scattering  saturation  mechanism  the  energy 
flowed  from  the  punp  to  the  waves  where  it  was  redistributed  in  k-space 
and  eventually  delivered  to  the  particles  via  linear  damping,  the  perturbed- 
°r>  T  t  saturation  mechanism  predicts  the  energy  to  flow  from  the  punp  to  the 
waves,  and  finally  to  the  particles  via  the  nonlinear  process  of  resonance- 
broadened  Landau-damping.  Finally,  it  is  interesting  to  note  that  (10) 
has  met  with  considerable  success  in  comparisons  with  numerical  simulations. 
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Kmer  and  Dawson*”  have  determined  fran  numerical  simulations  an  effective 
electron-ion  collision  rate,  v  =  of  precisely  the  same  form  as 

equation  (10).  Further  details  of  perturbed-orbit  saturation  theory  may 
be  found  in  Appendix  A. 

III.  Effect  of  density  inhcmogenelty  in  the  ionosphere;  applications  of 
nonlinear  saturation  theory  to  absorption,  scattering  and  airglow. 

a)  Anomalous  absorption  and  the  modifier  field  pattern.  We  have  taken 
into  account  sane  important  effects  of  density  inhomogeneity  on  parametric 
instabilities  in  the  ionosphere  through  the  following  single  model  :  The 
homogeneous  theory  of  the  previous  section  is  assumed  to  hold  at  each  altitude, 
provided  the  local  values  of  the  plasma  frequency  and  the  purrp  (modifier) 
electric  field  airplitude  at  that  altitude  are  used.  The  pump  electric  field 
is  not  assured  gi'.ven,  but  is  found  locally  in  an  assumed  vertically  plane¬ 
layered  ionosphere,  in  terms  of  the  incident  vacuum  intensity.  The  self- 
consistent  pump  field  solution  includes  the  effect  of  nonlinear  pump  depletion 
caused  by  the  flow  of  power  from  the  punp  to  the  parametrically  excited  plasma 

waves  at  each  altitude.  This  effect  is  contained  in  the  dissipative  nonlinear 
NL 

conductivity,  o^,  in  the  wave  equation  for  the  pump.  Effects  of  inhomogeneity 

on  the  linear  parametric  instability  theory  should  not  be  important  for  the 

21 

ionosphere  because  the  scale  height  exceeds  the  electron  mean  free  path  . 

We  are  also  ignoring  effects  of  inhomogeneity  in  the  saturation  theory  and  on 

the  plasma  wave  propagation.  (Sane  preliminary  estimates  of  inhomogeneous 

22 

plasma  wave  propagation  effects  have  been  carried  out  by  Amush  ,  neglecting 
saturation,  and  by  Perkins  ,  with  saturation  included.)  Basically,  our 
model  is  one  of  nonlinear  ordinary  wave  punp  propagation  vertically  upwards 
through  a  linear  density-altitude  profile  with  a  horizontal  magnetic  field 
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(or  no  magnetic  field).  Only  after  the  pimp  field  is  found  self-consistently 
in  this  manner  do  we  know  the  plasma  wave  intensity  spectrum  as  a  function  of 
altitude.  As  in  the  linear  solution,  the  local  pump  wavelength  is  much  larger 
than  the  wavelength  of  the  excited  plasma  waves  —  one  justification  for 
using  locally  homogeneous  theory  for  the  latter.  The  frequency  matched 
wave-number  defined  in  equation  (3)  is  now  a  function  of  altitude,  since 
Up  =  4rme  /hi  is  a  function  of  altitude.  Ignoring  u^a  relative  to  up,  the 
matched  wave  number  is 


km  1  Az 

=rjiT’ 


(id 


where  H  is  the  scale  height  and  Az  is  distance  below  the  reflection  altitude. 
Thus,  Langjnuir  waves  of  different  wavenumbers  are  excited  at  different 
altitudes;  long  wavelengths  are  excited  above  short  wavelengths.  Since 
the  shorter  plasma  waves  are  Landau-damped  and  the  longer  ones  collisionally 
damped,  the  induced-scattering  saturation  mechanism  should  be  valid  at 
the  higher  altitudes  and  the  orbit-perturbation  mechanism  at  the  lower 
altitudes.  The  situation  is  illustrated  in  Figure  3. 

r 

<e  have  solved  the  following  nonlinear  inhomogeneous  wave  equation  for 
the  modifier,  including  linear  damping  and  the  nonlinear  damping  represented 
by  of,  defined  in  equations  6  and  7: 


The  geometry  is  shown  in  Figure  4.  z  is  measured  in  units  of  c/w  .  z  =  0 

o 

is  the  reflection  point  and  z  =  zQ  is  the  point  of  incidence. 

The  first  set  of  computer  runs  have  been  performed  using  the  following 
set  of  parameters: 
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«5  ~A 

n  ■  5  x  l(r  cm 
o 

0  -  0.2eV  -  23L9°K  (daytime) 

f 

a  i  (ly^)1^3-  .00583  (oxygen) 

u  /2ir  =  modifier  (pump)  frequency  -  6.35 
o 

Derived  quantities  based  on  the  above  parameters  aie: 
(0  *  4  x  10*^  radians/sec 


ei. 


2yc0  -  350  Hz 


th 


■  °th.—  =  16  — 
4im  u 

o  e  p 


7  x  10 


,-5 


(13) 


I  =  E2  =  168  yW/m2  (locally)  (l1*) 

th  8?  oth 

The  peak  plasma  frequency  is  equal  to  the  pump  frequency ,  ty 

vei0  is  the  electron-ion  collision  frequency  at  the  reflection  altitude. 

'Die  threshold  for  the  parametric  instability  is  indicated  as  Ith,  and 

cE- 

should  be  understood  as  the  value  which  (including  swelling  due  to 

inhomogeneity)  must  exceed  locally  to  excite  the  parametric  decay  instability. 

wo 

We  consider  two  scale  heights:  H  ■  37-5  km  (zQ  =-£-  H  ■  5000)  and  H  - 
uH  h 

75  km  (z0  =  *  10  ). 

The  integrated  total  absorption  is  numerically  measured  by  the  quantity 
1  -  R,  where  R  ■  jE^j/E^J2  is  the  power  reflection  coefficient  at  the 
"source  point"  z  *  zQ.  We  also  compute  for  each  case  the  integrated  purely 
linear  reflectivity  R]_in,  obtained  by  Integrating  the  wave  equation  using 
only  the  linear  dissipative  conductivity.  These  results  conpare  favorably 
(within  5055)  with  the  value  of  R^  calculated  from  geometric  optics. 
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The  quantity  FL^  -  R  =  (1  -  R)  -  (1  -  R^)  is  a  measure  of  the  anomalous 
absorption  (excess  beyond  linear). 

In  Figure  5  both  ]  -  R  (left  axis)  and  R . j, n~R  (right  axis)  are  shewn  for 

2 

an  Incident  vacuum  Intensity  of  50  yW/m  ,  which  is  about  the  maximum 
presently  available  at  the  Plattevllle  transmitter^  in  Colorado.  (Approxi¬ 
mately  one  half  this  value  is  attainable  at  Arecibo,  Puerto  Rico.  The 
absorption  .s  plotted  as  a  function  of  the  multiplicative  scaling  parameter, 
n,  defined  in  equation  (7)  and  T&ble  1.  As  explained  in  Table  1,  the 
results  of  the  best  available  theories  indicate  n-values  between  h  and  3/^* 

For  a  scale  height  of  37*5  km,  the  linear  absorption  1  -  R  =  S%  for 
the  above  parameters  of  (13)  and  (14 ) .  The  excess  nonlinear  (anomalous) 
absorption  R]1n-R  is  between  9%  and  2355  depending  on  n.  The  behavior  of 
R|1n-R  with  n  is  roughly  linear. 

a  scale  height  of  75  km,  the  linear  absorption  is  1  -  R^n  =  17% 
for  the  same  parameters .  The  excess  nonlinear  absorption  -  R  Ls 
between  and  33%,  depending  on  n.  Longer  density  scale  lengths  are 
therefore  seen  to  lead  to  stronger  absorption.  The  physical  reason  Is 
that  the  absorption  region  is  "stretched,"  although  from  the  above  results 
it  appears  that  this  is  a  stronger  effect  for  the  linear  absorption  than 
the  nonlinear  contribution  (presumably  because  the  nonlinear  absorption 
region  ls  thinner  spatially). 

Tbs  above  results  are  for  relatively  high  daytime  electron  temperatures. 
Strictly  speaking,  the  effects  of  photoelectrons  on  the  plasma  wave  damping 
should  be  included  for  daytime  parameters.  We  have  crudely  estimated  the 
effects  of  photoelectrons  in  a  separate  data  run  by  adjusting  the  Landau 
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damping  to  be  three  times  larger  than  the  eollisional  damping  of  the  Langnuir 

w?ves  at  a  matched  wave  number  of  2  O.lk^.  (Ordinarily  the  Landau 

damping  is  much  smaller  than  the  collision  damping. )  The  factor  3  was 

chosen  as  an  average.  At  this  density  and  temperature  downward  propagating 

plasma  waves  have  a  photoelectron  damping  rate  equal  to  the  eollisional 

damping  rate  and  upward  propagating  plasma  waves  have  photoelectron  dairping 

five  times  eollisional  dairping.  (The  photoelectrons  originate  mainly  at 

lower  altitudes.)  Also,  there  should  not  be  many  photoelectrons  at  energies 

beyond  jK -j“)  ~  20eV.  This  has  the  effect  of  raising  the  threshold  for 

parametric  instability  by  a  factor  of  -4  at  altitudes  in  the  vicinity  of 

z/zQ  =  .03  and  by  more  below.  It  does  not  alter  the  nonlinear  conductivity 

at  higher  altitudes  because  the  phase  velocities  at  the  correspondingly 

smaller  matched  wave  numbers  are  too  high  to  allow  wave  resonance  with 

photoelectrons.  The  net  effect  is  to  condense  the  spatial  region  of  nonlinear 

absorption.  For  all  parameters  the  same  as  in  Figure  5,  but  photoelectrons 

included  in  the  above  manner,  and  a  scale  height  of  75  km,  we  find  for  n  * 

2/3,  Rj1n  -  R  =  11?  anomalous  absorption  compared  to  14?  anomalous  absorption 

without  photoelectrons;  and  for  n  =  2,  -  R  =  27?  anomalous  absorption, 

compared  to  33?  anomalous  absorption,  without  photoelectrons.  Photoelectrons 

thus  appear  to  reduce  the  percent  of  anomalous  absorption  by  a  factor  of 

about  20?.  When  photoelectron  damping  of  langnuir  waves  is  a  strong  effect, 

or,  indeed  even  if  the  normal  Landau  damping  region  is  to  be  treated 

properly,  the  absorptive  conductivity  arising  from  orbit  perturbation 
8  A 

saturation  *  must  be  used  rather  than  the  conductivity  associated  with 
weak  turbulence  saturation,  which  we  have  been  using  here.  We  are  presently 
putting  into  effect  a  program  for  doing  this. 

We  have  also  made  some  runs  at  lower  electron  temperature  (,13eV  = 

! 

{ 
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1500°K).  Since  the  linear  absorptive  conductivity  goes  as  \>ei  and  hence 
as  0~^,  the  linear  absorption  should  and  did  increase.  The  nonlinear 
absorption  also  increased,  although  not  as  strongly  as  the  linear  absorption. 
Making  some  numerical  comparisons  for  the  H  =  75  km  shown  in  Figure  5, 
but  with  all  other  parameter’s  except  temperature  the  same,  we  find  the 
linear  absorption  1  -  increased  from  172  to  282  as  0e  is  decreased  frcm 
2320°K  to  150 0°K,  whereas  the  anomalous  absorption  -  R  increased  from 
142  (.2eV)  to  192  (.13eV)  for  n=2/3  and  from  33%  (.2eV)to39%  (.13eV)  for  17=2. 
Nighttime  anomalous  absorption  would  thus  appear  to  be  stronger  than  day¬ 
time  anomalous  absorption,  although  orbit  perturbation  effects  must  be 
included  to  really  determine  whether  this  is  true.  Daytime  neutral  col- 
lisional  absor  tion  of  the  modifier  (pump)  at  lower  altitudes  (D-region) 
would  make  daytime  anomalous  absorption  still  weaker. 

In  Figure  6  we  have  plotted  the  absorption  1  -  R  as  a  function  of 
incident  intensity  for  the  scale  height  H  =  75  km  using  the  same  n-values 
aid  other,  fixed  parameters  as  in  Figure  5.  The  absorption  goes  up 
womewhat  more  slowly  than  linear. 

In  Figure  7  we  plot  the  pump  radiation  electric  field  amplitude  squarec’ 

versus  altitude.  The  pump  amplitude  is  normalized  to  the  incident  amplitude 

and  the  altitude  is  in  units  of  c/uQ.  For  the  parameters  of  (13)  and  (14) 

this  distance  is  7.5  meters.  The  spatial  wave  patterns  are  plotted  for  an 

incident  vacuum  intensity  of  50  uW/m  ,  a  scale  heigh, t  of  75  km  (zq  =  10  ), 

5  -3 

ard  a  nonlinear  scaling  parameter  n  =  2.  The  density  is  5  x  10  cm  and 
the  temperature  is  0.2eV.  Results  are  shown  both  for  a  purely  linear 
absorptive  conductivity  and  for  the  linear  plus  nonlinear  absorptive  con¬ 
ductivity.  The  peak  of  the  Airy  function  swelling  pattern  is  indicated  by 
X  for  comparison.  The  effect  of  pump  depletion  is  evident.  The  peak  swelling 
factor  at  the  first  maximum  is  reduced  from  82  with  no  absorption  to  71 
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with  linear  absorption,  to  56  with  linear  plus  nonlinear  absorption.  Ihe 
amplitudes  at  all  the  maxima  are  reduced  by  somewhat  smaller  percentages. 

For  the  case  shown  we  are  not  very  far  above  threshold  even  at  the  first 
maxima  of  the  field  pattern.  Based  on  the  actual  swelling  at  the  first 
maximum  of  the  curve  which  includes  linear  plus  nonlinear  absorption  we 
find  P  =  17.  If  we  had  used  the  Airy  function  swelling  factor  this  number 
would  have  come  out  to  be  P  -  24. 

In  Figure  8  we  have  plotted  the  negative  of  the  total  Poynting  flux 

In  units  of  the  incident  flux  for  the  same  parameters  as  in  Figure  In 

che  absence  of  any  absorption  this  net  flux  must  vanish.  Ihe  lower  curve 

in  Figure  8  includes  only  linear  absorption.  The  upper  curve  includes  both 

linear  and  nonlinear  absorption.  The  steeper  regions  of  the  oscillating 

flux  are  the  regions  where  the  absorption  is  strongest,  and  correspond  to  the 

maxima  of  Eq(z) .  Ihe  curve  including  both  linear  and  nonlinear  absorption 

shows  stronger  oscillations  than  the  curve  Including  only  linear  absorption 
wr  ?  4 

because  9J/3z  «  o  Eq  «  Eq.  Also,  the  curve  including  both  linear  and  non¬ 
linear  absorption  rises  to  a  higher  value  at  z  =  zQ  than  the  curve  with 
linear  absorption  alone  because  the  smaller  the  reflected  wave  at  that  point, 
the  higher  the  net  Poynting  flux. 

b)  UHF  scattering  near  the  plasma  frequency.  Ihe  peak  intensity  of  the 

p 

HF  modifier  probe  at  Arecibo  is  25  uW/m  .  Ihe  plasma  waves  probed  by  the 

430  Mhz  incident  wave  in  backseat tering  experiments  are  propagating  close 

to  the  vertical  and  have  wavenumber  k  =  4rr  x  430  Mhz/c  ■  0.18  cm-1.  Ihe 

pump  electric  field  is  assumed  to  lie  along  the  geomagnetic  field  direction, 

so  that  the  vertical  is  some  42°  from  the  puitp  direction  and  x2  =  cos2  42°  ■ 

0.55.  kp  depends  on  the  reflection  density  (i.e.,  on  f  )  and  temperature. 

It  ranges  from  k_.  =  3.4  cm”1  for  f  =  6  Mhz,  0  »  2400°K  to 
u  o  e 
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k  =  2  cm-1  f or  f  =  8.2  Mhz,  0^  =  1500°K.  The  second  value  is  for  the  set 
u  o  e 

of  parameters  cited  by  Kantor  ^ .  Thus,  k/k^  ranges  from  .05  to  .09 
depending  on  time  of  day  or  night.  The  frequency-matched  plasma  waves  are 
therefore  in  the  collision-damping  rather  than  the  Landau-damping  regime, 
so  the  saturation  theory  appropriate  to  cross-section  calculations  above 
threshold  is  the  one  based  on  induced -scattering  of  Langmuir  waves  off 
ions  rather  than  perturbed-orbit  theory.  At  these  wavenumbers ,  phctoelectrons 
cannot  play  a  very  important  role  in  the  Landau  damping  of  even  the  upward 
propagating  Langmuir  waves,  because  the  associated  phase  velocity  would 
Correspond  to  photoelectron  energies  of  at  least  25eV,  where  the  photo- 
electron  flux  rapidly  decreases.  From  equation  (11)  we  see  that  the  altitude 

_3 

from  which  the  scattering  return  is  expected  is  given  by  Az/H  =  7*5  x  10 
to  2.4  x  10"2,  or  from  about  0.5  to  2  kilometers  below  the  reflection  height. 

The  true  swelling  factor  to  these  altitudes  ,  for  a  modifier  with  15  uW/mc 
incident  intensity,  based  on  the  solutions  described  in  the  previous  section, 
is  -  40-60  provided  the  altitude  happens  to  correspond  to  a  maximum  in  the 
midifier  field.  In  general,  this  is  unlikely.  Kantor  has  verified  the 
standing  wave  pattern  by  allowing  the  density  profile  to  sink  without  track¬ 
ing  the  modifier  frequency  to  it.  For  plasma  waves  propagating  in  the 
pump  direction  at  this  altitude,  the  pump  peaks  would  be  between  2  and  8  times 
threshold  (P  =  2  to  8).  This  value  is  sensitive  to  f  ,  since  the  threshold 
power  goes  as  nQ^2.  P  =  2  corresponds  to  an  8  Mhz  pump;  P  -  8  corresponds 
to  a  6  Mhz  pump.  These  values  are  based  on  an  (optimistic)  75  l<m  scale 
height.  Even  the  most  detailed  theory  predicts  a  less  than  20  angular 

spread  of  plasma  waves  about  the  pump  direction  for  P  =  2,  and  effects  of 
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inhomogeneity  suggested  by  Perkins  would  probably  also  be  insufficient  to 
spread  the  saturated  spectrum  to  42°  with  aiy  appreciable  intensity. 
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Since  x2  =  0.55,  the  Langnuir  waves  probed  by  the  scatter  experiment  are 

probably  below  the  threshold  for  parametric  instability  for  this  case.  On 

the  other  hard  for  lower  HP  frequencies  if  one  happens  to  catch  a  pump  field 

18 

maxinum  it  is  possible  that  the  saturated  spectrum  has  spread  as  far  as 
U2°,  in  which  case  a  far  larger  cross-section  would  be  observed.  This  may 
account  for  some  of  the  variability  of  the  scattering  results. 

Ihe  detailed  comparisons  by  Perkins2^,  Harker2^,  and  Meltz  and  Tcmljano- 
vich2^  of  Arecibo  scattering  measurements2*^  with  the  anisotropic  and  iso¬ 
tropic  nonlinear  saturation  theories  of  the  scattering  cross-section  led  to 
the  conclusion  that  the  measured  waves  were  usually  stable.  Ihe  magnitude  of 
the  Arecibo  scattering  measurements  at  42°  to  the  punp  can  then  be  understood ty 
using  the  anisotropic  angular  spectrum  of  EuBois  and  Goldman6  and  Valeo, 
et  ai.9,  at  low  P  values.  At  this  angle  the  only  nonlinearity  entering  the 
cross-section  is  the  pump-enhanced  ion-acoustic  Cerenkov  emission  term 
first  predicted  in  1967  by  DuBois  aid  Goldman4b.  Ihis  arises  from  the  beat¬ 
ing  of  the  modifier  with  the  ion-acoustic  frequency  Cerenkov  emission  to 
produce  emission  at  the  Langmuir  wave  frequency.  Ihis  emission  is  larger 
by  a  factor  Pau/ui^  than  the  direct  Langnuir  bremsstrahlung  emission.  Ihe 
scattering  model  for  low  P  at  0  =  42°  is  called  critical  fluctuation  theory. 

It  predicts  that  the  cross-section  increases  linearly  with  P  at  the 
Langmuir  frequency.  For  large  P,  or  for  9  close  to  0  for  smaller  P,  the 
cross-section  should  increase  as  P2  times  a  large  coefficient.  A  dramatic 
increase  (by  as  much  as  three  orders  of  magnitude  for  an  increase  of  P  by 
5)  in  the  cross-section  should  be  characteristic  of  this  regime. 

Ihe  considerable  detail  of  the  Arecibo  backscatter  spectrum  presents 
many  spectral  features  besides  the  plasma  peak  to  be  explained  by  theory. 

Ihe  critical  fluctuation  theory  predicts  a  ratio  cf  peak  intensity  of  the 


enhanced  Stokes  Inngjnuir  line  (downshifted  from  the  pump  frequency  by  the 
ion-acoustic  frequency)  to  the  ion-acoustic  line >  which  is  in  agreement  with 
experiment.*  Likewise,  the  ratio  of  the  anti-Stokes  line  (Langnulr  line 
upshifted  from  the  pump  frequency  by  the  ion-acoustic  line)  to  the  Stokes 
line  is  in  agreement  with  experiment.^  A  spectral  peak  shifted  by  exactly 
the  pump  frequency  has  been  identified  as  the  high  frequency  component  of 
the  purely-growing  instability,  which  is  a  different  branch  of  the  dispersion 
relation.  In  Appendix  B,  the  critical  fluctuation  intensity  is  computed  for 
the  purely  growing  instability.  Assuming  values  of  P  which  appear  to  explain 
the  features  of  the  parametric  decay  instability,  we  compute  the  ratio  of 
the  purely-growing  line  to  the  decay  line  (Stokes  line) .  This  ratio  of 
cross-sections  turns  out  to  be  larger  by  a  factor  of  2  or  3  than  the  ratios 
observed  at  Arecibo.  Ibis  factor  could  arise  from  the  assumption  of  equal 
electron  and  ion  temperatures .  In  addition  to  Indicating  that  the  decay 
line  peak  and  the  growing  line  peak  should  be  of  comparable  order  of  magni¬ 
tude,  as  observed,  the  theory  predicts  a  variability  in  the  ratio  with 
heater  frequency  and  scale  height ,  since  the  observed  purely  growing  node  is 
excited  at  a  lower  altitude  than  that  of  the  decay  instability.  One  feature 
which  has  still  not  been  explained  consistently  with  the  critical-fluctuation 
theory  is  the  structure  of  the  so-called  broad  bump  in  the  incoherent 
scatter  data  near  the  stokes  frequency.  The  shape  is  qualitatively  similar 
to  that  predicted  in  detailed  spectral  shape  calculations'^*^;  however, 
these  theories  require  that  the  pump  te  well  above  threshold.  Ihe  scattering 
models  constricted  by  the  Boulder  group  predict,  in  addition  to  the  below 
and  above  threshold  behavior  indicated  above,  a  third  dramatically  different 
regime,  which  has  not  yet  been  probed  experimentally.  As  explained  above, 
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the  values  of  k^/k^  probed  by  the  4 30  Mhz  scatter  are  .05  to  .09.  A  factor 
of  3  or  4  Increase  In  the  incident  scattering  wave  frequency  would  bring 
the  frequency-matched  wave  number  into  the  Landau-danping  regime.  Thus,  the 
use  of  an  L-band  (1.3  GKz)  probe  should  give  cross-sections  characteristic  of  orbit 
perturbation  saturation.  The  detailed  predictions  here  are  given  in  Appendix 
A,  although  we  note  from  Figure  1  that  the  spectrum  will  be  broad  in  angle 
even  for  P  of  order  2,  so  large  scattering  cross-sections  may  be  expected. 

c)  Acceleration  of  electrons  and  alrglow. 

The  large  enhancements  of  6300  8  airglow  over  Boulde^  and  Arecibo  suggest 
that  the  electron  velocity  distribution  function  is  significantly  altered 
during  ionospheric  modification.  The  alteration  of  the  velocity  distribution 
can  have  important  consequences:  for  example,  the  formation  of  energetic 
tails  may  cause  plasma  wave  instability  or  change  the  saturation  results  of 
orbit  perturbation  theory.  In  Appendix  D,  a  theory  for  quantitatively 
predicting  the  change  in  the  electron  velocity  distribution  has  been  developed. 
Such  changes  are  primarily  due  to  acceleration  of  electrons  by  short  wave¬ 
length  (~  15  cm)  Langmuir  waves  propagating  at  Landau-damped  altitudes.  The 
electrons  mainly  affected  are  those  whose  velocities  are  close  to  the  wave 
phase  velocity.  The  altitudes  from  which  the  accelerated  electrons  originate 
are  illustrated  in  Figure  3-  The  airglow  measurements  provide  a  direct 
measurement  of  the  enhanced  photon  flux  from  atomic  excitations.  The  photon 
flux  is,  in  turn,  a  measure  of  the  enhanced  flux  of  electrons  responsible 

li 

for  the  excitation.  In  Appendix  D  it  is  found  that  5  x  10  electrons  per 
c.c.  are  produced  per  second  in  the  proper  velocity  range  to  excite  6300  8 

airglow,  when  P  *  3,  and  k^/kp  ~  0.2.  The  total  wave  energy  for  these  parameters 

-3  9 

is  10  nO,  according  to  perturbed-orbit  theory.  From  this  a  flux  of  10 

2 

electrons  per  cm  per  sec.  results.  This  flux  appears  sufficient  to  explain 
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the  observed  50  Rayleigh  enhancements  of  the  6 300  8  airglov;  over  Boulder. 
The  acceleration  mechanism  is  made  possible  by  the  resonance  broadening 
predicted  by  orbit-perturbation  theory,  which  accelerates  electrons  which 
are  as  much  as  a  thermal  velocity  away  from  the  wave  phase  matched  velocity 
(at  5  thermal  velocities). 


IV.  Effects  of  the  geomagnetic  field;  stimulated  diffusion  scattering 

a)  The  diffusion  mode 

In  addition  to  supporting  wave-like  responses  to  density  perturbations, 
a  homogeneous  plasma  can  respond  with  a  diffusion-relaxation  of  the  pertur¬ 
bation.  As  an  example,  electron-ion  collisions  can  cause  a  random  walk  of 
particles  across  the  magnetic  field.  This  leads  to  a  diffusion  mode  with 
wave  vector  very  close  to  orthogonal  to  the  magnetic  field.  The  situation 
is  illustrated  in  Figure  10.  Diffusion  along  the  magnetic  field  lines  occurs 
much  more  rapidly  than  diffusion  across  the  field,  provided  the  electron 
cyclotron  frequency  is  much  larger  than  the  electron-ion  collision  fre¬ 
quency,  a  condition  which  is  satisfied  .in  the  ionosphere.  The  cross-field 
diffusion  coefficient  in  the  simplest  approximation  is 
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where  v  .  is  the  electron-ion  collision  frequency  and  R  is  the  electron 
el  e 

gyro-radius.  The  linear  fluctuation-dissipation  theorem  predicts  an 
equilibrium  electron  density  fluctuation  spectrum  associated  with  this 
diffusion  node,  having  k  |  B,  and  a  bell-shaped  frequency  structure  around 
zero  with  half -width  given  by 

Au  =  D.k^  .  (15) 


The  integrated  spectrum  in  equilibrium  is  independent  of  vgi  or  B  and  is  of 
the  same  order  of  magnitude  as  the  integrated  low  frequency  spectrum  in  the 
absence  of  B.  The  high  frequency  scattering  cross-section  is  directly  pro¬ 
portional  to  the  electron  density  fluctuation  spectrum.  Baron  and  Petrie eks‘ 
may  have  observed  this  structure  in  L-band  scattering  experiments  orthogonal 
to  the  geomagnetic  field. 
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For  the  F-layer,  with  v  assumed  to  be  *100  Hz  and  R  =  2  cm,  D,  Is  of 
2  e  1 
the  order  3000  cm  /sec.  VHP ,  UHF ,  and  L-band  back-scattering  experiments 

probe  wave  numbers  given  by  k  =  hirf ±/c ,  where  q  is  the  incident  frequency. 
The  theory,  therefore,  predicts  central  line  returns  with  very  narrow  fre¬ 
quency  widths  for  back-scattering  orthogonal  to  magnetic  field  lines.  The 
width  is  given  by 

_li  p 

toi  s  5  x  10  q  , 

(16) 

q  in  MHz 

Thus,  for  example,  Aw  =  13  Hz  for  a  probe  at  157  MHz  and  Aw  =  1  kHz  at 
L-band.  The  above  results  do  not  take  into  account  the  modifier  or  its 
possible  effects  on  the  intensity  of  the  diffusion  mode.  In  Appendix  E  we 
derive  further  details  of  the  diffusion  mode,  using  more  refined  approxi¬ 
mations.  The  mode  is  found  to  propagate  only  within  10"6  radians  of  strict 
orthogonality  to  the  geomagnetic  field, 
b)  Stimulated  diffusion  scattering 

A  new  electromagnetic  parametric  instability  has  been  found  in  which 

an  O-wave  modifier  acts  as  a  pump  to  drive  unstable  the  diffusion  mode 

described  in  the  previous  note,  together  with  a  scattered  O-wave .  Figure  11 

shows  the  simple  geometry  for  which  the  calculation  has  been  performed. 

The  incident  and  scattered  wave  vectors  are  orthogonal  to  B  and  the  wave 

vector  of  the  diffusion  mode  is  also  orthogonal.  The  lov/est  threshold  is 

found  to  occur  for  back-scattering,  in  which  case,  k  ~  2k 

*  diffusion  pump* 

The  threshold  we  have  obtained  is  less  than  a  factor  of  two  higher  than  the 
threshold  for  the  parametric  decay  instability.  A  pump  with  local  vacuum 
intensity  of  5QiW/m2  incident  upon  the  F-layer  with  a  75  km  scale  length 
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should  be  an  order  of  magnitude  above  the  threshold  for  this  Instability. 

The  parametric  coupling  converts  the  diffusion  mode  Into  a  very  low  fre- 

quency  propagating  mode  of  real  frequency  u>  =  .  The  back- 

scattered  O-wave  differs  from  the  pump  by  this  frequency.  The  pump  wave 

number,  k  ,  near  the  reflection  point  is  much  smaller  than  the  free -space 
*  pump 

wave  number,  and  is  given  by 
“pump  "  -066  “o/c> 

for  a  100  km  scale  length.  For  a  6  MHz  modifier,  the  diffusion  node  fre¬ 
quency  is  on  the  order  of  10”^  Hz  and  the  parametric  growth  rate  is  on  the 
order  of  minutes,  which  appears  long  in  coqparison  with  cannon  spread-f 
onset  times.  V/e  have  also  calculated  the  growth  rate  of  stimulated  Brillouin 
scattering,  and  find  it  to  be  on  the  order  of  200  Hz,  which  gives  a  growth 
tine  of  several  milliseconds,  and  appears  short  compared  to  spread-f  times. 
The  details  of  stimulated  diffusion  scattering  can  be  found  in  Appendix  F. 

A  preliminary  saturation  theory  for  this  instability  is  given  in  Appendix  G, 

-2 

where  density  fluctuations  on  the  order  of  10  times  the  ambient  density 
are  found. 
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iH|cm 


Table  l 


Values  of  the  scaling  parameter  n  appearing  In  the  total  plasma  wave  energy 
E  or,  equivalently,  the  absorptive  nonlinear  conductivity 


n 

Theory 

l 

F 

Resonant  approx,  for  f(k),  cutoff  at  k  .  Analytic  saturated  theory 
of  Ref.  6. 

3 

F 

Resonant  approx,  for  f(k)  no  cutoff  at  k  .  Analytic  saturated  theory 
of  Ref.  9. 

Resonant  approx,  for  f(k),  cutoff  at  k  .  Time-dependent  numerical 
Integration  (B.  Godfrey,  D.F.  DuBois,  M.V.  Goldman16) 


g 

If  Exact  f(k).  No  cutoff  at  k  .  Time-dependent  1-D  numerical  inte- 

c17 

gration.  W.  Kruer,  E.  Valeo  . 


?  Exact  f (k) .  No  cutoff  at  k  .  Time-dependent  2-D  numerical  inte- 

ciq 

gration.  J.  Fejer,  Y.Y.  Kuo  . 
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Value  used  by  F.  Perkins  and  by  G.  Meltz  and  N.  Tbmljanttvlch  in 
recent  geometric  optics  approximations1^  to  the  nonlinear  absorption 
problem. 


Figure  1:  The  angular  spectrum  (equation  (8))  from  orbit-perturbation  theory 
as  a  function  of  the  angle  between  the  punp  field  and  k. 

Figure  2:  Total  wave  energy  E/nO  from  orbit-perturbation  theory  versus  the 

C 

pump  ratio  P  for  k^/kp  -  0.2,  or  r  *  1.  (see  equation  9) 

Figure  3:  Linear  Langmiir  wave  danping  mechanisms  as  a  function  of  altitude 
below  the  reflection  height.  The  electrons  producing  airglow  are 
accelerated  from  the  landau-damping  region  where  orbit-perturbation 
saturation  applies . 

Figure  Geometry  for  normal  incidence  of  electromagnetic  wave  on  a  plane¬ 
layered  inhomogeneous  plasma.  The  spatial  coordinate  z  is  in  the 
dimensionless  units  of  c/wq,  the  free-space  wavelength  over  2ir. 

Figure  5-'  Integrated  absorption  1  -  R,  R  =  |E  ^/E^l  ,  for  incident  intensity 

I.  =  cE?  /8ir  =  50tW/m2  (±25)  and  various  values  of  the  nonlinear 

me  me 

scaling  parameter  n  defined  In  equation  (6)  and  Table  1.  Both  linear 
and  nonlinear  absorption  are  included  and  the  following  parameters 
are  used:  n  *  5  x  10^  cirf-^,  0e  =  0.2eV  (2319°K), 
vgi  «  350  Hz.  The  ions  are  assumed  to  be  oxygen.  Scale  heights  of 

ii 

37.5  km  (zQ  =  5000)  and  75  km  (zQ  =  10  )  are  used.  The  right  side 
scale  on  each  plot  is  R  -  R11npar,  which  is  a  measure  of  the  rela¬ 
tive  percent  of  nonlinear  absorption. 

2 

Figure  6:  Integrated  linear  and  nonlinear  absorption,  1  -  R,  R  *  |Erefil  ' 

| Einc 1 2  as  a  function  of  incident  Intensity,  ^  3  cE^^/Stt .  The 

nonlinear  scaling  parameter  0  is  defined  in  equation  (6)  and  Table  1. 
The  scale  height  is  H  *  75  km  and  the  integrated  linear  absorption  is 


shown  as  well. 


Figure  7 


Figure  8: 


Figure  9: 


Figure  10: 


Figure  11: 


:  Pump  radiation  electric  field  amplitude  squared  (in  units  of  incident 
arplitude  squared)  versus  altitude  (in  units  of  c/u  =  7.5  meters). 

Ix»rer  altitudes  are  towards  the  right  and  the  lowest  altitude  shown 
is  *1  km  below  the  reflection  height.  Linear  and  linear  plus  nonlinear 
cases  are  plotted  for  an  incident  vacuum  intensity  of  50pW/m2,  a  scale 
height  of  H  =  75  km  (zQ  =  10**),  and  a  nonlinear  scaling  parameter 
n  -  2  (see  equation  (6)  and  Table  1).  The  density  is  5  x  I05cm"3 
and  the  temperature  is  0.2eV.  The  peak  of  the  Airy  function  swelling 
factor  is  shown  by  X  for  comparison.  The  integrated  linear  absorption 

is  1  -  R  =  17?,  whereas  the  integrated  linear  and  nonlinear  absorption 
is  an  additional  33?. 

Negative  of  Poynting  flux  in  units  of  incident;  flux  for  parameters  of 
Figure  6. 

The  component  of  the  diffusion  tensor  along  the  geomagnetic  field 
versus  velocity  along  the  geomagnetic  field.  u>e  is  the  plasma 
frequency  at  the  reflection  point  and  vg  is  the  electron  thermal 
velocity.  D(v)  is  a  maximum  at  v  near  5  thermal  velocities  where 
the  value  of  D(v)  is  about  1020  cm2  sec-3. 

Relaxation  of  a  density  perturbation  in  a  homogeneous  plasma  by 
diffusion  across  the  magnetic  field. 

Simple  geometry  used  for  stimulated  diffusion  back-scattering. 

The  diffusion  mode  has  twice  the  wave  number  of  the  incident  pump. 
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AFHSNDIX  A 


NONLINEAR  SATURATION  OF  PARAMETRIC  INSTABILITIES : 
8IECTRUM  OF  TUREULENCE  AND  ENHANCED  COLLISION  FREQUENCY 


Jerome  Wclnstock  and  Ban del  Bezzerldes 


ABSTRACT 


We  consider  the  nonlinear  saturated  state  of  parametric  instabilities 

driven  by  an  externally  applied  electromagnetic  pump  field  oscillating 

near  the  plasma  frequency.  For  »nort  wave  lengths  or  strong  pump  fields > 

modifications  of  electron  and  ion  orbits  are  the  dominant  saturation 

mechanism.  A  calculation  is  made  of  the  saturated  spectrum  of  Langmuir 

waves  excited  by  the  decay  instability  for  T  wl.  The  spectrum  is 

found  to  be  broad  in  angle,  in  the  absence  of  a  magnetic  field,  and  this 

has  been  confirmed  by  computer  simulations.  A  calculation  has  also  been 

* 

made  of  the  enhanced  collision  frequency  and  heating  rate  v  .  It  is 
found  that;  (a)  v  increases  rapidly  with  pump  intensity  when  the  intensity 
is  small,  (b)  v*  increases  very  slowly,  logorithmically,  with  intensity 
when  the  intensity  is  large,  and  (c)  v  significantly  depends  on  whether 
the  ratio  of  matching  wave  number  to  Debye  wave  number  is  greater  or 
less  than  0.2.  These  predictions  have  all  been  confirmed  by  computer 
simulations  and  by  a  laboratory  experiment.  A  previous  calculation  of 
the  authors  for  the  saturation  of  the  decay  and  oscillating  two-stream 
instabilities  with  Tfi  *  T^  was  limited  to  Just  the  total  fluctuation  energy 
at  saturation. 
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I.  IHTRCDUCTION 


The  theory  of  parametric  excitation  of  waves  in  a  plasma  is  of  consider¬ 
able  interest  because  of  its  role  in  laser  heating  and  ionospheric  modifications. 
A  problem  of  particular  current  interest  is  the  nonlinear  stabilization  o 
the  excited  waves  once  the  instability  threshold  is  exceeded.  Data  is  now 
available  from  a  wide  range  of  experimental  situations  which  provide  valuable 
information  about  the  final  saturated  state  of  the  plasma.  Computer  Simula- 

1  2  3 

tion  studies,  ionospheric  incoherent  scatter,  airglow  measurements,  and 

U- 

conventional  laboratory  plasma  experiments  all  suggest  the  existence  of  a 
saturated  state  of  highly  excited  waves.  Clearly,  any  attempt  at  a  theoretical 
understanding  of  this  state  must  provide  the  nonlinear  mechanism  which  finally 
limits  the  level  to  which  the  turbulent  fluctations  can  grow.  In  particular, 
the  saturated  state  of  the  parametric  decay  instability,  induced  by  an 
incident  electromagnetic  wave  EQsin  u)Qt,  bafl  been  the  object  of  intense 
current  research.  Several  nonlinear  mechanisms  have  been  presented  to  limit 
wave  growth.  Pustovalov  and  Silin  have  considered  the  nonlinear  damping 
of  acoustic  waves  by  induced  scattering  frcm  ions  for  the  case  Tg  »  T^,. 

Kruer  and  Dawson1  have  considered  the  effects  of  particle  trapping  for 

6  7 

strong  instabilities.  DuBois  and  Goldman,  and  Valeo,  Oberman,  and  Perkins' 

have  calculated  the  effects  of  scattering  of  the  parametrically  unstable 

g 

Langmuir  waves  by  ion  density  fluctuations.  The  authors  have  calculated 
the  role  of  the  deflection  and  heating  of  electrons  and  ions  by  unstable 
waves,  both  Langmuir  and  acoustic.  Saturation  is  achieved,  with  this  mechanism, 
when  the  rate  at  which  the  parametric  instability  feeds  energy  into  the  unstable 
waves  is  equal  to  the  rate  at  which  scattered  electrons  and  ions  remove  energy 

g 

from  the  unstable  waves.  It  is  interesting  to  note  that  in  the  former 
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mechanisms^'  it  i3  assumed  tfiat  the  particle  orbits  are  unperturbed  by  the 

presence  of  unstable  waves.  In  the  latter  mechanism,  on  the  other  hand,  the 

distortion  of  particle  orbits  by  unstable  waves  is  the  crucial  nonlinear  effect 

leading  to  saturation.  A  comparison  of  the  results  of  different  mechanisms 

will  be  presented  in  Sec.  V.  The  nonlinear  saturation  due  to  the  scattering 

0 

of  particles  by  unstable  waves  was  calculated  by  the  method  of  perturbed 
orbit  theory.  The  results  were  compared  with  recent  computer  simulation 
studies  of  Kruer  and  Dawson.^"  Quantitative  agreement  for  the  total  saturated 
wave  energy  was  obtained  over  a  wide  range  of  pump  powers  which  encompassed 
both  the  parametric  decay  and  the  oscillating  two-stream  instability  for  the 
case  Tg  »  ^ . 

In  this  paper  we  complete  the  task  begun  in  our  earlier  work,  i.e.,  we 
present  the  solution  to  the  more  comprehensive  problem  of  determining  the 
spectral  distribution  of  excited  waves  saturated  by  the  effect  of  perturbed 
orbits.  Here  we  treat  specifically  the  parametric  decay  instability  spectrum 
for  approximately  equal  ion  and  electron  temperatures.  As  will  be  clear 
subsequently,  the  approach  is  not  limited  to  this  case,  and  we  will  report 
the  detailed  spectrum  results  for  other  regimes  in  a  future  publication.  In 
fact,  the  general  procedure  is  applicable  to  a  wide  assortment  of  problems 
whenever  strong  wave-particle  interactions  contribute  to  the  saturation. 

We  find  as  the  principal  result  of  our  work  an  angular  spectral  intensity 
l(x)  of  Langmuir  waves  which  is  given  by 


where  x  is  the  cosine  of  the  angle  between  the  pump  and  the  wave  vector  k, 

P  is  the  ratio  of  the  pump  intensity  to  threshold  intensity  for  instability 
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at  the  matching  wavenumber  k^,  n0fi  Is  the  kinetic  energy  density  of  electrons, 

and  r  =  r(k  )  is  the  ratio  of  the  colllslonal  to  Landau  damping  decrements  for 
m 

the  matched  wave  number  k  .  The  value  of  k  is  obtained  from  the  matching 

m  m 

condition  u)q  ■  (k^)  +  u>a  (km) ,  where  ouq  is  the  pump  frequency,  is  the 

Langmuir  frequency,  and  uj  is  the  acoustic  frequency.  The  main  feature  of  this 

spectrum  is  that  it  is  broad  in  angle  with  a  rather  slow  dependence  on  pump 

2  2 

intensity  and  angle,  through  the  simple  factor  Rc  ,  with  a  cutoff  at  Px  si, 

Thus  the  spectrum  broadens  in  angle,  ultimately  reaching  isotropy,  with  increased 

driver  intensity  P.  The  k-spectrum  is  characterized  by  a  half  width  A k  on  the 

order  of  the  linear  width  Ak  of  excited  waves.  The  half  width  can  be  greater 

o 

or  less  than  AkQ,  depending  on  the  ratio  of  matching  wave  number  to  Debye 

wave  number.  These  results  are  in  sharp  contrast  to  those  obtained  by  the 

scattering  of  Langmuir  waves  by  ion  fluctuations.  There  one  finds  the  spectrum 

narrowly  confined  to  the  direction  along  the  pump  polarization  with  spreading 

6  T 

to  lower  k  values  as  the  pump  intensity  is  increased.  * 

A  basic  objective  in  any  calculation  of  the  saturated  wave  state  is  the 
increased  wave  damping  decrement  due  to  the  nonlinearity.  In  Sec.  II  we  derive 
the  nonlinear  susceptibilities  and  dielectric  function  in  the  presence  of  the 
external  pump  EQsin  u>ot.  Nonlinear  susceptibilities  are  determined  by  account¬ 
ing  for  the  perturbed  orbits  of  +>&  particles.  The  effects  of  perturbed  orbits 
are  characterized  by  a  deflection  frequency,  by  describing  the  rate  at 
which  an  electron  is  deflected  by  a  wavelength  due  to  unstable  waves.  This 
frequency  is  explicitly  calculated  in  terms  of  the  spectra  intensity.  The 
nonlinear  damping  decrement  due  to  waves  is  then  determined,  from  the  nonlinear 
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susceptibilities,  in  terms  of  and,  hence,  the  spectral  intensity.  In  steady 
state  the  parametric  growth  rate  is  balanced  by  this  nonlinear  damping  decrement, 
and  results  in  a  nonlinear  equation  for  the  spectral  intensity.  In  Sec.  Ill 
this  equation  is  analytically  solved  for  the  saturated  spectrum.  In  Sec.  IV 
we  calculate  the  anomalous  collision  frequency  and  heating  rate,  and  in  Sec.  V 
the  results  are  illustrated  for  specific  cases,  including  the  ionosphere.  A 
comparison  is  made  between  the  results  of  nonlinear  perturbed  orbit  theory  and 
nonlinear  Landau  damping.  For  small  values  of  P  we  find  that  the  nonlinear 
saturation  of  the  parametric  decay  instability  for  the  particular  case  in  which 
Tg  ~  is  governed  by  perturbed  orbits  for  k  ^  k^/5  and  is  governed  by 
nonlinear  Landau  damping  for  k  <  k^/5,  where  k^  is  the  Debye  wave  number  and  k 
is  the  wave  number  of  excited  Langmuir  wave .  This  demarcation  in  k  is  influenced 
by  the  pump  ratio  P.  For  large  enough  P  saturation  by  perturbed  orbits  is 
Important  for  all  k.  A  brief  summary  of  the  results  is  given  in  Sec.  VI. 


II .  NONLINEAR  DIELECTRIC  FUNCTIONS  AND  WAVE  DAMPING  RATE 
We  consider  an  isotropic,  homogeneous  plasma  in  the  presence  of  a  homo¬ 
geneous,  external,  sinusoidal  field  E^in  wt.  It  is  veil  known  that  such  a 
field  will  excite  parametric  instabilities  in  the  plasma  when  the  intensity  of 
the  field  exceeds  the  threshold  for  instability.  In  particular,  Langmuir  waves 

are  driven  unstable  when  E^/Unn9e  >  ^th^^e  =  ^  ^or  ^e  ^  ^i'  w^ere 

-the  linear  damping  decrement,  9^  is  the  electron  kinetic  energy  and  n  is 

the  electron  density.  These  waves  continue  to  grow  until  nonlinear 
effects  develop  which  saturate  this  growth.  Here  we  shall  calculate  the 
nonlinear  state  of  this  instability  at  saturation. 

The  ultimate  goal  of  our  nonlinear  calculation  is  the  spectral  density  of 
the  longitudinal  electric  field  fluctuations  l(k,u>)  (TV)  L<  |E(k,u>)  |2>, 

where  E(k.u>)  is  the  space-time  Fourier  component  of  the  saturated  field 
fluctuations,  and  T  and  V  are  representative  of  a  time  interval  of  observation 
and  the  volume  of  the  system,  respectively.  For  the  parametric  decay  instability 
l(k,tu)  is  dominated  by  resonances  in  the  neighborhood  of  and  u>q  - 


high 

frequency 


du> 

2ti 


I(k,u>)  . 


I  00 

a  — 


•i 


low 

frequency 


I 


(2) 


Since  the  low  frequency  wave  density  is  precisely  known  in  terms  of  the  high 
frequency  wave  density,11  it  is  sufficient  to  focus  on  the  high  frequency  density 
I^k).  The  spectral  density  I^k)  satisifes  the  kinetic  wave  equation 

5t  IL(-)  =  "  2Y5T  IL(-)  +  2  Sp  (-}  ’  (3) 
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where  Y^(k)  is  the  total  damping  rate  of  the  Langmuir  waves  (the  high  frequency 

waves  whose  spectral  density  Is  I^(k)  will  be  referred  to  as  Langmuir  waves). 

S  (k)  Is  the  rate  of  generation  of  wave  fluctuations  due  to  particle  dis- 
P  ~ 

creteness  effects  In  the  plasma  (thermal  noise).  The  detailed  form  of  S  Is 

Jtr 

not  important  except  close  to  instability  threshold.  To  obtain  we  derive 
the  nonlinear  dielectric  function  In  the  presence  of  the  pump  and  the  nonlinear 
waves.  In  the  Appendix  we  show  that  for  Tg  >  T^, 


YL  “ 


(M 


Here  the  last  term  on  the  right  hand  side  is  the  parametric  growth  rate  in 
which  is  the  total  linear  damping  decrement  with  +  vc>  where 


Debye  wave  number  k^,  and  is  the  collislonal  damping  decrement.  The  pump 
ratio  P  =  P(k)  is  the  ratio  of  actual  pump  intensity  to  the  threshold  intensity 
at  the  wave  number  k: 


\  VmBJ  ' 


(M 


It  is  important  to  bear  in  mind,  later,  that  P  varies  with  k  through  y^(k)  . 

The  resonance  function  f  measures  the  degree  of  frequency  mismatch  for  excited 

9 

waves  of  wave  number  k: 


f  =  f(k)  =  [1  +  (A<ju/2  y  )2  (bj  +  6)2/uj i  6]"1,  k  <  k  (5) 

B.  C 

=  0  ,  k  >  k,  , 

where  Au>  =  (juq  -  u^(k)  -  ««a(k),  and  5(k)  =  u>q  -  (^(k),  from  which  k£  is  defined 
as  6(kc)  =  0. 
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w 

In  Eq.  (4)  denotes  that  part  of  the  damping  decrement  which  is  due  to 

fill  i  o 

nonlinear  wave  interactions.  It  has  been  established  '  '  that  such  nonlinear 

effects  are  determined  by  the  method  of  perturbed  orbits  as 

\  ♦  Y  (|s  He  c®)'1  Tm  t®  (6) 

4 

*”'  =  1  ♦  J  xf  (k,  0.)  (7) 

st 

xf  (k,«j)  *  -i  -2|  Jdt  Jd^v  exp[i(u)  -k-v)t  ] 

k  o 

d(F  (v)> 

x  (exp  [-ik-Sr^^t )  ]>  k- - 57““  (8) 

where  e^1  s  c^(k,u>)  denotes  the  nonlinear  dielectric  function  at  frequency 
and  vavevector  k  and  x^(fc»®)  denotes  the  nonlinear  susceptibility  of  species 

““  S  “ 

s,  u>  is  the  plasma  frequency,  <F  (v))  is  the  background  velocity  distribution 
ps  s  — 

function,  and  fir  (-  t)  is  the  nonlinear  deflection  or  perturbation  of  the  particle 

"C 

orbit  at  time  t  due  to  the  action  of  the  fluctuating  6E.  fir  (t)  is  given  by 

“  1 0 

t 

6r  (t)  -  £-  Jdt •  t’  6E[r  (t  -t’),  t’J  (9) 

— s  m  °  - 6 

3  O 

where  6  Is  the  ch&rge  and  ni  Is  the  ni&ss*  It  should  be  noted  th&t  In 

s 

Eq.  (9),  is  evaluated  along  the  actual  particle  orbit  r^t)  givon  by 

r  (t)  »  r  -  tv  +  fir  (t) 

where  r  is  the  particle  position  at  t  =  o. 
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Equations  (6)  -  (9)  will  be  used  to  determine  ^  in  terms  of  (6E  6E),  and 

then  will  be  substituted  into  Eqs.  (4)  and  (3)  to  obtain  a  closed  equation 

for  (6E  6E).  It  is  clear  that  a  major  task  in  carrying  out  this  program  is  to 

evaluate  the  right  side  of  Eq.  (8)  for  in  terms  of  (6E  6E).  However,  it 

s  —  — 

has  been  demonstrated  previously  that  this  bosk  is  readily  accomplished  by 
means  of  a  cumulant  expansion  as  follows:^ 

(exp  [-ik-Sr^-T)  ]>  =  exp[-  |  ((k-Sr^C-T))2)  +  C^/n!),  (10) 

where  C  is  the  cumulant  of  the  nth  order  moment  of  ik*6r  and  is  assumed 
n  —  —e 

\ 

negligible  for  n  *  3.  From  Eq.  (9)  the  mean  square  displacement  term  in  the 

12 

exponent  of  Eq.  (10)  is  given  after  transformations  of  the  time  integrals  by: 

2  t  t 

!<(k.6r,(r)2>  =  ^  Jdt  t  Jdt'(t  -t '  )(k‘6E(r,T)  k-6E(r  (T  -t'),  t  ')>  •  (U) 

2  ~s  in  o  o  “s 

g 

Expanding  6E  in  terms  of  its  Fourier  components  E(k,u>)  =  k  E(k,u>)  and  using 

12 

the  spectral  definition  in  Eq.  (2)  we  find  that  Eq.  (ll)  may  be  written: 

dtt  Jdt'  (t-f)  Z  J^4  (k  k')2  Iv(k') 

o  (2tt)° 

x  exp[-i(u>  (k’)-k '  v)t  ’  -  i((k'-6r  (-t'))2>  (12) 

V  “  ^  c  “ 

where  k '  denotes  the  wave  vector  and  v  stands  for  either  the  Langmuir  or 
acoustic  mode.  This  equation  has  previously  been  investigated  in  the  diffusion 
limit  to  obtain  resonant  broadened  wave -particle  contributions.  Here  we  shall 
calculate  an  important  nonresonant  contribution  as  well.  Thus,  in  the  asymptotic 
limit  of  large  t,  compared  to  the  correlation  time  of  the  electric  field  along 


|<(h-6£s(-t))2>  =  4-Ee  / 

m  ~ 
8 


the  particle  orbit,  Eq.  (12)  becomes: 


!<Ql*«£8(-t)12>  -  |  t3  k*2(v)-k  +  |  t2  [^(v)]2  (13) 

where  the  diffusion  term  D(v)  and  the  frequency  (v)  are  defined  by 

<  »  0  A  (31c  *  *  A  ® 

—  —  S~  Re  S  Iv(k’)  Jdf  expf-ifa. Jk’Vk’.vJf.  i<fk‘-6r  (-t)]2» 

m  (2TT)-5  o  ~  -  d  -  ~e 


K  5  ■  §  J-^7C^'k’)2I  (k')  Jdt  t*  exp{-  i fo>  (k’)  -k’-v]  t 

m  (2n)J  o  “  ~ 


|<(k'-6r  (-t')]2» 


(lUb) 


The  resonant  broadened  diffusion  tensor  D(v)  has  been  previously  introduced  by 

Dupree  and  one  of  us.  The  quantity  U^  (v),  on  the  other  hand,  is  introduced 

in  Eq.  (13)  for  the  first  time.  A  short  discussion  of  u£r  and  its  relative 

d 

importance  compared  to  k-D(v)*k  is  thus  in  order  here.  To  understand  the 

nr 

quantities  k-D-k  and  u)d  it  is  useful  to  introduce  the  characteristic  deflection 

time,  r  ,  or  deflection  frequency,  uu.  »  t’^  defined  by 

a  d  d  * 

<^5is(-Ta),2>  s  1  us) 

so  that  Td(v)  is  the  time  it  takes  a  particle  to  be  "deflected”  the  distance 
of  a  wavelength  from  its  unperturbed  orbit,  and  u>d  *  ^(v),  is  the  frequency 

n 

ulth  which  such  deflections  occur.  J  ^(v)  is  a  function  of  v  since  the  orbit 
r(t)  is  a  function  of  v.  Now,  in  terms  of  U>d(v)  it  is  seen  from  the  t’  integral 
in  Eq.  (lUa)  that  D(v)  is  essentially  zero  unless  v  satisfies 


|w(k')  -  k*-v|  <  ^d(v) 


(16a) 


(l^a) 


ltenc»- ,  D(v)  Id  referred  to  as  a  resonant  broadened  diffusion  tensor.  On  the 
other  hand,  it  in  seen  from  the  t'  integral  in  Eq.  (lVb)  that  the  frequency 
a)^r(v)  is  large  and  positive  when  v  is  such  that 

|cu(k')  -  k'^vj  >  cOj(v)  «  w(k')  ,  all  k'  (l6b) 


in  which  case  (l4b)  reduces  to 


where  the  subscript  denotes  that  k1  is  restricted  to  |tu  -k'*v|  >  u)^.  Hence, 

t~  [uw“r(v)  ]“■  is  a  nonre.ionant  contribution  to  (Tk-6  r(-x)J^).  We  shall  thus 
a  -  - 

refer  to  u)”r(v)  as  the  nonresonant  deflection  frequency.  It  describes  diffusion 
due  to  the  action  of  waves  nonresonant  with  the  particle  when  | u>  -k'-vl  »  U)^. 
Similarly  one  could  describe  the  resonant  broadened  frequency  u£(v)  by 

«£(v)  2  [|  k-D(v).k]1//3  (18) 


in  which  case  Eq.  (13)  can  be  written  as 


(tk.fir^f-T)  ]2>  ~t3[u£(v)]3  +  T2(u)"r(v)]2  ,  (19) 

j* 

so  that  for  values  of  v  such  that  U>d(v)  is  small,  i.e.,  for  small  D(v),  we 

nr  -m  nT, 

have  ~  .  In  general,  both  resonant  and  nonresonant  terms,  uj*  and  10^  , 

must  be  included  to  calculate  the  susceptibility  in  Eq.  (8).  However,  it  is 

found  that  one  term  will  usually  dominate  the  other  depending  on  the  spectral 

distribution  Iv(k')  and  the  mode  frequency  u>(k).  A  more  detailed  analysis  of 

-£  and  u£r,  and  their  relative  importance  is  given  in  a  forthcoming  publication. 

(The  relative  value  of  and  W\j|r  is  discussed  in  Appendix  C  for  Langmuir  waves . 

It  is  demonstrated  there  that  >  u£  when  < 1 6E | 2 >  (lm^)"1  >  (k/kp)1*  for  the 

case  of  parametrically  excited  Langmuir  waves . ) 


All 


With  Eq.  (13)  and  (10)  substituted  into  (8)  ve  have  for  the  nonlinear 


electron  susceptibility 


KL 


(|)  w 

i  -Er  J  dT  Jd^v  expfi(<u-k*v)T  -  1/3  k'D-k'n 
k  0 


-  l/2(^r)2  t2}  k. 


a(F(v)> 

dv 


(20) 


Wc  can  now  evaluate  for  both  its  real  and  imaginary  parts.  The  real  part 
of  X,IL  determines  the  nonline^  frequency  shift  of  Langmuir  waves.  The  evalua¬ 
tion  of  Re  J®1  is  given  in  Appendix  B  where  it  is  found  that  for  Langmuir  waves 


(21a) 


where  <oE2/  is  the  mean  square  amplitude  of  the  Langmuir  waves .  The  nonlin^r 
Langmuir  frequency  uu^J  is  given,  in  the  presence  of  turbulence,  by 


(21b) 


Equation  (21b)  explicitly  shows  how  the  Langmuir  frequency  shift,  due  to  non- 

2  -1 

linearity,  scales  with  the  wave  energy  density  (5E  )  (4tt)  .  This  shift  is 

small  when  the  total  wave  energy  density  is  much  less  than  the  thermal  energy 
density.  It  should  be  noted  the  frequency  shift  appears  as  an  effective 
temperature  increase  v2  v2(l  +  (&E2)  (4th0)  )• 

We  now  turn  to  Im  x^ftrm  Eq.  (20),  which  determines  the  nonlinear  damping. 
First  we  must  perform  the  t  integration  in  Eq.  (20).  We  find  that  there  are 
two  cases  for  which  this  integration  is  elementary:  Case  (l)  is  for  (i»d  )  » 
(k.D.k)2^.  In  this  case  one  can  neglect  (Vc *D in  Eq.  (20),  so  that  the  t 
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NL 

Integral  ic  a  simple  Gaussian  and  Im  Xe  is  given  by 


In  X. 


NL 


»2 

JS. 

f  8<Fe> 

/tt/2 

(uj  -k*v)2 

k2  * 

/  -  -  “5v“ 

nr,  \  ^ 

"a 

a(.f)y 

(22a) 


Case  (2)  is  for  (k‘D«k)^3  >  ^ .  In  this  case  we  cannot  ignore  the  term 
(k*D.k)T'3  in  Eq.  (20)  since  this  term  dominates  (tu  )  t  .  However,  if  we  con¬ 
sider  the  r  integral  as  a  function  of  v  then  we  find  that  the  main  contributior 
comes  from  the  range  of  v  satisfying  |<u  -k-v|  ^  /2(l/3  k'D-kj'  .  Outside 
this  range  of  v  the  t  integral,  is  negligibly  small  and  even  slightly  negative . 

l/o 

Hence,  we  need  only  consider  v  in  the  range  |u)  -k.v'|  ^  /2(l/3  k’D-k)  . 

For  v  inside  this  range  the  imaginary  part  of  the  t  integral  is  accurately 
approximated  by 


00 

Im^dT  expfi(u)  -k-v)“  -  i  k -P -k  -  “(av^)2  T 2] 


1/3^(173) _ _ 

[(1/3  k-D-k)2'3  +  l/2  C^)2!172 


exp 


!- 


(u)  -k-v)‘ 


4 [1/3  k-D-k)2/3  +  l/2(^r)2]( 


(2a; 


((1/3  k-D-k)1/3  ^ 


where  f(l/3)  is  the  Gamma  function  evaluated  at  1/3.  The  two  cases  (22a)  and 


NL 


(22b)  can  be  combined  into  one  expression  for  Im  Xg  as  follows: 

(u)  -k  .v)£ 

j 


t  NL 
Im  *e  = 


k2  Tv-  3^7 


exp 


with 


«£(v) 


,3/2 

(f  k'D-k) 


♦  [^r(v)l2  , 


(23a) 


(23b) 
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Equation  (23)  reduces  to  case  (l)  for  <u”r  »  (1/3  k-D-k)1/^  and  reduces  to 
case  (2)  for  (1/3  k-D-k)1^3  >  u£r.  It  is  shown  in  Appendix  C  that  case  (l) 
obtains  when  the  total  wave  energy  <«E2)  roughly  satisfies  <6E2)  (4m  0  )_1  > 
(k/kp)*  and  that  case  (2)  obtains  when,  roughly,  <mj2>  (Unn0e)‘1  < 

Wc  now  consider  the  velocity  integral  in  Eq.  (23).  This  integral 
is  complicated  by  the  dependence  of  u>d(v)  on  v  through  the  Doppler  shifted 
frequency  -  k'-v  indicated  in  Eqs.  (14)  and  (17)*  A  significant  simplifi¬ 
cation  x3  possible,  however,  since  in  Eq.  (23)  u^(v)  can  be  approximated  by 
its  value  at  v^  »  o,  where  v  denotes  the  vector  component  of  v  perpendicular 
to  k,  with  v^  the  component  v  parallel  to  k;  that  is,  u^(v)  **•  ^(v^).  The 
simplification  arises  since  the  main  contribution  to  the  integrand  in  Eq.  (23) 

comes  from  v  <  v  =  (0  /m  )^2,  and  v  **•  tu/k.  5  v  ,  where  v  denotes  ••’he  wave 
X  e  e  e  II  P  p 

phase  velocity.  Consequently,  for  v  »  v  we  have  in  Eq.  (23)  v  »  v  . 

P  e  II  x 

and  therefore 


k'-v  (k'-k)  v 


(2k) 


in  'ju,  so  that  oj.(v) 
a  a  — 

reduces  to 


a)  (v  ).  Integrating  <F  >  over  v  ,  we  find  that  Eq.  (23) 
a  II  s  — x 


where 


\  _  (£  ^ 

00 

r 

1  If  V 

-G  (v 
e  tv  « 

~  2  2 
*  kv 

e 

1  V  V 

J  m  II  II 

»oo 

a)  (v  ) 
d  11 

/▼  \2 

2 

(  - 1  »  +  - 

(U 

(1 .  j  V 

W  ^(v,) 

(  vpj 

(25) 


(26) 
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The  integrand  in  Eq.  (25)  is  highly  peaked  in  the  neighborhood  of  v  =  v 
2  2  II  P 

^or  'JJiy^ud  >:>  ^^eed,  f°r  ^  “*  °  the  integrand  peaks  precisely  at  v  *  v  , 

reproducing  the  well-known  linear  result  of  Landau  dampings  In  the  present 
nonlinear  case  the  finiteness  of  ta^  extends  or  broadens  the  integrand  to  include 
velocities  with  v  <  v^  thereby  resulting  in  increased  damping.  (This  broaden¬ 
ing  of  the  integrand  to  wider  ranges  of  velocity  is  referred  to  as  a  broadening 
of  resonant  wave-particle  interactions). 

To  evaluate  Eq.  (25)  we  apply  Laplace's  method  of  asymptotic  integration 
to  obtain: 


Im  X^L(k,ou) 


v 

c 


V  CD  -  ( V  ) 
p  d  c ' 


exp-Ge(vc) 


(27) 


where  primes  denote  derivatives  with  respect  to  arguments,  and  v  is  determined 

c 

by  (v  )  =  o,  or  with  Eq.  (26), 


2\ 

d'v  =v 


(28) 


The  nonlinear  damping  decrement  can  now  be  obtained  by  substituting  (27)  into  (7)  and 
the  resulting  expression  into  (6).  We  thus  find  with  b  Re  e^/dui  =  2ud2  (ou^)"^, 

3  Sri  -  >< 

P  c 


W  NL 


or,  equivalently,  using  the  value  of  y.W, 

Y?  =  Vy(k)  (exp[A(k,v  )  -  a]  -  1) 


(29) 


where 


2 

ir 

P 


v  v2 

r) 

P/ 
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a  =  -  log 


(30) 


Equa  ion  (29)  gives  the  nonlinear  vave  dancing  yt  in  terms  of  v  and  tu  (v  ). 

c  d  '  c ' 
W 

To  obtain  vL  in  terms  of  the  spectrum  we  use  Eq.  (14)  to  express  U)j  in  terms 

of  the  spectrum  l(k)  and  then  use  Eq.  (28)  to  also  express  in  terms  of  the 

.pectrura.  Let  us,  then,  turn  to  Eqs.  (23b)  and  (l4)  so  as  to  evaluate  uj  in 

d 

terms  of  l(k).  To  obtain  u>d(v)  we  must  substitute  Eq.  (l4)  for  u/jr  and  k*D*k 
into  Eq.  (23b).  This  is  done  in  Appendix  D  where  it  is  demonstrated  that 
u>d(v)  is  approximately  given  by 


2  _  ef  f  <&'  M')2  ^  (k») 

d  “  m2  J  ( 2tt)3  (uj  -k'*v)2  +  ^(vc) 


(31) 


(Equation  (31)  is  a  synthesis  of  the  contributions  of  k-D-k  and  u£r  to  <s  .  From 
Appendix  D  one  sees  that  the  contribution  to  Eq.  (31)  from  the  region  of  k' 
for  which  |o)  -k'*v|  > /2  comes  from  a,^r>  whereas  the  contribution  from  the 
region  k '  for  which  |u>  -k '  -v|  < /2  tu,  comes  mostly  from  k-D*k  .)  To  perform 
the  k'  integration  ,  with  v  =  v^k,  in  Eq.  (31)  we  introduce  polar  coordinates: 


(k-k')g  I^k*) 

-  k  'v  ||(k*k ' )  ]2  +  3«d 


(32) 


where  0'  is  the  azimuthal  angle  that  k'  makes  with  the  pump  direction  e  and 

— o 

A  A 

x  ^e  direction  cosine  of  k  1  with  respect  to  the  pump.  We  next 

change  variables  k'  -  -  k'  in  the  second  term  on  the  right  of  Eq.  (32),  use 
I(k')  =  I(-  k ')>  denote  (k'/k)  ■  y  and  (kv^/tu^)  ■  z  to  obtain 
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J  -* 


(z  -  -z) 


[1  -  yz(xx'  +  J\  -x2  Jl  -  x  ,2  cos  0'))2  +  3^/^ 


(33) 


To  obtain  the  last  equation  in  (33)  we  have  explicitly  introduced  the  azimuthal 
symmetry  of  the  spectrum,  about  the  pump,  by  I(k)  **  l(k',x')j  and  we  have  defined 

tr  -r-  /  P  /  2 

the  direction  cosine  of  k  by  x  s  so  that  k-k '  »  xx'  +  /1-x  /1-x  '  cos  <t>' 

2/2 

The  0'  integration  in  Eq.  (33)  is  tabulated  and  yields,  with  (ju^/u^)  <  <  land  yz  5  1, 


w2(v  ) 

dv  r 
< 


ifc  y) 

4rm8 

© 


(yz)2 


(H 


+  yz[xx'(l-yz  xx')  +  yz(l-x2)(l-: 


B- 


(3*a) 


0), 


B2  =  (yz)2  (x  -x')2  ♦  (1  -yz)fl  ♦  yz(l  -2xx')J  +  3 

“l 


(3^b) 
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j.  -*■ 


and  we  have,  for  convenience,  normalized  u>d  to  <1^  and  l(k',x')  to  (Umej.  We 
have  also  reolaced  uWk')  in  Eq.  (33)  by  (^(k)  which  is  Justified  in  Sec.  III-C 
by  the  narrowness  of  the  spectrum  l(k '  ,x  ' )  in  k ' . 

To  perform  the  x'  integration  in  Eq.  (34a)  we  use  the  fact  that  l/B2  is 
sharply  peaked  about  x  =  x ' ,  f or  yz  ~  1 .  In  fact*.  B2  can  be  asymptotically 
expressed  by  B2  ~  (x  -x')2  +  « 

e  -  (1  -yz)  [1  +  yz(l  -2xx '  ]  +  3uo2/u^ 

c  «  1,  y  ~  1,  z  ~  zfi  ~  1,  (1  -yzc)  >  o 

z  ■  kv  /u L  ■  v  /v  (35) 

c  c'  1j  c'  p 

That  y  ~  1  follows  from  the  fact  that  l(k',x')  is  narrowly  peaked  about  k*  -  k 

m 

(see  Sec  III-C)  so  that  y  —  1  for  k  —  km.  Also,  we  have  (l  -zc)  «  1  for 
l(x)  (4nn6e)  1  «  1,  as  is  verified  by  Eq.  (38)-  Hence  (l. -yz)  «  1  for 
z  ~  zc,  k  ~  km-  Since  «<  l  it  follows  that  e  «  1  f or  z  ~  zc,  k  ~  km» 

With  Eq.  (35)  it  is  clear  that  l/B3  »  1  at  x  -  x '  and  that  the  main  contribu¬ 
tion  to  integration  in  Eq.  (34a)  comes  from  x'  ~x,  unless  l(k',x')  is  very 
narrow  in  x'.  In  fact,  l(k',x')  is  found  in  Sec.  Ill -A  to  be  rather  broad  in 
x'  so  that  Eq.  (34a)  becomes 


dk '  l(k',x)  _1 

(2n)2  l*"n9e  ** 


fax' 


Cxx'(l-yz  xx')  +  yz  (l-x2)(l-x ,2)J 


BJ 


(36) 
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The  /.  integral  is  straightforward  with  the  result  >  for  y  ~  1 , 


"d(V  _/'jef/'V\2  _IM  /  l\/z  (l’x2)  +  *B(1-  >\ 

r  vp/U/  (22)  ^(l-zia  ■+“&  x5j  J 

if5  A2  + 


,  (z  ~  zc) 


(37) 


where  l(x)  denotes  the  angular  spectrum  defined  by 

I(x) 


/,2  dkV  I(k  'x) 

2TT  2TT  ’  and 


^  ■  -§  (1-2C)(1  +  zc'av2)  <<:  1> for  « i 

c  e 

o 

Equation  (37)  expresses  «ij(v  )in  terms  of  both  the  angular  spectrum  l(x)  and 
V  However,  «>a(v  )  and  zq  are  related  to  each  other  by  Eq.  (28),  and  substitut¬ 
ing  Eq.  (37)  and,  its  derivative,  into  Eq.  (28)  we  find  the  following  explicit 

relation  for  l(x)  in  terms  of  z  , 

c 


or 


I(*L  _  1  +  1/2  P  _  fl  _  f 

4nn0e  q  zc'1  V  ’ 

(38a) 

■  =  3/2(l-zc)2  ,  for  zc  ~  1 

(38b) 

where  p  and  q  are  near  unity  and  vary  slowly  with  x  and  z  . 

c 

defined  by 

1  z  (l-x2)  +  x2(l-z  ) 

2  P  2  - - T~ 

1  +  z„  -  2z  x^ 
c  c 


They  are 


2  q  =  1  ~  x —  +  x 


(39) 
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?  2 

We  can  use  (38)  and  (39)  "to  eliminate  l(x)  or  z c  from  (37)  to  obtain,  with  !• 


and 


“d<2c> 


(l-z  ) 

•  ■  (1  +  1/2  p)  , 


ft! 


#  (-^r)  VSf  ' 


(UOa) 


(40b) 


which  explicitly  gives  u)^  in  terms  of  l(x).  Finally,  to  obtain  y^  in  terms  of 

l(x)  we  use  Eqe.(38&)  and  (40a)  to  eliminate  (l-zc)  and  uj^(zc)  from  Eqc.  (29) 

w 

and  (30).  The  result  for  y  ,  neglecting  the  small  quantity  a  compared  to  A,  is 


(Ma) 


(4lb) 


where  F  is  an  algebraic  function  of  l(x)  and  x,  whose  magnitude  is  less  than  1/2, 
implicitly  defined  by  Eqs .  (38a),  (40a)  and  (4la),  and  which  rapidly  approaches 
zero  as  l(x)  (4Tm0e)_1  -  o.  F  is  thus  a  minor  correction  and  can  be  ignored 
for  almost  all  purposes. 

Equation  (4l)>  to  our  knowledge,  nives  for  the  first  time  an  expression  for 
the  angular  dependence  of  the  nonlinear  wave  damping  rate  due  to  strong  turbulence 

w 

(perturbed  orbits).  We  note  that  is  a  highly  nonlinear  function  of  the  angular 

spectrum  l(x).  This  result  is  in  marked  distinction  to  nonlinear  Landau  -umping,  for 

example,  in  which  \  is  linear  in  l(x).  It  is  clear  that  a  result  such  as  Eq.  (4l) 

could  not  be  obtained  within  the  confines  of  weak  turbulence  theory.  Having 

obtained  y!  in  terms  of  l(x)  we  can  proceed  to  calculate  l(x)  from  the  wave  equation. 
L 
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III.  STEADY  -STATE  SIECTRUM  OF  WAVES 
A.  An <r j i a r  Spectrum 

The  object  of  this  section  is  to  calculate  the  saturated  angular  spectrum 
from  the  wave  equation.  To  do  this  we  substitute  given  by  (4l),  into  the 
wave  equation  (3)  at  steady  state. 

\  St  =  ■  (y£  +  VL)  1(h)  +  VL  Pf(k)x2  +  S(k)  -  O  (42) 

It  is  seen  in  Eq.  (42)  that  the  increased  damping  due  to  nonlinear  waves, 

W 

V  ,  plus  the  linear  damping  rate,  yL,  must  be  sufficiently  large  to  balance 
the  growth  rate  y^Pfx  for  each  mode  (k).  More  precisely,  y^  must  be  slightly 
larger  in  order  to  overcome  the  noise  source  S-  S  is  only  important  near 

14 

threshold  where  we  can  take 

S(k)  ■  ‘<"9.  yl  ji  ♦  5-^rrr  |  s  s(k,x)  (43) 

(The  first  term  in  brackets  is  the  usual  rate  of  generation  of  Langmuir  waves 
due  to  Cerenkov  and  Bremstrahlung  emission.  The  second  term  is  the  Cerenkov 
emission  at  the  low  frequency  coupled  through  the  pump  to  act  as  an  additional 
enh^'"ed  source  for  Langmuir  waves  and  significatly  dominates  the  first  terra.  ) 
At  steady-state  Eq.  (42)  can  be  formally  solved  for  the  spectura  as 

Kk.k)  -  -r, - - 5  (44) 

To  solve  Eq.  (44)  approximately  for  the  angular  spectrum  l(x)  we  first  add 
and  subtract  y^,  to  the  denominator  of  the  right  side  of  Eq.  (44)  and  divide 
l  the  numerator  and  denominator  by  y^  4  yL  (pfx2  -  l).  The  result  may 
be  written: 
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(45) 


W  J 

\  +  V 


Y,  +  YL(Pfx  -1) 


-  1 


WFi  ’ 


where  I  is  defined  by 
o 


I 


o 


S _ 

yt  +  yL(Pfx2-l) 


(46) 


and  is  of  the  order  of  magnitude  of  the  thermal  noise  level  of  spectral 

intensity.  For  unstable  modes  we  must  have  l(k,x)  »  Iq,  particularly  near 

the  peak  value  of  I(k,x)  at  k  «k  .  Hence,  for  k  «k  Eq.  (44)  loads  to 

in  m 


y£  (k,x)-YL  (k)  (Pfx2  -  1)  «o,  ^|  «  1  ,  (47) 

This  merely  states  that  at  saturation  the  total  growth  rate  of  the  strongly 
excited  modes  Ls  near  zero.  Let  us  define  r  3  r(k)  3  vc/v^(k)  to  be  the  ratio 
of  the  collision  damping  to  Landau  damping  decrements,  and 


Yl  (k)  =  (1  +  r)  (k)  . 


(48) 


Substituting  Eq.  (48)  in  Eq.  (47)  at  k  -  k  and  using  the  expression  for 


m 


\  (km,X^  given  Eci*  (4l),  wi'th  f(k^)  =  1»  we  obtain  from  Eq.  (47) 


exp  |(‘ 


f)i[f(r)  +  3]  ^  j-  <1+r>  f*2  + 


r  =  o 


(49) 


Solving  Eq.  (49)  for  l(x)  we  finally  have  the  desired  expression  for  the  saturated 
angular  spectrum  l(x): 

[log  f  (l+r)  Px2-r)]2  ,  ? 

t  I(x)  _  ET _  ^  - -  (l+F)-1  ,  Px2  2  1  (50a) 

‘  '  fe)'  ’ 


4nn0 
e 


72  i2 
*  1 
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(50b) 


|I  flog  ((l+r )P  x2  -  r}]2 


since  o  *  F  <  l/2 


p 

providing  P  x  >1.  Equation  (50)  describes  the  saturated  angular  spectrum 
I(x)  of  parametrically  excited  Langmuir  waves  in  terms  of  the  ratios  P, 

an<*  r*  Note  that  the  saturated  spectrum  is  spread  over  all  angles  for 

p 

which  P  x  >1  and  is  thus  quite  broad  in  angle  for  P  *  2.  Furthermore ,  this 
angular  spread  coincides  with  the  angles  that  are  linearly  excited  by  the 
instability  (P  x  >  l).  This  result  Justifies  our  previous  assumption  of  a  broad 
spectrum  which  was  used  in  Sec.  II.  Note,  too,  that  as  P  increases  the  spectrum 
tends  to  become  isotropic . 


2  Total,  Wave  Energy  Density 


The  total  saturated  wave  energy  density 


is  obtained  by  substituting  Eq.  (49)  in  Eq.  (5l)  with  the  result 


6E2) 

urn  6 


e 


which  becomes 


<6E2) 

e 


(2TA) 


[log(l+r)p]2  -  li[log(l+r)P  -  2)  .  1  (52) 

/(l+r)P  ) 


for  (l+r)P  »  r.  The  total  wave  energy  is  thus  seen  to  increase  rather  slowly 
with  the  pump  power  ratio  P  when  P  is  large,  and  to  be  quite  small  in  the  Landau 
damped  regime,  r  <  1,  and  large  in  the  collision  damped  regime,  r  >  1. 
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ctrum 


•  k -Dependence  of  Spe 


Thus  far  ve  have  determined  the  variation  of  the  saturated  spectrum  with 
angle,  given  by  Eq.  (50)  and  the  total  spectral  energy,  given  by  Eq.  (52). 

To  cauplete  the  determination  of  the  spectrum  we  now  consider  tin  dependence 
of  the  excited  wave  intensity  on  wave  number  k.  However,  to  simplify  matters 
we  shall  be  content,  in  this  paper,  to  limit  ourselves  to  the  gross  features 
of  the  k  dependence.  In  particular,  we  shall  only  calculate  the  half -widths 
of  the  spectrum  at  half  maximum.  The  values  of  k  at  the  half -width  shall  be 
denoted  by  k^  which,  in  general  will  vary  with  the  angle  cosine  x.  Hence,  if  we 
denote  the  value  of  k  at  the  peak  of  the  spectrum  by  kp,  then  by  definition 


we  have 


wr^y  =  2  •  1: 

P 

This  equation  car.  be  solved  for  k-  by  making  use  of  Eq.  (44)  for  l(ki  )  and 

*  t 

I(kp).  To  simplify  the  analysis  let  us  add  unity  to  each  side  of.  Eq.  (45) 
and  then  take  the  logarithm  of  both  sides  to  obtain 

(Y?  +  Y.)  r  Yt  P  I  I 

.  Icg[i  ♦  _  (Pfx2-1)]=  iced  ♦  -f)  -  .  (! 


The  first  term  on  the  eft  of  Eq.  (54)  is  determined  by  Eq.  (4l)  as 


~  *  P  (-)’  •  >] 


I(*c) 

mT 

e 


5  A(k) 


Equations  (54)  and  (55)  determine  l(k,x)  in  terms  of  l(x)  or  A(k).  Solving 

(54)  and  (5t>)  for  l(k^,x)  and  l(k^,x)  in  terms  of  A(k)  and  substituting  these 

solutions  into  Eq.  (53)  we  have, neglecting  the  variation  of  I  with  k. 

o  * 
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(56) 


A(k^)  -  log  f(l+r(k^))  P(k^)  f(k^)  x2  -  r(k|)} 

-  2[A(kp)  -  log  f(l+r(kp))  P(kp)  f(kp)  x2  -  r(kp) }  ]  , 
where  we  have  used  Eq.  (48).  Equation  (56)  may  be  rewritten  to  give 

(A(k^)  -  log  f(l+r(k^))  P(k^)  f(k^)  x2  -  r(k^)})  -  (k^  -  kp) 

-  A(kp)  -  log  f(l+r(kp))  P(kp)  f(kp)  x2  -  r(kp)} 

I 

=  fir;x)  (5T) 

P 

For  the  sake  of  simplicity  we  consider  the  solution  of  Eq.  (57)  for  (1+r)  P  »  r. 
Then  Eq.  (57)  becomes 


A(kj)  -  A(kp) 

J(l+r(k  ))  P(k  )  f(k  )  )  I 

4.  in£  /  - E - P.  P  i  0 

|(l+r(k^))  P(k^)  f(k^)  j  l(kp,x ) 

(58) 

with  (A9),  P(k  )/P(k.) 

P  e 

VV  Y^(k^)  (l+r(k^)) 

W  d«-(kp))  ’ 

(59) 

~o  that  Eq.  (58)  becomes 

A(k^)  -  A(kp) 

*  -  ^ 

(60) 

or  with  (55) 

tr  m 

■m 

(61) 

Since,  as  will  be  seen,  the  spectrum  is  relatively  narrow  in  k  we  can  approxi¬ 
mate  k  «  k  in  (6l)  and  expand  the  left-hand-side  about  k  «  k  «  k  .  For 
pm  pm 

l(x)Armee  «  2(3A)3,  we  can  then  write  for  (6l), 
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(62) 


where  we  have  used  from  Eq.  (5)  that 

f(k)  -  -  1  + 


1  .  .  (k  -  km)‘ 


k  -  k 


m 


Ak 


<  1  » 


(63) 


ard  AkQ  is  the  "linear  width"  given  by 

*0-  ■  (*) 

For  Ak/k  «  2/3  and  I  /i  (k  ,x)  «  1,  Eq.  (62)  is  solved  to  yield  two  solutions 
m  o  m 

for  Ak  corresponding  to  the  half -width  of  I  on  either  side  of  the  peak  value  of 
I.  Those  solutions  are. 


Ak(l)  «  0  ,  (65a) 

^(2)“-^o(i^)(x)3  (65b) 

Equations  (65a)  and  (65b)  show  that  the  k-spectrum  l(k,x)  is  skew-symmetric  about 
km  towards  lower  values  of  k.  For  a  hydrogen  plasma  with  k^/kp  ~  0.2,  we  note 
that  Ak,2j  from  Eq.  (65b)  implies  a  k  width  equal  to  the  linear  width  AkQ. 
Equations  (65a)  and  (65b)  apply  to  the  case  where  I  /I  «  1,  i.e.,  above  threshold 
for  instability.  At  threshold  we  expect,  on  physical  grounds,  that  the  width  of 
the  spectrum  Ak  to  be  equal  to  the  linear  width  AkQ. 
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D.  Fine  Structure  of  Angular  Spectrum 
Thus  fer  we  have  ignored  the  explicit  effects  of  F  on  the  spectrum  in 
Fq..  (4l).  We  have  been  able  to  do  so  because  |f|  <  l/2  and,  generally  speaking, 
F  will  not  alter  the  total  spectral  energy  by  more  than  about  30$.  However, 
for  large  enough  P  we  find  that  F  can  have  an  unexpected  and  surprising  effect 
on  the  shape  of  the  angular  spectrum.  Namely,  F  can  cause  the  spectrum  to 
peak  at  a  substantial  angle  to  the  direction  of  the  pump.  That  F  can  influence 
the  shape  of  l(x)  at  large  P  is  easy  to  understand,  since  for  large  enough  P 
it  can  be  seen  that  I(x)  as  given  by  Eq.  (4lb)  varies  very  slowly  with  x,  and 
indeed,  tends  to  become  isotropic  for  P  »  1.  Hence,  for  large  P,  the  x 
variation  of  l(x)  may  voii  be  influenced  by  the  subdominant  quantity  F  in 
Eq.  (4la).  That  F  causes  l(x)  to  maximize  at  x  /  1,  i.e.,  that  l(x)  peaks  off 
the  pump  direction,  is  illustrated  for  the  particular  case  of  P  =  10  and 

~  0.21.  For  this  case  we  evaluate  F  by  numerically  solving  the  set  of 
Eqs.  (29),  (30),  (38a),  (39),  (40a),  and  (4la)  to  obtain  an  approximate  value 
for  F.  This  numerically  evaluated  F  ij  then  used  in  Eq.  (50b)  to  evaluate 
l(x).  These  numerical  results  are  summarized  in  Fig.  4  where  we  plot  F, 
log(Px  )  and  l(x)  ~  (l  +  F)  1  (log  Px2)2  as  functions  of  0  =  cos"^x.  It  is  seen 
in  Fig.  4  that  l(x)  peaks  at  x  «  0.9  which  corresponds  to  an  angle  of  about  26° 
off  the  pump  direction.  For  weaker  pumping,  smaller  P,  the  peak  in  l(x)  tends 
to  lie  along  the  pump  direction.  It  is  interesting  to  note  that  in  a  recent 
computer  simulation,  by  Biskamp  and  Chodura,1^  of  the  current  driven  ion  acoustic 
instability  the  angular  distribution  exhibits  the  same  structure  as  Fig.  4 
with  a  peak  at  about  an  angle  of  40°  with  the  current  direction.  Although  the 
ion  acoustic  instability  is  quite  different  from  the  parametric  instability, 
it  should  be  noted  that  the  saturation  mechanisms  and  the  angular  dependence  of 
the  linear  growth  rates  are  both  similar  for  the  two  instabilities.^ 
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IV.  ANOMALOUS  COLLISION  FREQUENCY  AND  DISSIPATION  RATE 
With  the  saturated  spectrum  given  by  Eq.  (50)  ve  are  in  a  position  to 
calculate  the  anomalous  collision  frequency  and  dissipation  due  to  the  Langmuir 
waves.  The  dissipation  is  governed  by  the  temporal  evolution  of  the 

Ip 

averaged  electron  distribution  function  (F  (v,t) ) 

d<F  (v,t)>  d<F  (v,t)> 

- t; - -  —  E  sin  #t - «r-^ - 


+  5^  *  57  (66) 


where  D(v)  is  the  resonance  broadened  diffusion  tensor  defined  by  Eq.  (lUa). 
Equation  (66)  has  been  rigorously  derived  in  Ref.  12  for  any  external  field. 

The  kinetic  energy  density,  K.E.,  of  electrons  in  the  rest  frame,  oscillating 
with  pump,  is  defined  by 

K.JI.  =  n  Jdv  |  me(v  -  vj2  <Fe(v,t)>  (67) 

where  Iq  3  -  eE0  cos(tuot)/(imuo) .  Substituting  Eq.  (66)  into  the  time  derivative  of 
(67),  we  use  an  integration  by  parts  to  show  that  the  first  term  on  the  right 
side  of  Eq.  (66)  vanishes  and  we  are  then  left  with 


dK.E. 

"ST" 


-  f  “e  fel  '  Io)2  4  (sW  •  -ST") 

d(F  ) 

■  -  nme  -  v,)  •  D(v)  •  -5^ — 


*•  -  nm  Jdv  v  •  D(v) 


a(?e> 

5v~ 


(68) 
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since  the  main  contribution  to  Eq.  (68)  cooes  frcm  v  ~  v  and  v  <  v  «  v 

p  o  e  p 

Substituting  Eq.  (l4a)  for  g(v)  into  Eq.  (68)  ve  have 


x  k  . 


S<Fe> 


exp[-  i(o>  -  k-v)t  -  |  <[k  6r(-  t)l2>]  . 


But  the  main  contribution  to  the  v  integration  comes  frcm  k*v  sa  U4^  so  that, 
ignoring  tne  acoustic  energy  I  (k),  Eq.  (69)  becomes 


aie.E. 

~5t 


e2n  .  [  %  IOO  t  f 

1  me  i(2n)3  k2  I j  -J  - 


S<Fe> 


*  exp[-  l(»  .  k.v)t  -  |  <[k-6r(-  t)]2)) 


If  we  now  substitute  the  imaginary  part  of  Eq.  (8)  into  Eq.  (70)  we  have 


-1^1  =  f  %  ISlI  u u  I  v^fk  ML  ) 


and,  with  Eqs.  (6),  (7),  and  (47) 


JB  C&<\>  ■  *  Vj/o^ 

=  2fyPf(k)  x2  -  v  j/uj^ 


so  that  Eq.  (71)  becomes 


dK.E.  f  dl; 


Kfc)  ,  a 

Itn  tY,Pr(H)  x2 
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(73) 


^r1  aZjh I(x’  [\  i*2  -  w  ,  sx2  *i  , 

o 

for  the  many  cases  In  vhleh  I(k)  is  narrower  m  k  thln  f(k).  It  ls  nov  Eeen 
that  the  heating  rate  is  affect.!  by  the  angular  .tread  of  the  spectrum.  The 
(angu’ar)  Integral  over  x  with  I(x)  given  by  Eq.  (50)  1.  tabulated  art  can  be 
easily  performed.  Hovever,  it  ls  useful  to  express  the  resulting  integral  in 
terms  of  the  total  wave  energy  density  <|lE|2>  given  by  Eq.  (52).  This  can  be 
approximately  done  for  (l  +  r)P»  r  to  yield 


dK.E.  2  <|5eI2> 


(T^) 


where  <|5E|2>  ia  given  by  (52),  and  or  is  given  by 


«  -  1  +  (1  -  p"3/2)  ___ 


(74b) 


lnlll>-r)PJ  • 

The  quantity  a,  which  comes  from  the  angular  dependence  of  I(x),  varies  slowly 

between  1  and  about  2.  If  we  define  the  net  dissipation  rate,  or  enhanced 
collision  frequency,  v*  by 


dK.E. 


v*  E2 


ot  8rT 

then  we  have  with  (74),  (52),  and  (4a) 

“L 


(75) 


* 

v  x  or 


12  4nh 


(76a) 


ss  U  — 
12 


(^ff9/2)a  ((l0g[(1<r)P,)24Ogt(1+r)P)  *  2) 


8 


/(|l+r)p 


} 


(76b) 
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It  is  important  to  note  that  the  dissipation  rate  v  has  precisely  the  same 

form  as  that  proposed  by  Kruer  and  Dawson  to  describe  their  determination 
#  1 

of  v  by  simulation-  In  place  of  y^P  Kruer  and  Dawson  have  the  maximum 
crrovth 

growth  rate  y°  which  here  is  equal  to  y^P.  Two  qualitative  features 

of  this  formula  have  been  observed  by  experiment  or  by  simulation.  In  particu- 
lar,  v  increases  very  slowly  with  P  for  large  P.  This  result  may  be  viewed 
as  a  saturation"  of  the  enhancement  of  the  collision  frequency.  Such  a 

"saturation"  for  v  has  been  experimentally  observed  by  Ingraham;  secondly, 

* 

v  is  proportional  to  the  maximum  growth  rate,  approximately  given  here  by 
Y[P.  This  behavior  has  been  observed  by  computer  simulation.  Quantitatively, 

*  i 

the  value  of  v  observed  in  computer  simulations  is  found  to  be  in  agreement 
wit n  Eq.  (76). 
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V.  TOTAL  WAVE  ENERGY— EXAMPIZS 


A.  Ionosphere 

It  is  useful  to  consider  the  total  energy  of  waves  for  some  specific  cases 

since,  generally  speaking,  transport  phenomena  depend  on  the  total  energy.  In 

Table  1  we  show  numerical  results  for  the  total  energy  density  of  waves  given 

by  Eq.  (52)  as  compared  with  n9  .  Three  matching  wave  numbers  are  considered 

e  * 

corresponding  to  the  Landau  damped  regime,  the  colli  si  anally  damped  regime  and 

the  region  of  transition  between  the  two.  In  this  i. lustration  of  Eq.  (52)  we 

use  conditions  relevant  to  the  F  region  of  the  ionosphf'e  where  parametrically 

excited  waves  have  been  observed,2  viz.,  v  /u>  •*  2  x  lO-'*,  kT.a“  3  cm”1.  We 

c  pc  u 

note  that  the  total  wave  energy  density  is  very  sensitive  to  k^/^  in  the 
neighborhood  of  v  =  1,  corresponding  to  «  5-  Indeed,  the  energy  changes 

by  at  least  an  order  of  magnitude  for  k_A  =  k-5  to  k_A  *  6  5.  The 
saturation  is  quite  strong  for  r  «  1,  and  is  much  weaker  for  r  »  1. 


B.  Computer  Simulations 

g 

It  has  already  been  verified  that  the  total  saturated  wave  energy  predicted 
by  perturbed  orbit  theory  is  in  good  agreement  with  computer  simulations.1  For 


the  sake  of  completeness  it  is  useful  to  make  a  typical  comparison,  here,  between 
Eq.  (52)  and  the  total  wave  energy  found  in  computer  simulation.  [Computer 
simulations  are  often  dominated  by  the  oscillating  two-stream  instability. 

However,  it  can  be  proven  that  (52),  although  aerived  here  for  the  decay  instability 
also  applies  to  the  oscillating  two-stream  instability  for  Tg  »  T^  as  well  as 
at  Tg  rsT^  This  is  indicated  by  the  results  of  Ref.  8.]  Consider,  then  the 
case  Eo//(4nn0e)  =  0.5  given  by  Kruer  and  Dawson.1  In  this  case  (l+r)P  may 
lie  between  100  and  200,  and  kpAm  w  0.25 •  Substituting  these  values  into 
Eq.  (52)  we  easily  find  that 
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0.17  for  (l+r)P=100 
0.23  for  (l+r)P=20C 


«*21 

4  m  0 

C 

This  is  in  clear  agreement  with  the  simulation  value  of  (6E  )(Um0)  =0.2. 

[We  caution  the  reader  to  note  that  acoustic  waves  are  excited  in  the  simulations 

of  Kruer  and  Dawson  for  E  /Jh ml)  <0.^  and  will  contribute  to  saturation.  This 

o  e 

contribution  has  not  been  included  here,  since  Tg  «T^,  but  was  included  in 

Ref.  8.  and  found  to  decrease  <6E2)/4m0  by  about  an  order  of  magnitude 

e 

below  Eq.  (52).  ] 


C.  Comparison  of  Saturation  Mechanisms 
To  obtain  a  clear  idea  of  the  relative  strengths  of  the  saturation  by 
oerturbed  orbits  compared  to  induced  scatter  off  ions  for  the  parametric  insta¬ 
bility  wt  consider  the  total  wave  energy  predicted  by  the  two  mechanisms. 

In  Table  2  we  show  the  total  wave  energy  for  various  values  of 

P  =  E2/E2  .  where  E2  is  the  threshold  field  in  the  collision  damped  region, 

c  o  the’  the 

The  value-:  in  parentheses  are  the  total  energy  as  predicted  by  the  mechanism  of 
induced  scattering  off  ions.17^  It  is  apparent  that  the  perturbed  orbit  ir<*chanism 

is  dominant  for  r  ^  1.  On  the  other  hand,  for  r  »  1  the  perturbed  orbit 
mechanism  predicts  larger  energies  than  the  mechanism  of  induced  scattering 
from  ions.  Hence  the  transition  between  the  two  regions  of  linear  damping 
establishes  a  simple  demarcation  between  the  regions  where  the  saturation 
by  perturbed  orbits  is  dominant  or  not. 


— 

p 

0.22,  r  «  1 

*  0.19,  r  ■  1 

km^D  "  0,15 »  r  >>  1 

5 

0.0062 

0.011 

0.14 

0.017 

0.022 

0.18 

20 

0.035 

0.038 

0.21 

30 

0.049 

0.05 

0.23 

Table  1.  Total  wave  energy  <6E2>(Utt)'1  In  units  of  n6e  for  range  of  pump 
intensity  P  and  matched  wave  number  kffl  as  predicted  by  perturbed  orbit 
saturation  for  vc/uj  *  2  x  10 


pc 

-  «-a.  r  -  1.6 

A,  ■  °'23>  r  -  °'3 

10 

.023  (.024) 

.0028  (.007) 

~  0  (~  0) 

20 

.043  (.110) 

.010  (.076) 

.001  (.004) 

40 

.070  (.290) 

.024  (.279) 

.006  (.160) 

Table  2.  Total  wave  energy  (6E2>(Utt)’1  in  unites  of  n6fi  for  range  of  pump 
intensity  P  and  matched  wave  number  km  as  predicted  by  perturbed  orbit  (induced 
scattering  from  ions)  saturation  for  vc/u)  =  3  x  10  . 
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VI.  SUMMARY  AND  CONCLUSIONS 

A.  The  primary  aim  of  this  work  was  to  determine  the  saturated  spectrum 
of  the  parametric  decay  Instability  and  to  show  that  this  spectrum  is  broad 
in  angle.  The  nonlinear  mechanism  responsible  for  this  spectrum  was  the 
scattering  and  heating  of  electrons  by  Langmuir  waves  (perturbed  electron 
orbits).  An  analytical  solution  is  found  for  the  angular  spectrum  in  Sec.  Ill -A 
and  is  given  by  Eq.  (50). 

B.  The  width  of  the  spectrum  in  k  is  determined  in  Sec.  III-C,  where 

the  spectral  half -width  is  found  to  be  proportional  to  the  half -width  of  the 

linear  growth  rate.  The  total  saturated  wave  energy  density  is  given  by  Eq. 

(52)  and  is  in  agreement  with  computer  simulations.1 

* 

C .  The  anomalous  collision  frequency  v  ,  which  describes  the  dissipation 

rate  of  the  pump,  is  given,  in  Jej.  V,  by  (76).  It  is  pointed  out  th**t  the 

*  * 
enhancement  of  v  appears  to  saturate  for  large  pump  intensities  and  that  v 

is  proportional  to  the  linear  growth  rate  of  the  fastest  growing  mode. 

4 

This  behavior  has  been  confirmed  by  experimental  observations  of  Ingraham 
and  simulations-1  The  latter  observations  are  in  quantitative  agreement  with 
Eq.  (76). 

D.  The  angular  dependence  of  the  nonlinear  damping  rate  of  Langmuir 
waves,  denoted  by  y^,  has  been  analytically  obtained,  and  is  given  by  Eq.  (41). 

w 

This  is  the  first  time  (to  our  knowledge)  that  y_  has  been  calculated  fo  the 

il 

nonlinear  mechanism  in  which  Larmuir  waves  deflect  electrons  (perturbed 

W 

orbits).  To  obtain  the  variation  of  y^  with  angle  it  was  necessary  to  consider 
deflections  of  electrons  by  nonresonant  as  well  as  resonant  Langmuir  waves. 

E.  A  comparison  is  made  between  the  present  results  of  saturation  by 
perturbed  orbits  and  the  results  of  nonlinear  Landau  damping.  Th.'.s  is  '(escribed 
in  Sec.  V-C. 
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AHENDIX  A- 1 


To  obtain  v^1  we  consider  the  nonlinear  dielectric  function  in  the 
P 

presence  of  the  pump  for  the  decay  branch  of  the  parametric  instability.  The 
decay  branch  may  be  treated  separately  since  linearly  the  oscillating  two-stream 
instability  and  the  decay  instability  occur  in  disjoint  regions  of  k  space 
corresponding  to  6  <  o,  where  6  =  <i)  -  o^(k),  and  nonlinearly  the  orbit  per¬ 

turbations  do  not  result  in  any  transfer  between  these  regions.  The  decay 

13 

branch  of  the  parametric  instability  is  described  by  the  dielectric  function 


Here  denotes  the  dielectric  function  of  the  plasma  in  the  presence  of  a 

pump  Eq  sin  'JUQt,  where  €  Is  the  usual  linear  dielectric  function  in  the  absence 

of  the  pump,  e  =  1  +  E  x  >  with  x  as  the  linear  susceptibility  for  species  s. 

s  s  s 

The  zeros  of  determine  the  frequencies  and  growth  rates  of  plasma  modes 

Pump 

made  unstable  by  the  external  pump. 

It  has  been  demonstrated  in  a  previous  publication  that  the  nonlinear  dielectric 

function  with  the  inclusion  of  particle -erbit  perturbations  due  to  finite  amplitude 

unstable  waves  has  precisely  the  seme  form  as  Eq.  (Al)  with  the  linear  suscep- 

NL 

tibilities  x  everywhere  replaced  by  nonlinear  susceptibilities  x  >  that 

8  S 


NL  x  NL 

cp™p(-’  1  = 


(k»  -  J/  (xf'tk,®)  -  xf  (£.,*»-«>  ))  |  -jff-i - -  -  (A2) 

(  I  «  (£,u>-u>  ) 


where 


*  1  +  f  X® . 


(A3) 


and 


X^(k,tu)  * 

o 


(|£  09 


.  i  -°|  J  ax  J  a3v  <e-*i-Sie(-T>)k. 


d<F  _> 


k  o 


with  6r  (t)  the  wave  perturbation  of  the  particle  orbit  as  given  by  Eq.  (9). 
To  calculate  y^1  we  note  from  Appendix  B  that,  for  (u)^(o)/kVe)2  «  1, 

1 

Re  ~  Re  x 


so  that,  in  (A2), 


2  2 

fxf  (k,«)  -  xf  (k,»  -«0)f~  pj 


and,  hence,  (A2)  becomes 


.2 2.2  2 

wr  mt  k_/4k  A  x 

•wmp  -  ‘  (t’”)  -  -KT - r 

e  (k,ui  -u>o) 

where  A2  =  (0  /u>  E2  (4mQ  )-1  and  x  -  (k*e  )•  Equation  (a4)  accounts  for 

pe  o  o  e  —  — o 

nonlinearities  in  unstable  wave  fluctuations  through  the  perturbed  orbit  6r(-T). 
The  nonlinear  mode  frequencies  and  growth  rates,  of  the  parametric  instability, 
in  the  presence  of  finite  amplitude  wave  fluctuations  are  determined  by  setting 


€hjmp  *  0  > 


with  complex  solution 


NL  ,  NL 
<«  =  (U  -  iyTl 
P  R 


where  is  the  nonlinear  frequency  and  y^  is  the  desired  nonlinear  grvwth 
rate.  The  problem,  new,  is  to  carry  out  this  solution  so  as  to  obtain  y^  .  As 
one  might  expect,  the  solution  of  (a4)  for  y^  is  similar  to  the  solution, 
given  in  i\ef.  13  of  (Al)  for  the  linear  parametric  growth  rate  y^.  We  repeat 
the  main  steps  as  they  apply  to  (A4).  The  solution  of  (a4)  for  the  zero  of 


^  is  given  in  terms  of  the  zeros  of  e^k, on)  and  c**k(k,tt>-w  ).  These  zeros 

jump  ~ 

are  the  nonlinear  high  and  low  frequency  modes  that  are  coupled  by  the  pump. 

The  zeros  of  cNL(kf«n)  and  e^fk^  -a»o)  are  the  high  (Langmuir)  and  low  (acoustic) 

frequency  modes  denoted  by 


and 


NL 

“hlGH  =  “l 


“low  ■  •  iya 


NL 


Indeed,  the  calculation  of  and  is  the  subject  of  Sections  II  and  III 

NL 

of  this  paper.  In  the  neighborhoods  of  the  zeros  of  «  we  have 


«“■(*,*(•  -®0»  ~  l(r)  '*o  -<*  (A5) 


Hence  if  u>  ~  uu^  and  if  'J>o  is  chosen  such  that  (V  -u >q  ~  then  (AU)  becomes 


€RJMP  ^ 


1  _ ^ 


*  A2  x2 


o  NL.NL 

^  u>  -U)  +  W  +  iy 


(A6) 


If  we  now  solve  this  for  the  imaginary  part  of  «»  we  obtain 


NL  _  T  -  -  - _  -  _ 


(A7) 


where 


Au>  s 


..NL 


WL~  +  -  W 

Lao 


(A8) 


is  a  nonlinear 


generalization  of  the  "frequency  mismatch".1*4'  Finally,  with 


^  »  Ym'  (at  saturation  ^  ~  o),  Eq.  (AT)  becomes 
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NL  NL 


\  ~  \ 


m 


Y&^a  |  P  f  x2 


t  *  [i  +  (/Vtt/vf)2]-1 


V\ 


so  that  P  will  be  recognized  as  the 


>,  i.e.,  the  ratio  of  the  pump 


intensity  to  the  linear  threshold  value  of  the  pump  intensity.  Equation 
(A9)  expresses  the  nonlinear  growth  rate  of  the  parametric  decay  instability 
in  terms  of  the  nonlinear  damping  rates  of  Langmuir  and  acoustic  waves, 

and  ■  These  damping  rates  are,  of  course,  determined  by  e^k,®)  as, 
(e.g.,  see  (A5)) 


<v]'x 


lm  e^ui^) 


In  deriving  (A9)  we  have  temporarily  ignored  the  effects  of  collisional 
damping.  To  recover  collision  damping  we  note  that  in  the  presence  of  collisions 


v=+v 


M 


Im  e^Cu^) 


(A10) 


where  v is  the  collisional  damping  decrement  ar  '  -  ^(k)  is  the  Landau 

damping  decrement.  Substituting  (A10)  into  (A9)  we  thus  have 


c-\ 


2 

PI  X  , 


(All) 


vhere  k  satisfies  6(k  )  =  o,  and  separates  the  decay  instability  from  the 
c  c 

oscillating  two-stream  instability.  Clearly  for  Aw  «  Eq.  (A13)  reduces 
to  the  previous  form,  while  more'  generally  it  provides  the  necessary  cutoff 
between  the  two  branches  of  the  instability. 
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APffiNDIX 


The  real  part  of  the  susceptibility  is  given  by  Eq.  (2l)  as 

on 

..2 


Pe  X 


NL 


ns  A  ^<F  >  f 

(k,ui)  =  -*vj~  Id  v  k  •  J  dt  Re  exp{i(u)-k  *v)t 


|  k*D(v).k  t3  -  |[«f(v)J2  t2( 


(Cl) 


However,  the  main  contribution  to  the  v  integration  in  (Bl)  comes  from  v  ^  v  . 

This  is  because  the  t  integral  is  a  slowly  varying  function  of  v.  (Note  that 

NL 

this  is  not  the  case  for  1^  xg  since,  then,  the  time  integral  is  sharply 

peaked  about  v  =  v _  »  v  . )  Hence,  from  Eq,.  (lUa),  D(v)  will  be  very  close  to 

p  s  — 

zero  so  that  (Bl)  reduces  to 


Re  Xf,(k,»)  =  -%■ 
5  “  k" 


tt' 


h  -  ¥/ 


dt  Re  exp(i(u)  -k*v)t 


”  k  ll-E(l)-k  t3  "  |(“^r(v)]2  t2) 


w2 

--5 


(z  3<F, 

r  * 


s 1  (tt)  -fc..g) 

^ r 


(tt)  -kv\ 

TT^pT)  ’ 


(B2) 


_  -1  2  “  2 

where  Y(a)  =  or  exp(-  or  )  jdt  exp(t  )  ,  a  function  discussed 

o 

in  Ref.  18.  Since  the  integral  is  dominated  by  v's  which  satisfy  v  ^  v  ,  so 
that  U)  =  ua^  »  k>v,  we  use  the  asymptotic  expansion  for  small  k  -v/u>  in  (B2) 


to  obtain 


Re  X 


NL 


v2  u)2(o) 


tt) 


(03) 


where  tu^(o)  3  u)^(\_)  at  v  =  o.  From  (l4),  we  Bee  that 
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1  ,  2  x  L  2  Ve  <6E2) 

2  (X  3)lUpe  '2  “HIT 


(*) 


Therefore,  (B3)  for  electrons  can  be  written 


2  /  -  2 
NL  V  /,  . 

Be  - *5  I1  4  — 

tu  \  v 

\  P 


1/  2.  lx  <&E2> 


-  ,  X/  C.  ,  X\  ' 

1  4  2(x  4  5) 


) 


(P5) 


Corresr-ndlngly,  the  Langmuir  frequency  Is  shifted  by  the  fluctuating  waves  to 


,  NL^2  2 
(u^  )  =0) 


pe 


H 


.  1/2  lv  <6E  ) 

1  +  2  X  +  3  WS 


(b6) 


Hence,  ve  explicitly  have  a  turbulence  correction  to  the  Langmuir  frequency. 

Thin  turbulence  correction  appears  as  an  effective  temperature  increase  given  by 


2  _  2 
v  ^  v 
e  e 


^  4  If*'4  3> 


Furthermore ,  with  (Lo)  and  the  matching  condition, 

“‘o  “  “t"  (k)  ’ 

it  is  clear  that  increasing  turbulence  will  result  in  the  excited  waves  occuring 
at  longer  wavelengths.  This  prediction  has  been  observed  in  recent  computer 
simulations.1  In  these  simulations,  however,  the  transition  to  larger  wave¬ 
lengths  is  directly  correlated  with  electron  heating. 
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Equation  (Bj)  giver.  Re  x^(k,u)  -U)  )  in  the  neighborhood  of 

3  "*  O 


U)  -ttl  U) 


NL 

Wc  can  evaluate  Xs  (k »UJa)  in  the  same  manner  as  above.  For  present  purposes, 
ve  need  only  use  the  result  that  for  <6E2>/Um0e  «  1  ve  have 


Re  x'>,»a)  *=Re  Xs(K,»a),  ^  «  1 

e 
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APPENDIX  A- 3 


We  wish  to  show  that 


“dW  <  “fW* when  (e)U 


(Cl) 


We  are  only  interested  in  v  s»  kv  since,  as  is  evident  from  (27), 

"  >itf 

only  this  neighborhood  of  v  contributes  to  Im  and  damping.  The 
quantity  is  defined  by  (u£)3  =  (2/3)  k-D-k  with  D  given  by  (l4a) 
t '  integral  in  (l4a)  is  approximately  given  by 

Re  J**  exp^"  1^u1l  "  "  l^-M-t ')i2» 


The 


Vti/2  ou'1  ,  la^  -  k'-v} 


-  k’-v)  >/2 


(U 


d 


UL 


U). 


(C2) 


and  combining  (C2),  (l4a)  and  (l8)  we  have 


f..?\3  /  ti  geV  ft*  ^  C^')  rr_  ,  ,, 

4  ~  2  3»2  JJ^ji  - 5^ - Ht/5u’a  -  l\  -  H'-xH.  (C3) 


where  HfM]  is  the  Heavyside  step  function  defined  by 

!1,  0 

0,  M  <  (C4) 

Since  we  are  only  interested  in  u£  =  u£(v)  for  v  in  the  neighborhood  of  kv,  we 
have,  in  (C3), 

-  k'-v  =  <«L [1  -  k’-k  v/o^] 

*  "  I”  zc(xx'  +  A  -  x2  yT.  x  -2  cos  0  • )  ]  (C5) 
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and 


(k-k')‘ 


k2(xx' 


+  A  -  x2  A~^ 


cos 


rY 


(c  6) 


where  we  have  used  the  notation  of  (33)  in  which  x  and  x'  denote  the  cosine  of 
the  angle  that  k  and  k(,  respectively,  make  with  the  pump  direction,  and  0’ 
denotes  the  azimuthal  angle  between  k’  anu  the  pump  direction.  Substituting  (C5) 
and  (c6)  in  (C 3)  and  dividing  by  uu^,  we  have 


r\3 


K> 


zfi  k*  f  k’2ak'  fix '  f a»'  i(k *,x ' 

3'2  k2  J  2TT  /  2n  J  2n 


2tt 


’Ul 


—  H  |  72  —  -  |l  -  ~  z  (xx '  +  ^£x2  JT-x  *2 
®d  I  ®L  k  cV 


7  [xx  '  +  ;4-x2  Vl-X  '2  COS0’]2 


COS 


0)1 


(C7) 


In  arriving  at  (C7)  we  have  used  the  polar  coordinate  description  as  in  (33) 

Jdk '  =  Jk dk '  dx '  d0'.  Consider  first  the  k'  Integral  in  (C7)-  This  is 
governed  by  the  k'  dependence  of  l(k',x')-  It  is  seen  in  Sec.  IIIC  that  l(k',x’) 
is  zero  exce^  t  for  a  very  narrow  region  of  k ' .  That  is 


l(k  '  ,x  ' )  =0,  k'  s  k  ,  k '  ^  k  -  Ak 

P  P 


(C8) 


uh^re  k  is  the  value  of  k '  at  which  l(k',x')  is  a  maximum,  and  Ak  «  k  so  that 
P  P 


1  ■  ir  * ir  * 1  ’  for  I(k'»x’)  ^ 0 

p  p 

Hence,  if  k  =  k^  then  k'/k  in  (C7 )  is  unity  or  slightly  less  than  unity  and 
(C 7)  becomes  at  k  =  k 

P 

r  1  2tt 

(sj)  4Iv/T^2  fto'fjzr  fxx’  +'^  TlTP  coc0’f 


(C9) 


D  b  o 


~  H  j  £  ^  -  |l  -  zc(xx'  cos0 1 1  J  ,  k  =  kj 


(CIO) 
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To  obtain  (CIO)  we  have  integrated  over  k  '  talcing  advantage  of  the  narrowness  of 

l(k',x')  in  k’  to  set  k'  =  k  in  the  integrand,  and  we  have  used  the  definition 

P 

(2tt)  2Jd'2dk'  l(k',x')  3  I(x').  We  have  evaluated  (uj^/oj  )^,  in  (CIO),  at 

d  h 

k  =  kp  since  it  is  near  k  =  k^  that  the  spectrum  4-s  large .  To  complete  the 
evaluation  of  (CIO)  we  first  consider  the  case  of  x  =  1  to  obtain 


o 


since  the  integrand  of  (CIO)  is  independent  of  0'  when  x  =  1.  We  normalize  u)^  to 
by  substituting  (40b)  into  (Cll),  with  vg/v  =  k/k^,  to  obtain 


I(x-) 


%(1)  ,2 

x 


U) 


h{/2  - 

*  “L 


1+z  x 
c 


■} 


(C12) 


r  / 

Our  object,  now,  will  be  to  determine  under  what  conditions  «<  1.  From 

(C12),  and  the  definition  of  H,  it  is  seen  that  «  0  when 

/2  ~  -  1  +  zc  x’  <  0,  all  x  ’  (C13) 

Ij 

Substituting  (38a)  for  (l-z  )  and  (40b)  for  id  /u)T  into  (C13)  we  have 

c  a  Jj 


U). 


/2  — 


d.zc)  x 


-  (1-x  ’  ) 


(1-x') 


<  0  ,  all  x 1  . 


(Cl4) 


This  inequality  is  satisfied  for  all  x  1  when 


(C15) 
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lienee ,  (C15 )  in  satisfied  and 


Vi, 


H  f/2  ~  -  1  +  zcx'  |  *  0,  all  x'  (Clb) 


when  1(1)  (Unne)'1  >  (2'f/2)  (k/kjj)4.  Substituting  (Cl6)  into  (C12)  we  thus  have 


(4)  “0i 

'  X=1 


(C17) 


Furthermore,  when  cuj  «  0,  then  u>.  u)*?1*,  so  that  (Cl?)  becomes 


(C18) 


Although  (Cl8)  is  limited  to  x  =  1,  it  can  also  be  proven  to  hold  for  all  x. 
To  see  this  we  note  that  if  the  argument  of  H  in  (CIO)  is  negative  for  x  =  1> 
then  the  argument  of  H  is  negative  for  all  x.  Consequently,  if  (Cl8)  is  true 
for'x  =  1  then  it  is  true  for  all  x  that 


(C19) 


Roughly  speaking,  it  follows  from  (C19),  (50b),  and  (52)  that 

p  .  Ij. 

r  ^  nr  \6E  )  ^  /k  \ 

®d  **  ®d  ’  VhCn  >  )  ’  [  } 

since  when  l(l)  (W)"1  =  (27/2)  (k/kp)4  then  it  is  found  from  (50b)  and  (52) 
that  <SE2>  (Unne)"1  «  2(k/kI))4.  The  precise  value  of  <6E2)  at  which 
a)1  -  u>nr  ir.  rather  complex  to  determine. 
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To  derive  (31)  we  first  substitute  (lUa)  and  (lUb)  into  (23b),  and  drop 
the  stun  on  v  since  I^(k')  »  Ia(k')»  to  obtain 


e2  Rg 
n2  (2n)3 


Jd k'  (k-k ' 


>  *,<*•> 


x  expf-  k  ' *v)t '  -  |<[k  '  «5r(-t)  ]2>] 


(Dl) 


In  arriving  at  (Dl)  we  have  assumed  that  varies  slowly  with  k,  compared  to 
the  variation  of  I_(k)  with  k.  The  quantity  ( [k  '  -6r(-t)  ]2)  is  given  by  (19) 
and  to  perform  the  t '  integrations  in  (Dl)  we  shall  use  the  fact  that  the 
orbit  function  expf-  -( [k ' -6r(-t) ]  >}  approaches  zero  very  rapidly  with  t'  for 
t'  >  u>d".  The  t'  integration  can  be  divided  into  two  parts  according  to 
whether  k  '  satisfies 


u)  -  k  '  *v  |  ^  /2  , 

k'eyi 

(D2) 

u)  -  k  '  •  v|  </2  10^  , 

k 'eV2 

(D3) 

For  convenience  of  notation  we  denote  the  regione.  of  k '  space  for  which  either 
(D2)  or  (D3)  obtain  by  V1  and  --espectively.  In  the  region  k'eV.^,  we  have 


j 
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+ 


(U)  -  k'-v)2  -  (9/5 )a^ 

f(ou  -  k’-v)2  -  u >2]2  '  "  *  Vl 

—  —  a 


(u)  -  k'-v)2  -  (9/5  )uj2 
V(.  -  k '  *v)2  -  uj2]2 


k'  e  Vx 


0*) 


where  0(p)  mean  le3S  than  F.  This  result  is  exact  in  the  asymptotic  limit  of 
2 

large  (uj  -  k'*v)  ,  and  is  a  good  approximation  for  any  conceivable  orbit  func¬ 
tion  for  all  k'e  that  is,  for  |uj  -  k'*v|  i/2^.  Note  that  the  number 
(9/5)  in  (D4)  reflects  the  fact  that  when  exp(-  |<[k ' *6r(-t) ]2>)  is  Gaussian 
in  (uj.t')  then  the  second  integral  in  (D1*)  vanishes  for  (uu  -  k'*v)2  (9/5  )uj2. 

To  perform  the  integration  in  (D4)  we  have  used  the  fact  that  the  orbit  function 
approaches  zero  very  rapidly  with  t'  f or  t  ’  >  m"1.  This  integration  does  not 
vary  much  with  the  precise  form  of  the  orbit  function  providing  k'e  and 

the  orbit  function  approaches  zero  rapidly  with  t'  on  the  time  scale  m!1. 

d 

Note  that  the  surviving  contribution  to  the  t  ’  integration  in  (D^)  comes  from 
nr 

the  uj  term.  The  contribution  from  k-D-k  for  k’e  is  negligible. 

For  the  k'  region  |u>  -  k'-v|  < /2  u>d,  the  t’  integration  in  (Dl)  is, 
approximately, 


CD 


(^d  "  t  )  eXP^"  i^<U  "  t')2>) 


g~l/6  r(T/3) 

a,? 


exp 


[- 


(uj  -  k'  'vY 


] 


-1/3 


UJ 


r  (5/3)  r 

I  x  “ 


(uj  -  k '  *v)‘ 


au j! 


] 


,  k'e  V„ 


(D5) 
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To  arrive  at  (D5)  we  have  used  the  fact  that  the  orbit  function  approaches 

zero  very  rapidly  for  t '  >  id"1.  We  have,  in  fact,  used  expf-  |<0c'-6r(-  t')]2>  « 

exp(-  -  ou^  tD).  This  form  for  the  orbit  function  is  adequate  in  (D5)  for  the 

k'  region  |cu  -  k'-v|  ^  /2  u>d.  For  this  region  (D5)  does  not  vary  much  with 

the  precise  form  of  the  orbit  function  providing  the  orbit  function  rapidly 
,  -1  * 

approaches  zero  on  the  time  scale  uj.  .  It  is  of  some  interest  that  both 

a. 

k*D.k  and  furj1")2  contribute  to  (D5) >  and  that  the  contribution  from  (tu^r)2  can 

be  seen  to  be  negative  and  thus  tends  to  oppose  and  weaken  k*D*k  when 

|tu  -  k  *  - v |  <  /2  (u  .  Note,  too,  that  the  contribution  from  (u)^)2  vanishes  for 

(u)  -  k  '  *v)  as  y 2ll).. 

—  —  a 

Substituting  (d4)  and  (D5)  into  (Dl)  we  obtain  u>,  as 

d 


But,  for  k'e  V?  the  following  is  a  good  approximation 
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This  approximation  is  correct  to  within  5#  for  (ou  -k'-v)  =0  and,  at  worst, 

is  correct  tr.  within  a  factor  of  2  for  (ou  -  k '  -v)  =  Sz  ou„ .  In  addition,  for 

“  “  a 

|'i)  -  k'-vj  >  j2  ou.,  the  following  is  a  good  approximation. 


(ou  -  k'-v)2-  (9/5 )u>2 

[(u)  -  k'-v)2  _  w2]2  7 k '  -  v)2+  3u/j 


k'e  v 


(D8) 


Thi3  is  asymptotically  exact  for  IL&rge  (ou  -  k'-v)  and  is  correct  to  within  a 
factor  of  2  for  the  worst  case  of  (ou  -  k'-v)  =  /2  oufl.  Finally,  we  substitute 
the  approximations  (A7)  and  (A8)  and  use 


to  obtain  the  desired  result  for 


uu. 


2,  v  e 

uJd(v)  =  -s - - 

d  m2  (2n)3 


(k-k’)2  I^k') 
(ou  -  k ' -v)2+ 

“  u 


(D9) 


where,  as  has  been  shown,  the  integrand  of  (D9)  is  very  close  to  the  integrand  of 
(D6)  except  for  k'  in  the  vicinity  (ou  -  k'-v)2  m  2uu2  where  they  agree  to  within 
a  factor  of  2. 
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Fig.  1  Plot  of  (log  ePx2)2  (loge  Px2)2  (l  +  F) 
ond  F  for  P=IO  and  kmXD-0.2l 
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CRITICAL  FLUCTUATION  LEVEL  FOR  THE 
PURELY  GROWING  PARAMETRIC  INSTABILITY 


D.F.  DuBois 
B.  Bezzerides 


CRITICAL  FLUCTUATION  LEVEL  FOR  THE 
PURELY  GROWING  PARAMETRIC  INSTABILITY 

I.  Introduction 

Recent  detailed  calculations  of  Harker1,  Meltz2  and  Perkins3  have  led 

to  the  conclusion  that  the  observed  incoherent  backscatter  spectrum  at 

Arecibo  from  the  modified  ionosphere  can  be  explained  in  terms  of  pump- 

enhanced  critical  fluctuation  theory.  This  theory  was  first  developed  by 

4 

DuBols  and  Goldman  who  calculated  the  steady-state  wave-fluctuation  spectrum 

for  a  pump,  near  but  below  the  threshold  for  the  parametric  decay  instability. 

It  appears  that  the  Arecibo  experiments  have  been  carried  out  with  local 

pump  amplitudes  just  at  or  slightly  below  threshold  for  plasma  waves 

propagating  in  the  pump  direction  with  the  wave  number  observed  in  the 

backscatter.  Therefore  in  the  geometry  at  Arecibo  the  waves  observed  in 

the  backscatter,  which  propagate  at  an  angle  of  about  40°  from  the  pump 
2 

polarization  (cos  6  K  ,  have  a  threshold  about  twice  that  of  the  pump. 

The  considerable  detail  of  the  Arecibo  backscatter  spectrum  nresents 
many  spectral  features  to  be  explained  by  theory.  The  critical  fluctuation 
theory  predicts  a  ratio  of  peak  intensity  of  the  enhanced  decay  line  (i.e. 
the  Stokes  or  downshifted  component  of  the  plasma  line)  to  the  ion  line 
which  is  in  agreement  with  experiment5.  Likewise,  the  ratio  of  the  anti- 
Stokes  component  of  the  plasma  line  is  consistent  with  experiment6.  A 
spectral  peak  shifted  by  exactly  the  heater  frequency  has  been  identified6 
as  the  high-frequency  component  of  the  purely  growing  parametric  instability 
which  is  a  different  branch  of  the  dispersion  relation. 

In  this  note  we  compute  the  critical  fluctuation  level  for  the  purely 
growing  instability.  Then  assuming  values  for  the  pump  which  appear  to 


12  3 

explain  the  features  of  the  decay  instability  '  *  we  compute  (.he  ratio  of 
the  purely  growing  line  to  the  decay  line  (Stokes  line).  While  this  report 
was  in  preparation  we  became  aware  of  the  results  for  a  homogeneous  plasma 
which  have  recently  been  computed  numerically  in  detail  for  equal  electron 
and  ion  temperators  by  Hagsfors  and  Gieraltowski  who  used  the  complete 
plasma  dispersion  function.  The  approximate  analytic  result  given  here  seems 
to  agree  with  the  more  complete  numerical  results  and  allow  us  to  take  into 
account  the  altitude  integration  of  the  backscattered  signal. 


II.  Critical  Fluctuation  Theory 

A  general  discussion  of  the  parametrically  excited  wave-fluctuation 

7  4 

spectrum  has  been  given  elsewhere  ’  . 


I(k,u)  -  nu>o) 


lim 


nm' 


m 


|D(k,ai) 


N (k,0)  - 


mo)  ) 
o 


(2.1) 


where  C  is  the  cofactor  of  the  parametric  coupling  matrix  M  and 

nm 

D(k,d))  -  det  (M  ).  The  source  term  K(k,a))  is  the  spectral  density  of  the 
fluctuating  densities 


N(k,U))  -  lim  _  -1  ^  fs  Im  XS(k,u)  (2.2) 

T,V-h»  TV  2ir  s  a) 

where  we  have  used  the  fluctuation-dissipation  theorem  to  relate  this  to  the 
imaginary  part  of  the  electrostatic  susceptibility  arising  from  particle 
species  s  with  temperature  dj .  In  deriving  this  expression  it  is  assumed 
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th«  the  off-diagonal  components  of  ere  reel  so  that  no  mixed  correlation 
terms  of  the  form  <$p<k,(a  -  na>o>«p*<k.«  -  n'^)) occur  for  n*.\ 

lo  treat  the  purely  growing  mode  It  Is  well-known8  that  the  three 
Fourier  components  at  a,  a,  -  «  and  u  -  2U  are  coupled.  The  matrix  H 

U  WQ 

in  this  case  has  been  shown  to  be*’7 


e(k,w)  +  I2 


M  - 


-Kkp/k) 

-I2 


-Kkjj/k) 

e(k,uj  -  a)  )-2 12 

O 

ICkjj/k) 


where  ^  aV;  A2 


47m©  ;  y 
e 


COS0 


(2.3) 


KOCjj/k) 

e(k,w  -  2w  )  +  I2 

O  ml 


(2.4) 


HereJSo  is  the  pump  amplitude  and  y  is  the  cosine  of  the  angle  between  E 
and  k.  We  assume  that  Afr/tyn  «  1  so  that  the  Bessel  functions  which  ° 
appear  in  the  general  expression  for  can  be  Taylor  expanded7.  If  we  use 
this  explicit  form  for  we  can  evaluate  (2.1)  to  obtain 


I(k,w) 


- 


k  |D(k,u)  |  ‘ 


.  4 _  2 

(— )  lJiJA~^2!  N(k»ai> 

\k  / 


(2.5) 


$ 


<7  |JAI2  A2  Nft.c  -  <so)  +  (£>)  |jt  .  i|2  k.u  -  2c  ) 


where 


i  Js  -  Js  +  7^ 


Js  -  e(k,aj) 

JA  -  e(k,<u-2a)o)  ;  JA  -  J  +  A2 


(2.6) 


*  y~~  e(k,a)-a>  ) 

V 

e 
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In  this  notation  the  determinant  is^ 


D(k’W)  *  [JSJAJi  +  +  0(  ^  A  )] 

*1) 


(2.7) 


For  both  instabilities  the  high  frequency  waves  are  near  enough  to  the 
Langmuir  frequencies  i^OO  to  use  the  expansions 


JS  “  ^  (wA+1yl) 

(2.8) 

JA  “  "  (0>*2a)o+0JL+iYL) 


we  will  show.  We  will  chose 

this  case  is  N(k,u)--<i)  )  which 

o 

use  the  relation 

Yv  -  Oej/3wv)  1  2  lm  xS(k,w  ) 

3 


We  must  be  more  careful  in  treating  as 
<i)  to  be  near  The  dominant  noise  source  in 
is  predominantly  ion  Cerenkov  emission.  If  we 


for  the  damping  of  a  longitudinal  wave  of  frequency  we  can  compute  the  ratios 


N(k,w-u>o) 

N(k,oj) 


N(k,a)-OJ  ) 
_ o 

N(k, w-2a)  ) 
o 


Yr  2 
L  (0 


(2.9) 


For  0e  -  0i,  we  have  Ya  *  and  this  ratio  is  (k^/k2)  (w^)  (w^) .  With 


.  1.  2 

A  <^Ac 


•nd  Jj  s  1,  we  see  that  the  second  term  exceeds  the  first  term  by  the  large 
factor  The  last  term  is  similarly  small. 
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III.  Decay  Instability 

Under  the  conditions  u)^  »  y  it  is  readily  seen  that  J.  »  J  and  the 
dispersion  D(k,w)  -  0  reduces  to 


D(k,o)) 


-  V  [Vi  +  *  <V»]  -  0 


(3.1) 


In  this  case  it  is  valid  to  assume  J±  is  near  one  of  its  zeros  and  approximate 
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-2  (w-w  +o)  +iv  ) 
j  a  0  a  a 

1  (0 


(3.2) 


where  (0^  and  y *  are  the  frequency  and  damping  of  an  ion  acoustic  wave. 
The  dispersion  relation  is  then  a  quadratic  in  (0  with  the  solutions 
(provided  A2  «  l)*’7 


“+  •  “o  -  -n+  '  “l  -  "o  - 


(-$) 


(3.3) 


CO  ““CO  =  — n  ■  CJ  —  iy 

o  -  a  s 


(3.4) 


where  it  has  been  assumed  that  A  «  1,  yL  «  y *  and  (A2/A  2)(Yy/Y  )  «  1. 
The  threshold  value  of  A  2  for  the  decay  instability  in  this  cate  is 


V - ^~2  ;  A  2  ■  16  —  — 

T  f(Aw)y2  c  “l  w« 


(3.5) 


f(Aw)  -  - r - =-  ;  Aw  -  a)  -  ul  -  u 

(Aw)z+y  o  ~L  a 


(3.6) 
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Then  if  we  define  ft  ■  co  -  U)  we  can  combine  these  results  to  write 


D(k,b))  -  D(k,w  -ft) 
o 


V  A 

JA  ^  «+>»«-) 


(3.7) 


Substitution  of  this  result  into  (2.5)  yields 


I(wo-ft) 


(4tt)2  £ 
,  2  16 


(ft-^o-HUi)2  +  YL2(l  - 


A2\2  (ft-w  )2  +  y  2 

rr  a  Ta 


-T  (3.8) 


Since  «  Y a  the  first  ft  dependent  denominator  is  sharply  peaked  at 
^  m  ^  (especially  near  threshold)  and  we  can  write  approximately 


T/..  N  _  (Att)2  A2  “l  “a . 

I(w^)  1  16  2  2  7  f  (Aw)N(k,w  -OL}TT6(ft-a)  4ul  ) 

V  ^a-AVV^2 

These  are  the  familiar  results  for  the  decay  instability.* 

IV.  Purely  Growing  Instability 

Using  the  definitions  (2.6)  and  the  approximations  (2.8)  we  can  write 
the  dispersion  relation  as 


k2  «2-<S!-1y  )2 

— 7  B(k.u-n)  "  - 5 - —  J.  +  A2  (J  -1) 

V  °  fo^M)  1  1  “l 


(4.1) 


where 


6  =  wo  - 


(A. 2) 
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N«xt  we  assume  that  |fl|  «  kv^  when  is  the  ion  thermal  velocity.  Then  it 
can  be  shown  that  we  can  write  approximately 


J 


1 


e(k,0) 


(4.3) 


Then  we  can  write 


where  0  is  the  ion  to  electron  temperature  ratio 


(4.4) 


(4.5) 


With  this  approximation,  the  determinant  is  again  quadratic  in  ft  and  can 
be  factored  into  its  roots 


— V  D(k,iu0-n)  .  CJ.n+)  (SM!_) 

\  “l 


where 


fl  -  iy  +  16  J  ^ - 

1  L  f  (-6)  A  2 


and  where 


A  2  -  -8-<1+e)  \  _  (1+B)  ^ 

8  y2  “l  2  Ya  y2 


(4.6) 


(4.7) 


(4.8) 


The  root  fi+  attains  a  zero  imaginary  part  for  the  minimum  value  of  A2y2  when 
6  -  -yL  and  this  minimum  threshold  is  just  A  2  No  instability  is  possible 

o 

unless  6  <  0  i.e.  (0q  <  (^(k).  These  results  agree  with  those  of  Nishikawa8 


and  Sanmartin^. 


Inserting  eqn.  ^4.7)  into  (4.6)  we  find  for  the  spectral  contribution 


of  the  growing  mode  (for  6  <  0) 


_ _ , 

°  i^2  16  (i+8) 2  [n2+(YL-6s)2][n2+(YL+<5s)< 


N(k,n) 


(4.9) 


where 


2\__l 

<-«>  A  2 


(4.10) 


*9  *)  tx  0  0 

where  we  have  replaced  |  J .  (ft)  |  by  |j.  (0)|  -  — r  (6  +YT  )  since  the  other 

V 

factors  peak  at  ft  ■  0. 


This  function  peaks  at  a  0  with  the  peak  value 


K2«V> 

0  l^2  16  (1+8)2  (yl2-62S2)2 


N(k,0) 


(4.11) 


Defining  6  -  6 /Yy  we  can  rewrite  this  using  (4.10)  as 


,/  x  (4tt)2  A2  82  ^  T  1+6 2 

I(w  )  -  J — i— tt - x - x  - 3 — 3 - x 

10  (1+8)  Yl  (146^+26P  V 


N(k,0) 


(4.12) 


where  P  =  A2/A  2. 

g  g 

2  2 

The  factor  in  square  brackets  is  a  maxinun  at  6  “  YT '  lor  P  “  1.  The  width  of  this 

L  g 

—  2  2 

function  is  A6  -  1  (or  A6  ■  YT ) •  Assuming  8  ■  1  we  have  A  »  A  .  If  the 

g  c 

nodes  propagati..*,  in  the  y  -  1  direction  are  taken  to  be  just  at  threshold 
2  2  2  o 

then  A  ■  A^  .  For  y  ■  *5  (0  =  40  )  the  factor  in  square  brackets  reaches  a 
maximum  at 


6  a  -.68  (for  P  -  J5) 

g 


(4.13) 
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* 


* 


¥ 


h 


and  I(o)q)  has  the  maximum  value 


> .  mii.  v  (2 


I(w. 


2  16 


4)  N(k,0) 


40 


(1+6)' 


(4.14) 


V.  Comparison  with  Experiment 

The  ratio  of  the  maximum  decay  line  to  the  maximum  of  the  purely 

2  2  2 

growing  line  for  a  homogeneous  plasma  with  A  y  /h^  -  and  3  -  1  Is  then 

(for  0  *  0.  we  take  y  -  w  ) 

ft  X  ft  ft 


I(k,(i)o)  (2. A)  N(k,0) 

I(k,uL(k))  "  A  NCk.^-u^Xk)) 


(5.1) 


The  ratio  of  the  Ion  fluctuation  factors  Is 

2 


N(k.O) 


N(k»aj0-wL(k)  ) 


(w^Ck))- 

+  2  2 
-  e  2k 


(5.2) 


2 

where  ■  ©j/m^.  At  the  spectral  peak  considered  here  u^-to^k)  *  C^. 
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Using  *  1.6(0e/m^)  for  0^  *  0^  we  have  (assuming^m  /m^  «  1) 


JLObPI 


0 

2T  (1.28) 
91 


N(k,0)Q-<i)L(k)) 


(5.3) 


-  3.6  for  0-0.  (0-1) 

ft  1 

Therefore  for  equal  electron  and  Ion  temperatures 


I(k,u>  ) 

-  ~  2.2  (5. A) 

Kk.Uj  (k)) 
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The  actual  backscattered  signal  observed  at  Areclbo  is  proportional 

12  3 

to  the  Integral  of  I  over  a  range  of  altitudes  *  *  .  In  the  case  of  the 
decay  line  we  find  from  (3.9)  that  I  is  proportional  to  )  where 

w^(k,2)  - 

so  that  a  given  frequency  arises  from  a  definite  altitude.  Assuming  that 

the  pump  amplitude  is  constant  over  the  altitude  range  in  question  the 

maximum  frequency  return  arises  for  the  altitude  z,  where  Au)  -  0 

d 


«0  "  “L<k.2d)+Wa(k) 


In  this  case  we  have 


(M2  A2(*)  m>2(*,) 

■^4"  -TZ-9- - ~ - - - N(k,<i)  (k))  — Z®- 

W(1-pd>  Vzd> 


(5.5) 


2  2  2 

where  ■  A  p  /A^  -  the  pump  ratio  relative  to  the  minimum  decay  instability 

threshold  and  where  we  have  taken 


d  (^(k.z) 
dz 


W 

2H 


where  H  is  the  density  gradient  scale  height. 

In  the  case  of  the  zero  frequency  mode  the  scattered  return  from  each 

altitude  is  centered  at  a  shifted  frequency  of  U)  .  The  most  sensitive  altitude 

o 

dependent  quantity  in  (4*12)  is 


BIO 


(5.6) 


m  -*l  (k,z)  w  (z  ) 

^(*)  ■  — v - "  -1  -  (z-z_)  8 


q  2HYt 


where  we  have  taken  z  to  be  the  altitude  of  best  matching  for  which 

ft 


Wk>i«)  ■  -v 


(5.7) 


The  altitude  integrated  spectrum  is  then 


,  (4tt)2A2  B2  (20>  ?.  [l^2U)]N(k.O) 

/  dz  1(0)  ,z)  -  2~~  16  ~2  2 -  dzT~^> - 7  (5.8) 

^  k^  (1+8)2  YL2(z  )  J  [l-H52(z)+26(z)P?J2 


(6<0) 


_  (4tt)2  A2  B2  H  (*o>  .  2H\ 

Kjj2  16  (1+B)2  Yl2(*g)  “p(Zo) 


(5.9) 


■  p  2  .  .  “ 
_ g _ + _ 1  nr/2  +  tan  *  ] 

a-*g2)  a-pg2)3/2  l  X7W_ 


where  P  -  A2jj2/A  2‘. 

8  g 

We  then  obtain  the  ratio  of  maximum  backs cat tering  cross  sections 


Hax(da/do))growint 
Max  (da  / da))  de  cay 


fiz  1(Vz)  «  B2 

/dz  Ko^.z)  V  A2(z.)  (1+B)2 


Kf‘*  ^ 


^tt/2  +  tan 


1  PR  |  N(k,0) 


fi-p  ; 

g 


(5.10) 
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With 


our  assumption  P  ■  %  for  the  observed  waves  and  assuming 


2  2 

A  (f  )  -  A  (z.)  this  ratio  turns  out  to  be  (for  $-1) 
8  d 


la) growing  =  2<16 


(5.11) 


Max  (  d<7  /  dw)  d  e  c  ay 


This  ratio  Is  larger  by  a  factor  of  2  or  3  than  those  observed  at  Areclbo. 
A  factor  of  2  or  so  could  arise  from  the  crudeness  of  our  approximations. 

In  particular  our  treatment  of  the  decay  instability  suffers  when  0g  ■  0^ 
because  the  Ion  acoust'c  mode  Is  poorly  defined.  The  theory  does  clearly 
Indicate  that  the  decay  line  peak  and  the  growing  line  peak  should  be  of 
comparable  order  of  magnitude  as  observed. 

In  the  calculation  above  we  assumed  that  the  pump  amplitude  was  the  same 
at  the  altitudes  z d  and  at  z  .  The  altitudes  are  defined  respectively  by 


“o^L  “L(k'V  \  -  *.  2H 


(5.12) 


V^a00  "  "  ,UL(k’0)  '  Zd  2ff 


(5.13) 


so  that 


2YT3 


As  “  Z  x  -  Z  ■  [Yt+W  (k)l  —  «  — — 
tf  g  L'L  a  J  u  fi) 


P  P 


(5.14) 


For  the  parameters  at  A^-ibo 


^  a1.2  z  10”3 


(5.15) 
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Thia  small  range  Justifies  our  approxi"u»tinn  of  using  constant  values  of 
in  our  calculations. 

The  spacing  of  the  Airy  function  maxima  of  the  pump  field  is  given  by 


As  ~  . 
“*  '  2 


(5.16) 


For  a  scale  height  of  H  -  100  km  and  c/u  ~  10m  this  is 


max  ~  4.10  ^ 

H  - 


(5.17) 


2  2 

Because  of  this  the  assumption  that  A  (z.)  *  A  (z  )  may  not  be  valid.  The 

d  g 

2  2 

ratio  is  reduced  if  we  assume  A  (z  )  <  A  (z.)  ,  i.e.  the  altitude  at  which  the 

g  d 

observed  growing  mode  is  excited  is  at  a  lower  point  on  the  Airy  pattern  than 

the  altitude  z^  at  which  the  observed  decay  mode  is  excited.  This  implies 

that  a  variability  in  this  ratio  should  be  observed  with  charging  scale  height 

H  and  changing  modifier  frequency  w  . 
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APPENDIX  C 


NONLINEAR  WAVE  OPTICS  OF  PARAMETRIC  PUMP  RADIATION 
IN  AN  INHOMOGENEOUS  PLASMA 


Donald  F.  DuBois t  Martin  V.  Goldman 
Dean  KcKinnls 


I.  Introduction 

Intense  electromagnetic  radiation  incident  on  a  plasma  can  excite 
parametrically  unstable  electrostatic  plasma  waves.1  This  flow  of 
energy  from  electromagnetic  to  electrostatic  waves  can  be  described  as 

2  3 

an  anomalous  absorption  of  the  electromagnetic  wave . '  Recent  calculations  * 
of  the  anomalous  absorptive  conductivity  (or  equivalently  the  nonlinear 

resistivity)  based  on  a  theory  of  saturation  of  the  electron-ion  decay 
instability  indicate  a  conductivity  which  is  linear  in  the  radiation  power 
when  the  latter  exceed  the  instability  threshold.  In  this  paper  we 
extend  this  calculation  and  find  that  when  the  radiation  intensity 
sufficiently  exceeds  threshold  the  nonlinear  conductivity  saturates 
at  a  constant  value.  This  model  of  the  nonlinear  conductivity  is 
added  to  the  linear  conductivity,  arising  from  electron-ion  collisions, 
in  the  wave  equation  for  the  radiation  as  it  propagates  into  a 

ty  _ 

plane-layered  inhomogeneous  plasma.  We  compute  the  effects  of  anomalous  absorp  - 
tion  on  the  reflectivity  and  the  spatial  distribution  of  the  radiation 
field,  flux,  energy  density  and  force  density.  The  nonlinear  resis¬ 
tivity  is  scaled  with  a  variable  multiplicative  parameter,  n,  to  enable 

4  5 

comparisons  with  various  more  detailed  theories.  *  This  form  of 
resistivity  may  also  be  of  some  interest  in  other  contexts.  For  example, 
numerical  simulations  in  other  parameter  regimes  (e.g. ,  unequal  temper¬ 
atures)  are  often  interpreted  using  similar  forms  for  the  resistivity, 
based  on  different  saturation  arguments. 

In  Section  2  we  review  the  results  of  the  theory  of  nonlinear  con¬ 
ductivity  plus  the  necessary  approximations.  The  geometry  and  para¬ 
meters  associated  with  the  radiation  wave  equations  are  set  forth. 


Cl 


In  Section  3,  some  of  the  approximations  of  previous  theory  are 
relaxed  and  the  nonlinear  resistivity  Is  generalized. 

Section  k  describes  the  method  of  Integration  of  the  radiation  wave 

equation. 

Section  5  Is  devoted  to  applications  of  the  theory  to  Ionospheric 
modification  experiments.7*8  We  conclude  that  anomalous  absorption  of 
-6  Wiz  radio  waves  with  vacuum  Intensity  -50  u  W/m2  at  the  F-layer 
of  the  ionosphere  (200-300km)  is  on  the  order  of  10?  -  20?. 
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II.  Review  of  'Theory 

It  has  been  shown  recently  * 3  how  to  relate  the  hl^- frequency 
electron  resistivity  to  the  spectrum  of  enhanced  Langnuir  waves ,  driven 
parametrically  unstable  by  Incident  radiation  and  saturated  by  induced 
scattering  off  ions.  The  conditions  for  validity  of  the  calculation 
were 

a)  The  electron  and  Ion  temperatures  do  not  differ  significantly 
(e.g.,  by  more  than  a  factor  of  two). 

b)  Ihe  total  Langmuir  wave  energy  density  must  be  less  than 
the  electron  particle  energy  density. 

c)  Ihe  radiation  field  must  sufficiently  exceed  the  linear 
parametric  decay  instability  threshold. 

Ihe  power  flows  from  the  radiation  into  the  Iangmuir  waves  and  finally 
into  the  particles.  The  first  condition  guarantees  that  the  acoustic 
waves  are  insignificantly  enhanced  relative  to  the  Langmuir  waves  and 
therefore  play  no  role  in  the  enhancement  of  the  resistivity.  The  higi- 
frequency  resistivity  formula  obtained  in  Reference  2  thus  differs  from 
the  usual  Dawson-Cberman  enhancement  formula.^  The  second  condition 
wa3  necessary  to  Justify  the  application  of  weak  turbulence  theory1*  to 
the  saturation  of  the  Langmuir  wave  intensity  spectrum.  The  third 
condition  guarantees  that  spontaneous  emission  plays  no  large  role  in 
the  energy  balance  considerations.  The  result  obtained  in  Reference  2 
for  the  dissipative  part  of  the  high-frequency  electronogpetic  (trans¬ 
verse)  conductivity,  a  !JL,  is 
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CD 


A 
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where  is  the  local  radiation  field,  I,  (k)  is  the  spectrum  of  electro- 
static  field  fluctuations  integrated  over  frequencies  near  the  Langnuir 
frequency,  and  y(k)  is  the  total  linear  damping  rate  of  electrostatic 
waves  of  wave  nuntoer  k. 

The  theory  developed  *  for  the  spectrum  1^  (k)  required  the 
constraints  a),  b),  and  c)  above;  in  addition,  neglected  the  zero-frequency 
(or  oscillating  two-streem)  instability  of  Nishikawa^;  and,  finally,  was 
valid  only  in  a  certain  region  of  k-space  to  be  described  below.  The 
spectrum  was  found  to  be  highly  peaked  in  angle  about  the  punp  polarization 
direction  far  wave  nurtoers  which  were  not  too  small.  As  a  function  of 
wave  nuntoer,  it  was  found  to  cut-off  sharply  for  wave  nuntoers  larger 
than  the  frequency-matched  wave  nuntoer,  k^,  determined  by 

o>  •  w  (1  +  3k2)*4  +  a  k  w  ,  (2) 

op  m  m  p 

where  u)Q  is  the  punp  frequency,  «p  the  plasma  frequency,  a  *  (n^A^), 
and  k  is  measured  in  units  of  the  Debye  wave  nuntoer.  The  analytic 
expression  found  in  Reference  2  for  the  angle-averaged  spectrum  was 


kc 

^<y-^rk  (PfMd7 

p 

Here  ft  is  the  electron  temperature  in  energy  units ,  d  -  2  x  10  a  , 

k  is  the  critical  wave  nuntoer  separating  the  decay  from  the  zero- 
c 

frequency  instabilities. 


i 
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•o  "  up(1  +  3kc>h 

and  f  Is  the  frequency-matching  resonance  function  given  In  Reference  2 
which  peaks  at  k^,  Is  zero  at  k^  and  tends  to  zero  as  K  decreases 
away  from  P  is  the  square  of  the  ratio  of  the  pump  field  to  the 
minimum  threshold  field  for  the  equal  temperature  three-mode  decay 
instability : 

P  '  Eo/e4  (5) 


4 

4ito0  id 

P 


(6) 


The  spectrum  given  in  Eq.  (3)  spreads  towards  lower  wave  numbers  with 
a  width  given  roughly  by 


kc  -  k 


,  2 


■In  •  3  Pa  £°r  kmln  >  0. 


(7) 


Poimila  (3)  for  the  spectrum  I  ±  (k)  is  valid  only  for  P  »  1,  for  k^  in 
the  collision-damped  regime  rather  than  the  Landau-damped  regime,  and 
for  wave  numbers  such  that  ok  >  y  (k).  When  this  last  condition  is 
violated,  a  second  (antistokes)  Langnuir  wave  upshift ed  frcm  the  punp 
frequency  is  dragged  in  and  increases  the  linear  decay  instability  threshold, 
thereby  reducing  the  spectral  intensity. 
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Under  the  above  conditions,  and  under  the  further  condition  that 

the  spectrum  not  he  spread  to  zero  in  wave  number  (i.e.,  2/3)Po  <  k  ), 

c 

expressions  for  the  total  energy  and  dissipative  nonlinear  conductivity 
were  obtained  in  Reference  2: 


E  = 


=  1  I 


u 


P2  n6 


(8) 


,  _  NL  .  1  . 

4»Ct  .-Py 


(9) 


In  the  course  of  the  present  paper,  it  will  be  found  necessary  to 

lift  all  of  the  above  restrictions,  with  the  exception  of  a),  b),  and  c) 

10a 

above  and  the  neglect  of  the  zero- frequency  instability. 

In  an  inhomogeneous  plasma,  we  cannot  regard  the  punp. field  as 
explicitly  known  locally  near  the  reflection  point  without  solving  Maxwell's 
equation  for  the  pump  field  and  including  any  nonlinear  contributions  to  the 
conductivity  such  as  that  of  Eqs.  (1)  or  (9).  If  we  adopt  the  plane  geometry 
illustrated  in  Fig.  1,  we  must  solve  the  following  second-order  nonlinear 
equation  for  the  (conplex)  punp  field  E  *  Eq  eiC^  in  the  presence  of  the 
absorptive  part  of  the  transverse  conductivity,  o^,: 


if?  +  r .  ^  +  (e„.»>  ij 

-2  2  u  I 

dx  a)©  o  I 


E  -  0 


(10) 


Here  the  spatial  coordinate  z  has  been  made  dimensionless  by  scaling  it 
in  units  of  c/g>0  ,  the  free-space  wavelength  over  2tt.  z  is  measured  from 
the  reflection  point  towards  lower  densities.  For  the  linear  density 
profile  we  are  considering  in  the  present  paper. 


C6 


(11) 


where  zq  is  a  dimensionless  density  scale  length,  related  to  the  customary 


scale  length  L  by 


(12) 


In  the  absence  of  any  dissipation  (a-T  -  0)  the  solution  to  Eq.  (10) 

which  is  damped  on  the  overdense  side  of  the  reflection  point  is  the  well- 

known  Airy  function,  sketched  in  Pig.  2.  Notice  the  first  maximum  occurs 
1/3 

at  z  *  zq  and  represents  an  enhancement  over  the  incident  power  by  a 

(swelling)  factor,  3.8  z^3.  A  value  of  zq  =  1000  corresponds,  for  the 

Nd- laser  value  of  c/L^  =  O.lji,  to  a  density  scale  length  L  »  100 p.  =  0.1  mm. 

For  conparlson,  in  radio-wave  ionospheric  experiments,  c/u>o  is 

~10  meters,  and  a  zq  of  1000  yields  L  «  10^  meters,  compared  to 

li  q 

measured  scale  lengths  of  10  -  10  meters. 

The  dissipative  conductivity  contained  in  Eq.  (10)  consists  of  the 
linear  conductivity  plus  a  nonlinear  contribution  suitably  evaluated  in  the 
inhomogeneous  plasma.  Perkins  and  Flick^  have  shown  that  the  linear 
stability  analysis  for  an  inhomogeneous  plasma  in  the  geometry  we  are 
considering  is  unchanged  from  the  analysis  of  an  infinite  homogeneous 
plasma  usually  considered,  provided  the  density  scale  length  exceeds  the 
electron  mean  free  path.  We  shall  consider  only  such  cases  here  and 
assume  that  o^,  can  be  obtained  from  the  Eqs.  (5)  (6)  and  (IS)  by  using  local 
values  of  n,  0,  and  Eq.  The  evaluation  of  at  a  given  position  in  the  inhomo¬ 

geneous  plasma  from  Eq.  (1)  requires  a  knowledge  of  the  angle-averaged  spectrum 
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(k)  for  the  appropriate  frequency-matched  wave  number,  k^.  Prom 
Eqs.  (2)  and  (11)  it  is  easy  to  see  that  when  1^  «  1,  it  is  (in  units  of 
the  Debye  wavenumber)  given  by  approximately 


f!3) 


Thus,  as  the  purrp  wave  propagates  inwards  towards  the  reflection  point, 
it  first  drives  Langmuir  waves  with  km  in  the  Landau-damped  regime, 
then  in  the  collision-damped  regime,  and  finally  it  encounters  the  range 
very  close  to  the  reflection  point  where  four-mode  effects  raise  the 
threshold  (k^  ^  Yc/ a ) .  In  order  to  evaluate  the  resistivity  in  all 
these  regions,  the  results  of  Reference  2  must  be  generalized  to  deal  with 
landau  damping,  four-mode  effects,  and  extreme  broadening  of  the 
spectrum.  These  effects  will  be  considered  next. 


HI.  Generalization  of  the  Nonlinear  Resistivity 

In  this  section  we  shall  discuss  the  qualititave  and  quantitative 
inplications  of  Landau  damping,  four-mode  effects,  and  spectrum 
k-broadening  for  the  nonlinear  resistivity  derived  in  Reference  2  and 
described  in  Section  II. 

As  we  shall  soon  see,  the  shape  of  the  punp  amplitude  squared 
obtained  from  solving  the  nonlinear  equation  (10)  will  not  differ  in 
quality  substantially  from  the  Airy  function  solution  illustrated  in  Fig.  2. 
That  is,  the  pump  amplitude  squared  will  generally  have  a  peak  at  z  * 
and  lower  peaks  near  3z^^,  5z^^,  etc.  Four-mode  effects  are  important 
when  ak  <C7+  7  ^(k)  /<up.  Whenever  they  are  important,  7^  is  reduced, 
so  we  may  qualitatively  look  for  the  solution  to 


C8 


(14) 


to  determine  the  wave  number  at  which  four-mode  effects  are 
inport  ant.  In  the  small  k  limit  y(k)  is  the  collision-damping  rate 
Yc,  and  independent  of  k,  so  one  solution  of  Eq.  (14)  is 


k  - 

4 


1 

a 


(15) 


We  shall  take  y  /w  to  be  of  the  order  4  x  10  for  the  ionospheric  F- 
c  p 

layer  at  ng  -  5  x  105  Cm  and  ©e  »  0.2  eV  (2319°k)  and  of  the  order 

~4 

3  x  10  for  a  laser-produced  plasma  at  n  s  10  and  ©  «  1  keV, 

6  6 


using 


12 


(k^/n)  in  0 


a) 


y  /a  * 
c  p 


(16) 


48 


For  the  ionosphere  the  ions  at  the  F-layer  are  mainly  oxygen,  so  a  =  5.83  x  10”^, 
whereas  for  laser  applications,  a  =  .0234  (for  hydrogen).  Four-mode  effects 
are  important,  therefore,  for  k  <  k^,  where,  for  the  two  cases, 

kIonosphere  .  x  10-H_  klaser  .  -0129  (17) 

If  we  use  Eq.  (13)  to  relate  this  frequency-matched  wave  number  to  a 
distance  from  the  reflection  point,  we  obtain 

^onosphere  „  ^  x  10\^er  .  5  x  ^  m 

For  a  typical  ionospheric  scale  height,  zQ  -  5  x  10 3 ,  whereas  for  a  typical 
laser-produced  plasma,  z  =10^,  so 


.Ionosphere  .  ?  x  1q-3j  .laser  =  >5 


(19) 
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which  are  both  on  the  low  pump-field  side  of  the  first  maximum.  The  situation 
illustrated  in  Pig.  3  •  In  this  region  the  nonlinear  conductivity  will 
usually  not  be  very  effective  in  the  nonlinear  Maxwell  equation  since  the 


punp  field  is  far  below  its  first  maximum.  Another  solution  to  Eq.  (14) 
is  in  the  Landau-damped  region  where  y(^)  becomes  large  as  increases. 
In  this  region, Eq.  (14)  becomes 


1  ,-3/2  ,-l/2kJ 
K 


(20) 


Which,  for  a  *  5.83  x  10”^  has  the  approximate  solution  kJonosphere  *  .240, 
and  for  a- =.0233  has  the  approximate  solution  k^61,  =  .265  .  Ibur-mode 
effects  therefore  again  begin  to  be  important  for  wave  numbers  exceeding 
.24  to  .265  in  the  Landau-damping  region.  Spatially,  this  wave  nunfcer 
corresponds  to  a  distance  from  the  reflection  point  found  from  Eq.  (13) 
to  be 


Z|j  =  210 

•3 

for  a  scale  length  of  zQ  =  10 J.  At  this  point  the  modified  Airy  pattern 
will  have  undergone  over  10  oscillations  and  decreased  in  amplitude  at 
the  peaks.  We  shall  soon  see  how  Landau  damping  reduces  nonlinear 
resistivity  in  this  region  even  more  drastically.  Thus ,  we  conclude 
that  four-mode  effects  are  of  marginal  importance. 

For  frequency -matched  wave  numbers  k^  in  the  collision-damped 
regime  and  sufficiently  large  pimping  power,  P,  the  expressions  in 
Eqs.  (8)  and  (9)  for  the  total  wave  energy  and  the  resistivity  should  be 
modified  to  be  consistent  with  a  spectrum  that  has  broadened  all  the  way 
down  to  zero  wave  number.  The  approximate  expression  for  the  angle- 
averaged  spectrum  I^k)  in  Eq.  (3)  is  still  valid.  From  Eq.  (7)  we  see 
that  for  punp  power  P  exceeding  P»,  where. 
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3k 
_ c 

2a 


(21) 


P* 


the  spectrum  will  be  totally  broadened.  The  total  energy  E  can  generally 
be  found  from  the  approximate  spectrum  of  Eq.  (3)  as  follows: 


k 

c 

dk  k 
k 

min 


Ix(k) 


48ir3d 


k 


min 


dk' 


(22) 


Integrating  by  parts  and  noticing  that  the  boundary  terms  are  negligible 
if  the  spectrum  is  small  at  kc  and  k||rin  or,  in  the  case  of  ®  0, 
may  even  diverge,  provided 


11m  k3!  (k)  -  0 
k-*0  -  1 


Then 


E 


48rr3d 


(23) 


Using  the  peakedness  of  f(k)  at  k^,  and  the  approximate  result  *bkat  the 
integral  of  f(k)  ■  2a  /3,  we  obtain 


-  ki® 

E  :  3 

48ir  d 


k 

m 


-  (k.-k 

I  ,  a  -  k  (k  -  k  .  )  +  — V-®10 
3  c  c  min  2 


>■1 

J 


(2*0 


It  is  easy  to  show  xkat  kQ  -  k^  *»  a/6.  The  exact  expression  for  d  ls^a 
d  *  (a  /**32ir  )  (k^/n)  (u^/y).  We  thus  obtain  from  Eq.  (2*0  two  expressions 
for  the  energy  E.  For  the  case  k^  >0,  kp  -  k^^  ,  p(2/3)a  and  the  ex¬ 


pression  in  Eq.  (8)  results  when  P  »  1.  For  the  case  k  *  0,  or 

(2/3  Pa  >  k  ,  we  obt*iin 
c 
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_  ,  k 

E  '  n©  —  X-  — 
*  u  a 
P 


,  k 

for  P>  ^  — 
2  a 


(25) 


For  sufficiently  large  P,  therefore  the  total  energy  goes  as  P  rather 
2 

than  P  and  the  corresponding  conductivity  is  eventually  independent  of 
the  punp  power.  From  Eqs.  (1)  arid  (25),  in  the  collision-damped  region 

for  V 


tooj1  .  i«S£L 
I  .2 

o 


,v  K  3k 

=  2L  _E  /i _ E\ 

2 it  a  u  4aP'  * 


p  >  2  »  i 

—  2  a 


Another  effect  not  included  in  the  above  analysis  but  of  possible  importance 
at  low  wave  numbers  and  high  P  is  the  inception  of  mode-coupling  of 
Langnuir  and  acoustic  waves  as  a  stronger  saturation  mechanism  than 


ir.duced  scattering  of  Langmuir  waves  into  Langmuir  waves.  ■  Mode¬ 


coupling  is  discussed  in  Reference  4a,  Eqs.  (25)  through  (28),  where  it  is 
concluded  that  it  is  relatively  weak  except  at  small  wave  numbers.  When 
it  dominates,  the  spectrum  can  be  shown  to  be  rigorously  isotropic  in 
angle  as  a  consequence  of  the  balance  between  spontaneous  emission  and 
mode-coupling  damping.  The  spectrum  can  then  be  shown  to  be  about  one- 
third  as  wide  in  k  as  the  narrow  spectrum  discussed  above  in  the  absence 
of  mode-coupling  and  the  total  wave  energy  is  doubled.  The  ratio  of  the 
Langmuir  wave  nonlinear  mode-coupling  to  induced  scattering  damping 
rates  at  k^  can  be  shown  from  the  spectral  form  to  be 


r!c  I  p  y 

ysp'  5  7  «D 

m  v 


(27) 


When  this  ratio  is  large  compared  to  one,  the  total  wave  energy  and, 
hence,  the  conductivity  is  multiplied  by  a  factor  2  in  comparison  with  its 
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value  In  the  absence  of  mode-coupling.  Since  the  isotropic  spectrum  does 

not  spread  to  zero  wave  number  until  P  is  three  times  larger  tlian  before  , 

the  total  wave  energy  remains  proportional  to  P  for  higher  values 

of  P  than  previously.  All  of  these  features  increase  the  resistivity  for 

NC  SP 

those  values  of  P  and  for  which  y  /y  is  large  compared  to  one. 

Next,  we  explore  the  effects  of  linear  landau  damping  on  the 
spectrum  and  resistivity  for  the  case  of  k^  exceeding  about  0.2.  A 
ccqplete  treatment  of  saturation  in  the  landau-damped  region  requires 
the  addition  of  particle  orbit  perturbation  corrections  as  well  as  the 
induced-scattering  saturation  mechanisms  we  have  included.  If  we  ignore 
the  orbit  perturbation  corrections  (which  may  be  weaker) ,  then  a  linear 
Iandau-danping  term  must  be  added  to  y^1  in  the  equation  (4)  of  Reference 
2.  Ibis  leads  to  an  approximation  for  the  angle-average  spectrum  which 
is  of  the  following  form: 


6  fkc 

V10  ■  6di?  /dk'  <pcf(k,)  -  i  -  r(k-) 

_  k 

^  • 


(28) 


P  in  this  expression  is  defined  in  accordance  with  Eqs.  (5)  and  (6), 
c 

but  with  y  in  Eq.  (6)  redefined  to  be  the  collision  damping  rate  yc 
cited  in  Eq.  (16).  T(k)  is  the  ratio  of  the  landau  damping  rate  to  the 
collision-damping  rate, 


r(k)  i  'f  /5f  3/2  i,  e-  1/2  k2 


(29) 
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The  integral  over  r(k* )  is  trivial,  leading  to 


*i  00  •  i k 


[/r 


(Pcf(k')  -  i)  _  k3  r(k  )  +  k3  r(k) 


The  width  of  the  spectrum  is  now  given  by 


(30) 


kc  "  kmin  =  Pc  3  °  "  kc  r(kc>  +  *Ln  r(kmin> 


(31) 


Unlike  the  result  in  the  collision-damped  regime,  this  width  depends  on 

k  and  is  narrower  than  before.  The  pump  power  P  at  which  the  spectrum 

c 

is  Just  spread  from  k  to  k  ■  0  is  given  by  the  solution  to  the 

transcental  equation  (31)  with  k^  =  0.  For  a  -  1/^2  and  ycA>p  =  3  x  10“* 

some  solutions  are  given  in  Table  2.  From  the  values  of  T(k  )  we  see 

c 

that  Landau  damping  is  just  equal  to  collision  damping  for  k  ~  -0.215. 

This  result  is  not  changed  very  much  for  the  ionosphere.  For  snaller 

kc,  Pc  remains  about  14.  In  Figs.  H  through  0,  the  spectrum  of  Eq.  (30) 

is  plotted  numerically  and  integrated  to  obtain  the  total  energy, 

for  a  =  1/H2  and  y  /u  =  3  x  10“\  In  Figs  H  and  5,  P  is  kept  fixed 

at  10  and  the  matched  wave  number  is  increased  from  .21  where  Landau 

dairying  is  Just  becoming  Important  to  .25  where  it  is  very  strong.  In 

agreement  with  Table  1,  the  spectrum  has  not  spread  to  zero  for  P  =10. 

c 

As  k^  increases,  the  spectrum  is  seen  to  narrow  and  diminish  in  inten¬ 
sity  until  it  is  flat  at  zero  for  k^  =  0.25  (not  shown).  In  spatial 
tenia  the  spectrum  diminuation  occurs  over  one  or  two  cycles  of  the 
modified  Aliy  pattern  in  the  region  where  the  field  peaks  are  already 
diminished.  In  Pigs.  6  through  8,  P  is  fixed  at  20  and  the  natched 

w 

rave  number  is  again  increased  from  .21  to  .25.  At  k^  =  0.21  and  0.23 


1 
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the  spectrum  has  broadened  to  zero,  in  agreement  with  Table  1,  whereas 
fork^-0  •25  the  spectrum  has  narrowed  and  diminished  in  intensity 
considerably  due  to  heavy  Landau  danping.  The  cross-marks  in  Pigs.  6  and 

p_ 

7  are  values  of  k  I^(k),  which  is  seen  to  go  to  zero  as  k-K).  Thus, 

there  is  no  large  contribution  to  the  total  energy  from  small  wave 

numbers  (where  four-mode  effects  should  really  be  included  anyway) . 

In  Table  2  the  total  energy  is  integrated  numerically  from  the  spectra 

of  Pigs.  4  through  8  and  also  for  the  case  P  *  40.  Except  for  P  values 

c  c 

of  the  order  of  40  or  more,  the  total  energy  of  waves  decreases  sharply 
in  the  Landau-damped  region.  For  P  values  less  than  -40,  the  contributions 

C 

to  the  resistivity  in  the  Landau-damped  region  are  also  decreasing  with 

larger  (farther  from  the  reflection  point)  and  we  shall  adopt  the 

point  of  view  that  a  crude  Interpolation  of  the  resistivity  from  the 

collision-damped  Vc  lue  to  zero  deep  in  the  Landau-damped  region  is  adequate 

at  these  P„  values, 
c 

To  conclude  this  discussion  of  generalization  of  the  nonlinear  re¬ 
sistivity  we  next  examine  modifications  due  to  the  new  findings  of  more 
exact  theories 4,5,5a.  Numerical  integrations  of  the  time-dependent 
1-dimensionalized  nonlinear  plasma  wave  kinetic  equations4,5  yield  total 
plasma  wave  energies  from  2  to  T  times  larger  than  the  values  obtained 
fban  the  approximate  analytic  theories  (e.g.,  eqn.  8).  2-dimensional 
integrations  have  also  been  performed  recently511  but  the  resulting  total 
energy  has  not  yet  been  given.  In  all  of  these  integrations  the  form  used  for 
the  resonant  mismatch  factor  f(k)  also  can  vary  the  coefficient  of  the  total 
energy  by  a  factor  of  2  or  3*  DuBois  and  Goldman  use  a  form  of  f(k)  based 
on  a  resonance  approximation  for  the  ion-acoustic  dielectric  function 
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"hich,  however,  is  truncated  to  zero  beyond  k  ,  the  critical  wave  number 

V 

separating  the  decay  from  the  zero-frequency  instability  (k  is  defined  in 

0 

equation  4  of  the  present  paper).  In  references  3,5  and  5a,  the  exact 
equal-temperat ure  ion-acoustic  dielectric  function  is  used,  but  no  truncation 
at  kc  is  included.  The  treatment  of  Reference  3  appears  to  yield  a  total  energy 
about  9 A  times  larger  than  that  of  Reference  4,  due  to  the  different 
choice  for  f(k). 

In  view  of  the  fact  that  the  total  wave  energy  (and  hence  the  nonlinear 
absorptive  conductivity  which  is  proportional  to  the  total  wave  energy)  are 
known  only  to  within  a  constant  which  may  vary  by  as  much  as  an  order  of 
magnitude,  we  shall  use  a  variable  multiplicative  parameter  n  definsd  in  terms 
of  the  total  plasma  wave  energy  E  as  follows: 

E.  Unions,  r<|£” 


E 


,  k 


(32) 


Comparing  with  eqns.  (8)  and  (25)  we  see  that  n  =  for  the  theory  upon  which 
those  equations  rest.  In  table  3  we  indicate  values  of  n  which  may 
arise  from  various  theories. 
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In  accordance  with  the  above  discussion,  the  absorptive  conductivity 
we  shall  use  in  the  integrations  of  the  radiation  wave  equation  (10)  is 
2  .. 


4ira„ 


u 


u 


2  Y„ 


*•*  o  th  o 


(33) 


The  first  term  on  the  right  side  is  the  linear  collisional  damping  tenn. 

vei  is  the  electron- ion  collision  frequency,  which  we  take  to  be  2y  where 

0 

Yc  Is  defined  in  equation  (16).  The  second  term  on  the  right  side  is  the 
nonlinear  contribution  to  the  absorptive  conductivity  which  includes  four- 
mode  effects,  an  interpolation  from  the  collision-damped  region  into  the  heavily 
Landau-damped  region,  and  broad-spectrum  effects,  as  well  as  the  scaling  factor 
n  defined  in  equations  (32)  and  Table  3  ■  The  various  quantities  appearing 
in  o^are  defined  below: 


?  Y(kJ 

E7,  *  64ttti(z)0 - — 

th  fa)_ 


*  +  4^ 
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u 


a 


(34) 


Here  y(^)  is  the  total  damping  rate.  Including  collision-dairping  plus  Larriau- 
danplng  when  significant.  YO^)  assures  that  the  conductivity  decreases  in 
the  landau-damped  region.  The  factor  in  square  brackets  in  eqn  (34)  guar¬ 
antees  the  increase  in  threshold  associated  with  the  onset  of  four-mode 
effects,  although  it  is  quantitatively  correct  only  for  electrons  not  compared 
to  ions  and  then  more  precisely  should  Involve  y1^,  as  described  above  in 
eqn  (14).  The  reduction  of  the  conductivity  in  the  landau-damped  regime  is 
likewise  only  a  quantitative  interpolation. 

The  factor  W(E)  accounts  for  the  broad  spectrum  effects  described 
by  equations  (8),  (25),  and  (32): 
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where  P  is  defined  in  Eq.  (5). 

IV.  fritegration  of  the  Nonlinear  V/ave  Equation  for  the  Pump  Field  in  a  Linear 
Density  Profile 

Here  our  object  is  to  integrate  the  wave  equation  (10)  for  the 
radiation  field  in  a  linear  density  profile,  building  into  the  nonlinear 
conductivity  the  effects  described  by  equation  (33)  of  the  last  section. 

Even  in  the  presence  of  linear  and  nonlinear  damping  the  wave 
equation  generally  has  a  solution  which  grows  in  the  overdense  region  in  the 
direction  of  higher  densities  and  one  which  decays.  We  begin  our  numerical 
integration  with  an  arbitrary  boundary  value  of  E  deep  in  the  overdense  region 
and  integrate  towards  lower  densities.  In  this  direction  the  unphysical  and 
therefore  unwanted  "growing  solution"  decays  spatially  and  we  have  verified  it 
is  negligible  by  the  time  it  reaches  the  oscillatory  region  (z  >  0). 

At  z  =  zq  the  density  vanishes  and  E  consists  of  an  incident  and 
reflected  plane  wave: 

E  :  a  e1(vz)  ♦  b  e-«Vz)  =  E^,  M  +  (36) 

Since  we  are  assuming  time-dependence  as  e-*wt  and  decreasing  z  is  increasing 
density,  |a|  is  the  amplitude  of  the  incident  wave  and  |b|  is  the  amplitude 
of  the  reflected  wave.  Ihe  reflection  coefficient 


is  then  extracted  at  zQ  along  with  the  incident  intensity  at  zQ,  both  corresponding 
to  the  chosen  boundary  value  for  E  in  the  deep  overdense  region.  Let  us  separate 
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(38) 


the  field  E  Into  Its  real  and  Imaginary  parts, 

E"Ereal(z)  +  1  W(z) 

Hie  wave  equation  (10)  with  the  appropriate  is  then  Integrated  from  the 

deep  ovendense  boundary  value  at  z^^  <  0  towards  positive  z  by  a 

Runge-Kutta  method  with  error  control.  When  Ereal  and  near  zq  are 

known  f!ran  the  numerical  integration  which  started  at  z.  ,  <0  the 

initial  * 

reflection  coefficient  at  z  =  zq  is  determined  as  follows:  Take  the  derivative 
of  Eqs.  (36)  and  (38) 


3Et 

3z| 


i  a  +  i  b  ■  E  ,+E. 

real  imag. 


(39) 


Equations  (36)  -  (39)  are  now  algebraic  equations  which  may  be  solved  for  a 
and  b  in  terms  of  Ereal(zo),  Elmag>(zo)  and  their  first  derivatives  at  zq. 
Die  reflection  coefficient  is  then: 


R  ■  |— 
1  a 


If  °W<zo>  +  <z.»2  * 

|_*Ereal^zo^  E  lmag/zo"  +  ^Elmag/zo')  +  E'rea]/zo^  _ 


(40) 


To  understand  the  radiation  field  patterns  produced  by  this  nonlinear 
resistivity,  we  first  examine  the  zero-resistivity  (a  =  0)  Airy  function 

pattern.  As  a  function  of  position  and  time,  the  electric  field  may  be 
written  as 

E(z)  ■  |A(z)  I  cos ($  (z)  +  ut)  , 

where  A(z)  is  an  Airy  function  solution.  The  Airy  function  solution  A(z) 
normally  oscillates  both  positive  and  negative,  and  the  phase  <P  is  constant, 
so  E(z)  is  a  standing  wave.  We  have  instead  chosen  to  use  the  absolute  value 
|A(z)|  and  represent  the  negative  and  positive  character  of  A  by  allowing  the 
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phase  <J>  to  change  discontlnuously  by  tt  at  the  nodes  of  A(z)  for  reasons  to 
become  evident.  Equation  (4l)  of  course  still  represents  a  standing  wave  and 
with  the  choice  of  zero  phase  at  z  =  0  is  identical  with  A(z)  cos  wt.  The 
phase  <J>( z )  and  amplitude  squared  A2  are  plotted  in  Fig.  9>  with  A  normalized 

to  an  incident  wave  of  unit  amplitude. 

In  the  following  sections  we  consider  ionospheric  and  laser  cases  of 
nonlinear  absorption.  For  each  set  of  values  of  n,the  incident  intensity, 

t  inc 

Iinc=  Bn  * 

and  the  scale  length  z  ,  we  have  obtained  a  value  for  the  reflectivity  and  plots 

c 

of  the  radiation  amplitude  |A(z)|2  (in  units  of  the  incident  intensity)  and 
the  phase  <Kz),  defined  by  equation  (4l).  The  amplitude  maxima  will  be  seen 
to  exhibit  a  reduced  swelling  factor,  due  to  the  absorption,  but  otherwise 
show  the  characteristic  oscillatory  structure  of  Figure  9.  The  reduction  of 
the  relative  intensity  at  the  maxima  is  an  example  of  the  saturation  mechanism 
we  have  called  punp  depletion,  and  which  is  the  analog  of  bump-on-tail  flat¬ 
tening  in  quasilinear  theories,  such  as  the  gentle  burp  plasma  instability. 

That  is,  the  unstable  plasma  waves  are  reacting  back  on  their  source  of  energy 
(in  this  case  the  pimp),  tending  (through  absorption)  to  reduce  its  effect¬ 
iveness  in  driving  the  waves  unstable.  Also,  the  amplitude  minima  lift  from 
zero  and  the  phase  in  their  vicinity  will  show  smoothing  from  the  step-like 
behavior  of  Figure  9,  so  the  wave  begins  to  show  sane  of  the  character  of 
a  travelling  wave  rather  than  a  standing  wave.  (Note  a  perfect  travelling  plane 
wave  would  have  a  linear  phase  and  constant  amplitude.) 

These  effects  are  more  pronounced  closer  to  the  incidence  point  zQ 
rather  than  the  reflection  point  z  =  0.  The  reason  is  as  follows: 
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One  half  of  the  absorption  of  a  wave  reflecting  off  an  overdense  plasma 

occurs  on  the  return  trip.  Even  In  a  case  of  total  absorption,  the  wave 

patterns  very  near  the  reflection  point  would  not  look  like  a  travelling  wave 

because  the  reflected  wave  has  not  damped  to  zero  until  it  returns  all  the  way 

back  to  zQ.  Thus,  even  for  cases  of  strong  linear  or  nonlinear  absorption 

the  oscillatory  structure  of  the  radiation  wave  pattern  near  the  reflection  point 

is  preserved.  It  is  only  after  several  oscillations  (at  least  1255  zQ)  away 

fron  the  reflection  point  that  the  nonlinear  absorption  will  have  had  its  full 

14 

effect.  Geometric  optics  approximations  which  only  contain  the  envelope 
of  the  wave  pattern  are  inherently  inaccurate.  Although  they  can  be  scaled 
to  give  reasonably  good  reflect i\ ities,  they  cannot  show  the  large  variations  in 
pump  field  and  consequently  in  plasma  wave  energy  near  the  reflection  point. 

This  structure  is  of  interest  in  the  ionospheric  and  laser  applications,  in  the 
foiraer  case,  because  of  scattering  measurements  from  specific  altitudes  near  the 
reflection  height,  and  in  the  latter  case  because  of  possible  pondermotive- 
f  ore e- induced  motions  of  the  density  profile. 

These  quantities  associated  with  the  solution  to  the  nonlinear  pump 
radiation  wave  equation  which  are  of  interest  to  us  are  the  Poynting  flux, 
the  energy  density,  and  the  pondermotive  force.  In  our  geometry  and  units, 
with  a  plane  polarized  pump  field,  the  time-average  Poynting  flux  is 

2TT/0) 

Jz  ‘  W  2?  \  dt  I?  (z>t)  x  S'2'*01 


0 


C 

5? 


p2  14 

%  3z 


(*2) 


where  E(z)  =  E  e'^,  and  z  is  in  units  of  c/u>,  as  before.  By  plotting  -J 
o  z 

we  can  get  seme  idea  of  where  the  absorption  is  occurring,  since  the  loss 
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tenn  in  the  conservation  equation  is  proportional  to  -3J  /3  .  With  this 

z  z 

definition  of  J_  as  the  total  Poynting  flux  it  must  vanish  for  the  case  of 

Z  - 

total  reflection.  This  is  evident  from  equation  (42)  and  Figure  9,  since 
the  Airy  function  solution  has  steps  in  <j>  only  at  the  nodes  of  Eq. 

Ihe  time-average  energy  density  of  the  radiation  in  the  plasma  is 

given  by 


This  expression  includes  magnetic,  electric  and  coherent  particle  kinetic 
15 

energy  densities  . 
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V.  Nonlinear  Absorption  of  Ionosphere-Modifying  Radio  Waves 

8 

In  recent  experiments  low  frequency  (5-10  Mhz)  ordinary  radio  waves 
with  local  vacuum  intensities  of  up  to  50  Moj/m^  are  reflected  from  the  F- layer 
of  the  ionosphere  (200-300  km  altitudes).  The  electric  field  of  these  waves 
near  the  reflection  point  should  be  approximately  plane-polarized  along  the 
direction  of  the  geomagnetic  field  and  is  believed  at  its  peak  to  exceed  the 
threshold  for  excitation  of  parametric  instability.^  Ionospheric  measure¬ 
ments  show  the  reflectivity  of  lower  intensity  diagnostic  radiation  in  the 
same  frequency  range  drops  within  seconds  when  the  high-intensity  radiation 
pump  is  turned  on  and  cones  back  up  again  within  seconds  when  the  pump  is 
turned  off.  Since  the  growth  time  for  the  decay  parametric  instability  is 
on  the  order  of  tens  of  milliseconds  this  behavior  may  be  due  to  scattering 
ft*an  associated  slow  hydrodynamic  instabilities  rather  than  from  anomalous 
absorption  or  scattering  off  enhanced  Langrmir  or  ion— acoustic  waves.  However, 
the  issue  is  still  not  resolved. 

Incoherent  back-scattering  experiments  performed  at  Arecibo,1^  using 

430  Mhz  radio  waves,  have  probed  approximately  vertically  propagating 

4irf 

Langmuir  and  ion-acoustic  waves  of  wave  number  k  =  — —=■  ~  .05  kp  (for 

"  0.2  eV) .  According  to  equation  (13) »  this  corresponds  to  scattering 
from  an  altitude  z/zQ  =  7.5  x  10  neglecting  geomagnetic  field  effects. 

At  present,  there  is  some  question  as  to  whether  the  scattering  is  from 
nonlinear ly  saturated  plasma  waves  or  from  the  pump  wave  beating  with  ion- 
acoustic  fluctuations  to  produce  Langnuir  fluctuations  below  the  parametric 
instability  threshold. 

Against  the  background  of  these  experiments  it  is  especially  important 
to  determine  the  extent  of  pump  nonlinear  absorption  and  the  spatial  behavior 
of  the  punp  electric  field.  The  existence  of  significant  nonlinear 
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ab  orption  bears  on  reflectivity  experiments  for  the  pump  or  a  secondary 
diagnostic  wave,  as  well  as  the  entire  question  of  anomalous  heating  and 
possibly  airglow  production.  The  local  values  of  the  pump  electric  field 
and  Its  modified  swelling  pattern  are  needed  to  understand  the  scattering 
experiments. 

The  first  set  of  computer  runs  have  been  performed  using  the  following 
set  of  parameters: 

nQ  ■  5  x  105  cm- ^ 

0 *  0.2  eV  «  2319°K  (daytime) 

I  m 

a  -  e  *00583  (oxygen) 

0>o/27T  *  modifier  (pump)  frequency  =  6.35  Mhz 

Derived  quantities  based  on  the  above  parameters  are: 
w  =  x  10'  radlans/sec 


Veio  *  ^co  *  350 


4' 


th 

4Tm  0 
o  e 


=  16  -^°  =  7  x  10"5 


*th  *  §?  Eo  =  i^W/m2  (locally) 
th 


(*»5> 


The  peak  plasma  frequency,  w  is  equal  to  the  pump  frequency,  y  is 

corputed  using  equation  (16).  The  threshold  for  the  parametric  instability 

cE2 

is  indicated  as  Ith,  and  should  be  understood  as  the  value  which  -g^- 

(including  swelling  due  to  inhcmogeneity)  must  exceed  locally  to  excite  the 

parametric  decay  instability.  We  consider  two  scale  heights:  H  =  37.5  km 
to  H  j. 

(z^  =  -£  H  =  5000)  and  H  =  75  km  (zn  =  =  10  ). 

The  integrated  total  absorption  is  numerically  measured  by  the 
quantity  1  -  R,  where  R  =  |  /E^c  I  ^  is  the  power  reflection  coeffi¬ 
cient  at  the  "source  point"  z  =  zq.  We  also  compute  for  each  case  the 
integrated  purely  linear  reflectivity  obtained  by  integrating  the 

wave  equation  (10)  using  only  the  linear  dissipative  conductivity  (first 
term  on  right  side  of  equation  (33)).  These  results  compare  favorably 
(withir  50$)  with  the  value  of  calculated  from  geometric  optics. 

-8  Vel  z 

_  3U  0 

Rlln  =  e  0  '' 

The  quantity  Rn  -  R  =  (1  -  R)  -  _  R^^)  is  a  measure  of  the  anomalous 

absorption  (excess  beyond  linear). 

In  Figure  10  both  1  -  R  (left  axis)  and  -  R  (right  axis)  are 

2 

shown  for  an  incident  vacuum  intensity  of  50uw/m  ,  which  is  about  the 

Q 

maximum  presently  available  at  the  Platteville  transmitter  in  ( olorado. 
(Approximately  one  half  this  mlue  is  attainable  at  Arecibo,  Puerto  Rico.) 
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The  absorption  is  plotted  as  a  function  of  the  multiplicative  scaling 
parameter,  n,  defined  in  equation  (33)  and  Table  3.  As  explained  in 
Table  3,  the  results  of  the  best  available  theories  indicate  n-values 
between  h  and  3/*J-  'IVie  recent  results  of  Fejer  and  Kuo^a  indicate  a 
saturated  spectrum  of  plasma  waves  which  has  seme  angular  extent,  although 
they  do  not  compute  the  total  plasma  wave  energy  or  conductivity.  The 

ii  iio 

isotropic  angle  spectrum  approximation  of  Goldman  and  DuBois  *  predicts 
a  rough  doubling  of  the  total  plasma  wave  energy  as  compared  to  the  extreme 
anisotropic  spectrum  solution,  so  we  also  look  at  n  values  as  high  as  2. 

For  a  scale  height  of  37*5  km,  the  linear  absorption  1  -  R  =  9%  for 
the  parameters  of  (44)  and  (45).  (Equation  (46)  gives  1255.)  The  excess 
nonlinear  (anomalous)  absorption  FL^  -  R  is  between  955  and  23/5  depending 
on  n.  The  behavior  of  -  R  with  n  is  roughly  linear. 

For  a  scale  height  of  75  km,  the  linear  absorption  is  1  -  R^n  =  17$ 
for  the  parameters  of  (44)  and  (45).  (Equation  (46)  gives  2355.)  The 
excess  nonlinear  absorption  R|ir|  -  R  is  between  14%  and  33%,  depending  on 
n.  Longer  density  scale  lengths  are  therefore  seen  to  lead  to  stronger 
absorption.  The  physical  reason  is  that  the  absorption  region  is 
"stretched",  although  from  the  above  results  it  appears  that  this  is  a 
stronger  effect  for  the  linear  absorption  than  the  nonlinear  contribution 
(presumably  because  the  nonlinear  absorption  region  is  thinner  spatially) . 

The  above  results  are  for  relatively  high  daytime  electron  tempera¬ 
tures.  Strictly  speaking,  the  effects  of  photoelectrons  on  the  plasra 
wave  damping  should  be  included  for  daytime  parameters.  We  have  crudely 
estimated  the  effects  of  photoelectrons  in  a  separate  data  run  by  adjusting 
the  Landau  damping  to  be  three  times  larger  than  the  collisional  damping 
of  the  Langnuir  waves  at  a  matched  wave  number  of  km  =  O.lkp.  (Ordinarily 
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the  Landau  damping  is  much  smaller  than  the  collision  damping. )  The  factor 
3  was  chosen  as  an  average.  At  this  density  and  temperature  downward  propa¬ 
gating  plasma  waves  have  a  photoelectron  damping  rate  equal  to  the  collisional 
damping  rate  and  upward  propagating  plasma  waves  have  photoelectron  damping 
five  times  collisional  damping.  (The  photoelectrons  originate  mainly  at  lower 
altitudes.)  Also,  there  should  not  be  many  photoelectrcns  at  energies 

beyond  ra  /Jd.2  ~  20eV.  This  has  the  effect  of  raising  the  threshold  for 
2  'k  ' 

parametric  instability  by  a  factor  of  ~4  at  altitudes  in  the  vicinity 

of  2/zq  =  .03  and  by  more  below.  It  does  not  alter  the  nonlinear 
conductivity  at  higher  altitudes  because  the  phase  velocities  at  the 

correspondingly  smaller  matched  wave  numbers  are  too  high  to  allow  wave 
resonance  with  photoelectrons.  The  net  effect  is  to  condense  the  spatial 
region  of  nonlinear  absorption.  For  all  parameters  the  same  as  in 
Figure  IQ,  but  photoelectrons  included  in  the  above  manner,  and  a  scale 
height  of  75  km,  we  find  for  n  =  2/3,  FL^-R*  11$  anomalous  absorption 
compared  to  14$  anomalous  absorption  without  photoelectrons;  and  for 
n  *  2,  R]1n  -  R  =  27$  anomalous  absorption,  compared  to  33$  anomalous 
absorption, without  photoelectrons.  Photoelectrons  thus  appear  to  reduce 
the  percent  of  anomalous  absorption  by  a  factor  of  about  20$.  When  photo¬ 
electron  damping  of  Langmuir  waves  is  a  strong  effect,  or,  indeed  even  If 

the  normal  Landau  damping  region  is  to  be  treated  properly,  the  absorptive 

13 

conductivity  arising  from  orbit  perturbation  saturation  must  be  used 
rather  than  the  conductivity  associated  with  weak  turbulence  saturation, 
which  we  have  been  using  here.  We  are  presently  putting  into  effect  a 
program  for  doing  this. 

We  have  also  made  some  runs  at  lower  electron  temperature  ( .  13eV  = 

1500°K).  Since  the  linear  absorptive  conductivity  goes  as  and  hence 
-3/2 

as  0  ,  the  linear  absorption  should  and  did  increase.  The  nonlinear 

w 
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absorption  also  Increased,  although  not  as  strongly  as  the  linear  absoiption 
Making  sane  numerical  comparisons  for  the  H  =  75  km  shown  in  Figure  IQ  but 
with  all  other  parameters  except  temperature  the  same,  we  find  the  linear 
absorption  1  -  R^  Increased  from  11%  to  28%  as  Ge  Is  decreased  front 
2320  K  to  1500°K,  whereas  the  anomalous  absorption  R^  -  R  Increased 
fran  HIS  ( .2eV)  to  19%  (.13eV)  for  n  =  2/3  and  from  33%  (.2eV)  to  39% 
(.13eV).  Nighttime  anomalous  absorption  would  thus  appear  to  be  stronger 
than  daytime  anomalous  absorption,  although  orbit  perturbation  effects'^3 
must  be  included  to  really  determine  whether  this  is  true.  Daytime  neutral 
collisional  absorption  of  the  modifier  (pimp)  at  lower  altitudes  (D-region) 
would  make  daytime  anomalous  absorption  still  weaker. 

In  Figures  11  and  12  we  have  plotted  the  absorption  1  -  R  as  a 
function  of  incident  intensity  for  the  two  scale  heights  H  *  37.5  km  and 
75  km  using  the  same  n- values  and  other,  fixed  parameters  as  in  Figure  10. 
Ihe  absorption  goes  up  somewhat  more  slowly  than  linear. 

In  Figures  13  and  14  we  plot  the  pump  radiation  electric  field 

amplitude  squared  versus  altitude.  Ihe  punp  amplitude  is  normalized  to 

the  incident  amplitude  and  the  altitude  is  in  units  of  c/w  .  For  the 

o 

parameters  of  (44)  and  (45)  this  distance  is  7.5  meters.  The  spatial  wave 
patterns  are  plotted  for  an  incident  vacuum  intensity  of  50  yw/m^,  a  scale 
height  of  75  km  (zq  =  10  ),  and  a  nonlinear  scaling  parameter  n  =  2. 

The  density  is  5  x  105  cm"3  and  the  temperature  is  0.2eV.  Results  are 
shown  both  for  a  purely  linear  absorptive  conductivity  and  for  the  linear 
plus  nonlinear  absorptive  conductivity  of  equation  (33).  The  peak  of  the 
Airy  function  swelling  pattern  is  indicated  by  X  for  comparison.  Ihe  effect 
of  pump  depletion  discussed  in  section  IV  Is  evident.  Ihe  peak  swelling 
factor  at  the  first  maxirum  is  reduced  from  82  with  no  absorption  to  71 
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with  linear  absorption,  to  56  with  linear  plus  nonlinear  absorption.  The 
amplitudes  at  all  the  maxima  are  reduced  by  somewhat  smaller  percentages. 

Far  the  case  shown  we  are  not  very  far  above  threshold  even  at  the  first 
maxima  of  the  field  pattern.  Based  on  the  actual  swelling  at  the  first 
maximum  of  the  curve  which  incudes  linear  plus  nonlinear  absorption  we 
find  P  =  17.  If  we  had  used  the  Airy  function  swelling  factor  this 
number  would  have  ccme  out  to  be  P  -  24. 

For  comparisons  with  scattering  experiments,  z  must  be  multiplied  by 
7.5  meters  -to  convert  it  into  distance  below  the  reflection  point.  If  the 
scattering  experiment  measures  a  wave  number  k^/k^  =  .05,  this  would  corre¬ 
spond  to  z  =  75  in  Figure  13 .  Since  the  density  scale  height  is  never  known 
precisely  and  since  we  have  left  out  effects  of  the  geomagnetic  field,  this 
information  is  not  enough  to  determine  whether  the  pump  field  is  at  a  trough 

or  a  peak  for  the  measured  plasma  waves.  (Also,  Incident  intensities  of  order 

2 

25hW/m  should  be  used  for  Arecibo  scattering  experiments.)  Pimp  depletion  is 
seen  to  reduce  the  magnitude  of  the  peaks  from  the  values  they  would  other¬ 
wise  have  had  from  linear  swelling  and  linear  absorption  alone.  Note,  the 
pump  field  amplitude  still  shows  the  nodes  characteristic  of  a  standing  wave 
pattern,  in  accordance  with  the  discussion  at  the  end  of  Section  IV,  due  to 
reflection  over  short  distances  from  the  reflection  point  z  =  o. 

Figure  14  shows  the  pump  field  amplitude  squared  down  at  lower 
altitudes.  The  node-lifting  associated  with  the  increased  absorption  and 
partial  conversion  from  a  standing  wave  to  a  travelling  wave  is  apparent 
here. 

Figure  15  shows  the  pump  phase  versus  distance  for  the  same  para¬ 
meters  and  distance  range  of  Figure  14.  The  phase  $  is  defined  by 
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E  ■  Eq(z)  cos  (<Kz)  +  w t),  where  the  amplitude  Eq(z)  Is  defined  to  be  positive. 
Figure  8  shews  the  phase  in  the  absence  of  absorption  for  zQ  =  1C^.  The 
smoothing  of  the  corners  of  the  step  like  phase  is  also  indicative  of  the 
transition  from  a  standing  wave  to  a  travelling  wave. 

In  Figure  16  we  have  plotted  the  negative  of  the  total  Poynting  flux 

in  units  of  the  incident  flux  for  the  same  parameters  as  in  Figure  13. 

The  total  Poynting  flux  is  defined  in  equation  (42).  In  the  absence  of  any 

absorption  this  net  flux  must  vanish.  The  lower  curve  in  Figure  16  includes 

only  linear  absorption.  The  upper  curve  includes  both  linear  and  nonlinear 

absorption.  The  steeper  regions  of  the  oscillating  flux  are  the  regions 

where  the  absorption  is  strongest,  and  correspond  to  the  maxima  of  Eq(z) . 

The  curve  including  both  linear  and  nonlinear  absorption  shows  stronger 

oscillations  than  the  curve  including  only  linear  absorption  because 

3J/3z  «  ciNlj  E  “  EH.  Also,  the  curve  including  both  linear  and  nonlinear 
o  o 

absorption  rises  to  a  higher  value  at  z  =  zq  than  the  curve  with  linear 
absorption  alone  because  the  smaller  the  reflected  wave  at  that  point,  the 
higher  the  net  Poynting  flux. 

Figure  T7  is  another  plot  of  the  pump  radiation  electric  field  amplitude 
squared  versus  distance  z,  with  all  parameters  the  same  as  in  Figure  13, 
except  for  a  value  of  the  nonlinear  scaling  parameter  n  *  2/3  (compared  to 
n  *  2  in  Figure  13).  Again,  the  points  around  z  =  50  to  100  are  where  the 
incoherent  scatter  return  fretn  a  430  Mhz  diagnostic  might  be  expected  to 
originate  since  these  altitudes  correspond  to  the  frequency  matched  wave 
number. 
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IV .  Conclus  ion 


Ihe  nonlinear  transverse  absorptive  conductivity  which  is 
associated  with  the  steady  state  flow  of  energy  from  an  electromagnetic 
pump  rave  to  parametrically  excited  electrostatic  Langnuir  waves  has  been 
generalized  beyond  previous  1-dimensional  saturation  calculations.  New 
effects  which  have  been  built  into  o^  include: 

a)  Effects  of  a  broad  spectrum  in  |k|  at  high  pump  powers.  This 

changes  the  dependence  of  of  from  linear  in  the  punp  power  to  independent 
of  the  pump  power. 

b)  Pour  mode  effects  which  decrease  for  very  small  and  for 
^  in  the  Landau-damped  region. 

O  Qualitative  decrease  of  o^  when  is  in  the  Landau-damped  region. 

d)  Qualitative  effects  of  more  exact  1-D  and  2-D  theories  are  included 
through  a  multiplicative  scaling  parameter  n,  to  which  o^  is  proportional. 

Ihe  resulting  expression  for  o^  is  added  to  the  linear  conductivity  in 
the  wave  equation  for  the  pump  radiation  as  it  propagates  towards  its 
reflection  point  in  a  plane-layered  inhomogeneous  plasma.  The  pump  wave 
equation  has  been  integrated  to  find  the  total  absorption  (or  reflectivity) 
at  various  incident  power  levels  and  frequencies,  and  various  density  scale 
lengths,  temperatures,  and  n-values.  In  addition,  the  spatial  pattern  of 
the  pump  amplitude  and  phase  has  been  displayed,  am  the  local  total  Poynting 
flux  hao  been  plotted  as  a  function  of  position  for  the  various  runs.  The 
method  used  does  not  rely  in  any  essential  way  on  a  linear  density  gradient 
and  can  be  generalized  in  future  computations. 

For  the  ionospheric  cases  run,  the  effects  (a)  and  (b)  above  do  not 
matter  in  of  because  the  power  levels  are  so  low  (typically  no  more  than 
20  times  the  parametric  instability  threshold  power) .  Anomalous  absorption 
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Is  defined  as  total  absorption  minus  linear  absorption.  For  the  current 

maximum  available  vacuum  pump  power  of  50  nw/m2  at  200-300  km  in  the  5-10 

Mhz  radio-wave  range  in  no  case  do  we  find  anomalous  absorption  of  more  than 

40$.  Without  the  effects  of  orbit  perturbation  built  into  nighttime 

anomalous  absorption  is  stronger  than  daytime  anomalous  absorption.  Typical 

daytime  anomalous  absorption  (with  photoelectron  effects  built  in  very  crudely 

as  a  limitation  on  the  spatial  extent  of  the  nonlinear  region)  is  10-25$ 

2 

at  50  uw/m  .  At  night  it  is  20-40$.  These  values  correspond  to  a  scale 
height  of  75  km.  The  absorption  which  occurs  at  higher  pump  powers  is 
computed  and  may  be  of  interest  if  self- focusing  effects10  can  intensify  the 
local  pump  power. 

An  interesting  feature  of  the  pump  radiation  field  as  a  function  of 

position  is  that  it  is  essentially  a  standing  wave  at  altitudes  within  a 

few  percent  of  the  density  scale  height  below  the  reflection  point  even 

when  there  is  strong  anomalous  absorption.  This  kind  of  behavior  appears  to 

17 

have  been  bom  out  by  incoherent  scatter  measurements.  We  find  it  leads 

to  strong  oscillations  in  the  Poynting  flux  pattern  and  hence  to  bards  of 

altitudes  at  which  local  anomalous  heating  would  be  strongest.  In  general, 

the  maxima  of  the  pump  amplitude  oscillations  are  lower  than  they  would  be 

with  no  absorption  or  with  linear  absorption  alone.  This  phenomenon  of 
i| 

pimp  depletion  is  evidence  of  the  transfer  of  energy  from  the  pump  to  the 
plasma  waves. 
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Table  1. 

Pump  Power  Ratio  Pc  for  Spreading  Spectrum  to  k  =  0  as  a 

Function  of  k  .  (a  =  1 /42  and  V  /co  =3xl0"4). 
c  c  p 


k 

c 

r(kc) 

Pc  for  spreading  to  zero  k 

.21 

.  70 

14 

.  23 

2:  S3 

16 

.25 

8.  95 

25 
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Table  3 


Values  of  the  scaling  parameter  n  appearing  In  the  total  plasma  wave  energy 
E  or,  equivalently,  the  absorptive  nonlinear  conductivity  ol^1. 


n 


Theory 


Resonant  approx,  for  f(k),  cutoff  at  k  . 

c 

of  Ref.  2,  iJa 


Resonant  approx,  for  f(k)  no  cutoff  at  kQ. 
of  Ref.  3 


Analytic  saturated  theory 


Analytic  saturated  theory 


Resonant  approx,  for  f(k),  cutoff  at  k  .  Time-dependent  numerical 

c  h 

Integration  (B.  Godfrey,  D.P.  DuBois,  M.V.  Goldman  ) 


a 

x-  Exact  f(k).  No  cutoff  at  k  .  Time-dependent  1-D  numerical  inte- 
gration.  w.  Kruer,  E.  Valeo' 


?  Exact  f(k) .  No  cutoff  at  k  .  Time-dependent  2-D  numerical  inte¬ 
gration.  J.  Fejer,  Y.Y.  Kuo^ 


Value  used  by  F.  Perkins  and  by  G.  Meltz  and  N.  Tbmljairovich  in 
recent  geometric  optics  approximations111  to  the  nonlinear  absorption 
problem. 


2 


Figure  Captions 


Figure  1 

Geometry  for  normal  Incidence  of  electromagnetic  wave  on  a  plane-layered 
inhomogeneous  plasma.  The  spatial  coordinate  z  is  in  the  dimensionless 
units  of  c/u>o>  the  free-space  wavelength  over  2tt. 

Figure  2 

Airy  function  solution  (zero  resistivity)  for  field  amplitude  squared  in 

units  of  incident  field  amplitude  squared.  zq  is  the  dimensionless 

density  scale  length  (in  units  of  c/u  ). 

o 

Figure  3 

Frequency-matched  wave  number  km  (in  units  of  kp)  as  a  function  of  spatial 
position  z  (in  units  of  wq/c),  illustrating  the  regions  of  low  and  high 
resistivity  corresponding,  respectively,  to  low  and  high  Iangnuir-wave 
intensity. 


Figure 

Angle-averaged  intensity  of  the  Langrulr  spectrum  versus  wave  number  for 
Pc  -  10  and  matched  wave  number  k^  =  0.21.  Wave  number  is  in  units  of  the 
Debye  wave  number. 

Figure  5 

Angle-averaged  intensity  of  the  Langjrulr  spectrum  versus  wave  number  for 

P  ■  10  and  matched  wave  number  k  =0.23.  Wave  number  is  in  units  of  the 
v  m 

Debye  wave  number. 
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Angle-averaged  Intensity  of  the  Langnuir  spectrum  versus  wave  number  for 
Pc  *  20  and  matched  wave  number  k^  *  0.21.  Wave  number  is  in  units  of 
the  Debye  wave  number. 


Figure  7 

Angle-averaged  intensity  of  the  Langnuir  spectrum  versus  wave  number  for 
Pc  *  20  and  matched  wave  number  k^  *  0.23.  Wave  number  is  in  units  of 
the  Debye  wave  number. 


Angle-averaged  intensity  of  the  Langnuir  spectrum  versus  wave  number  for 

*  20  and  mat  died  wave  number  k  =0.25.  Wave  number  is  in -units  of  the 
c  m 

Debye  wave  number. 


Amplitude  squared  and  phase  of  the  Airy  function  (zero  resistivity) 
solution,  E(z)  =  |A(z)|  cos  ($(z)  +  off).  The  sign  of  the  Airy  function 
A(z)  is  represented  as  m  phase  changes  at  nodes  of  A(z).  The  amolitude 

•3 

is  normalized  to  incident  intensity  unity;  zq  =  10  . 
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Figure  10 


Integrated  absorption  1  -  R,  R  -  lE^/E^J2,  for  Incident  Intensity 
^Inc  *  cElncy/®1T  3  5&iW/m  02%)  and  various  values  of  the  nonlinear 
scaling  parameter  n  defined  In  equation  (33)  and  Table  3.  Both  linear 
and  nonlinear  absorption  are  Incited  and  the  following  parameters  are 
used:  n  a  5  x  105  cm'\  0e  -  0.2eV  (2319°X>,  f^,^  =  6.35  mhz, 
vei  ~  ^50  IIz*  ^ie  *ons  a:re  assumed  to  be  oxygen.  Scale  heights  of 
37.5  km  (zq  =  5000)  and  75  km  1zq  =  104)  are  used.  The  rigit  side  scale 

on  each  plot  Is  P  -  which  Is  a  measure  of  the  relative  percent 

of  nonlinear  absorption. 


Figure  11 

Integrated  linear  and  nonlinear  absorption,  1  -  R.  R  =  Ie  I2/If  I2 

1  refl1  X|  lnc1 

as  a  function  of  Incident  Intensity,  1^  =  cE^c/8tt.  The  nonlinear 

scaling  parameter  n  Is  defined  in  equation  (33)  and  Table  3.  The  scale 

height  is  H  37.5  km  and  the  integrated  linear  absorption  is  shown  as 
well. 


Figure  12 

Integrated  linear  and  nonlinear  absorption,  1  -  R,  R  =  |E  ^J2/jjr  |2 
as  a  function  of  Incident  Intensity,  1^  -  cE^/Stt.  The  nonlinear 
scaling  parameter  n  is  defined  In  equation  (33)  and  Table  3.  The  scale 

height  is*  H  =  75  km  and  the  Integrated  linear  absorption  is  shown  as 
well. 
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Figure  13 


Pimp  radiation  electric  field  amplitude  squared  (In  units  of  Incident 
amplitude  squared)  versus  altitude  (In  units  of  c/u  =  7.5  meters). 

Lower  altitudes  are  towards  the  right  and  the  lowest  altitude  shown  is 
*1  km  below  the  reflection  height.  Linear  and  linear  plus  nonlinear  cases 

p 

are  plotted  for  an  Incident  vacuum  Intensity  of  50uW/m  ,  a  scale  height 
of  H  =  75  km  (zq  =  10^),  and  a  nonlinear  scaling  parameter  n  =  2  (see 
equation  (33)  and  Table  3).  The  density  is  5  x  lO^cm” ^  and  the  tempera¬ 
ture  is  0.2eV.  The  peak  of  the  Airy  function  swelling  factor  is  shown 
by  X  for  comparison.  The  integrated  linear  absorption  is  1  -  R  =  17?, 
whereas  the  integrated  linear  and  nonlinear  absorption  is  an  additional 
33?. 


Figure  14 

Radiation  amplitude  squared  versus  distance  for  the  same  parameters  as  in 
Figure  12,  but  at  lower  altitudes  (in  the  vicinity  of  3  km  below  the 
reflection  point).  Node  lifting  is  beginning  to  become  apparent  here. 


Figure  15 

Radiation  phase  versus  distance  for  the  same  parameters  and  distance 
range  of  Figure  13.  The  phase,  <f>,  is  defined  by  E  =  Eq(z)  cos  (<J>(z)  +  wt), 
inhere  the  amplitude  Eq  is  defined  to  be  positive.  Figure  8  shows  the 

•3 

phase  in  the  absence  of  absorption  for  zq  =  10  . 
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Figure  16 


Negative  of  Poynting  flux  In  units  of  incident  flux  for  parameters  of 
Figure  12. 


Figure  17 

Radiation  amplitude  squared  (in  units  of  Incident  amplitude  squared) 

versus  distance  (in  units  of  c/w  «  7.5  meters).  The  incident  Intensity 
2 

is  50UwAi  ,  the  scale  height  is  75  km,  and  tlie  nonlinear  scaling  para¬ 
meter  n  *  2/3  (see  equation  (33)  and  Table  3).  Tbe  density  is  5  x  10^ 
cmi”^  and  the  tenperature  is  0.2  eV.  The  peak  of  the  Airy  function 
swelling  factor  is  shewn  by  X  for  comparison. 
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IM4-1 


ANISOTROPIC  ANGLE-AVERAGED  INTENSITY 


ANISOTROPIC  ANGLE'/VERAGEO  INTENSITY 


Pc  =  20.00  KM  =  .23 


FIGURE  7 
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FIGURE  9 
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INTEGRATED  LINEAR  +  NONLINEAR  ABSORPTION, 
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|Eo/Eincl 


SCALE  HEIGHT  =  75  km  (Z0=I04)  77  =  2 

v  AIRY  SWELLING  FACTOR 
80 1~  *  NO  ABSORPTION 


DISTANCE  Z  IN  UNITS  OF  c/« 


FIGURE  13 
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lEo'ElncI 


SCALE  HEIGHT  =  75km  {Z0=I04!  17  =  2 


400  405  410  415  420  425 

DISTANCE  Z  IN  UNITS  OF  c/cu 


FIGURE  14 
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RADIATION  PHASE,  *(z),  IN  RADIANS 


i 


SCALE  HEIGHT  =  75km  <Z0  =  l04)i7  =  2 


400  405  410  415  420  425 

DISTANCE  Z  IN  UNITS  OF  c/w 


FIGURE  15 
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NEGATIVE  POYNTING  FLUX,- J/J|„c 


SCALE  HEIGHT  s  75  km  (Z0*IC4)  77=2 


DISTANCE  Z  IN  UNITS  OF  c/cu 


FIGURE  16 
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SCALE  HEIGHT  =  75 km  (Z0  =  I04)  77  =  2/3 
X  AIRY  SWELLING  FACTOR 


DISTANCE. Z  IN  UNITS  OF  c/<u 


FIGURE  17 
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APPEND IX  D 


ELECTRON  ACCEI2RATION 


Jerome  Weinstock  and  Bandel  Bezzerides 


I  I  W1.4 


HECTRON  ACCE1ERATI0H: 


V 


The  electron  velocity  distribution  function  is  eubstentlally  changed  during 
modification.  This  change  la  evidenced  by  large  enhancements  of  natural 
63OQA  (red)  and  5500^  (green)  airglov  observed  over  Boulder  and  63OQA  observed 
over  Arecibo.  These  alterations  of  the  electron  velocity  distribution  can  have 
important  consequences.  For  examples ,  the  formation  of  energetic  tails  may 

(1)  create  new  plasma  wave  instabilities  such  as  the  ion-acoustic  instability, 

(2)  change  the  saturated  spectrum  of  short  scale  turbulence,  and  (3)  change  the 
radar  scattering  cross-section  of  the  turbulence. 

We  have  developed  a  theory  for  quantitatively  predicting  the  change  in  the 
velocity  distribution  of  electrons .  Such  changes  are  primarily  due  to  accelera¬ 
tion  of  electrons  by  short  wavelength  plasma  waves  propagating  in  the  Landau 
damped  region,  where  typical  wavelengths  are  about  15  cm  (about  5  Debye  lengths). 
Such  short  scale  waves  are  excited  in  the  bottom  of  the  modified  region  as 
indicated  in  Figure  1. 

The  electron  velocity  distribution  function  f  »  f(v,z)  satisfies  the 
diffusion  equation 


df  .  df  e 

St  +  51’5 


(1) 


where  v  is  velocity,  is  the  pump  field,  D(v,z)  is  the  diffusion  tensor. 

Horizontal  stratification  is  assumed  so  that  f  and  D  depend  upon  height  z. 

'f  we  average-  (l)  over  a  small  enough  vertical  range  A  within  the  heated 

z 

region  such  that 


1 


r 


df 

57 


(2) 


D1 


then  (l)  becomes 


df  d  9f  ' 

5t  5v  ‘  I  2  ’  57 


(3) 


where  the  superbar  denotes  an  average  over  the  vertical  range  A^,  and  we  have 
neglected  the  small  effect  of  Eq  and  of  collisions  upon  f . 


RED  AIR  GLOW 

The  airglow  measurements  provide  a  direct  measurement  of  the  enhanced 
photon  flux  of  atomic  excitations.  This  photon  flux  is,  in  turn,  a  measure 
of  the  enhanced  flux  of  electrons  causing  the  excitations.  The  flux,  and 
lumber  density,  of  exciting  electrons  is  determined  by  the  quantity 

^ dv  o(v)  f(v) 

where  o(v)  is  the  excitation  cross-section  as  a  function  of  v.  We  can  solve 
(l)  for  d/dt  Jdv  ct(v)  f(v)  and  then  determine  the  enhanced  flux  of  exciting 


electrons.  We  multiply  both  sides  of  (l)  by  a(v)  and  integrate  over  v  to 

obtain 


(«0 

(5) 


The  diffusion  tensor  D  is  given  by 

09 


2  ?  r 

D-(-’z)“?  Reai/^ 


dt  exp  [i(o)-k *t)  t  -  k  :DtJ  ] 


(6) 


o 

and  has  been  calculated  in  detail  for  modification  conditions.  The  component 
of  D  along  the  pump  direction  is  described  in  Fig.  2. 


D2 


Note  that  D(v)  is  fairly  narrowly  peaked  about  the  wave  phase  velocity.  (This 
phase  velocity  is  about *5  vfi  in  the  Landau  damped  regime.)  Hence,  the  main 
contribution  of1  the  v  integral  in  the  right  side  of  (5)  comes  from  v  »  w/k 


st  =  -  faz  i 


Substituting  (6),  for  D,  into  (7)  we  have 


si  ct(v)  f(v) 


ri  -  2  t  dk  C  C  k*v 

-&k  ^ 


x  k  .  exp[-i(to  -  k*v)t  -k2:  D  tj 


Ae  now  note  that  velocity  integral  on  the  right  side  of  (8)  is  proportional  to 

NL 

the  imaginary  part  of  the  nonlinear  susceptibility  -  x  (k,  to)  so  that,  with 
k*v  ra  u>,  (8)  can  be  written 


ct(v)  f(v) 


H  /■%  *t>  *L.V  , 

fJ  Ji&ip  ~*r 


Im  X  (k,u>)  w2(yJL  +  y^)A> 


where  y  is  the  nonlinear  damping  rate  of  Lan^nuir  waves  and  y.  is  linear 
L  X 


Landau  damping.  Furthermore,  at  saturation 


NL  (linear  growth  rate) 
\  "  YL 


Yl(P  f(k)  cos^e  -  1) 


where  P  is  the  pump  ratio,  y^  is  the  total  linear  damping  (“^  =  v  +  Y^) 
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collisions!  +  Landau  damping,  f(k)  is  the  width  function,  and  9  is  the  angle 
between  pump  and  k.  Substituting  (10)  and  (ll)  into  (9)  we  have 


3 

St 


J dv  o(v) 


<4)  2  . 

-j~  Pf(k)  cos  9  -  v] 


(12) 


Using  the  saturation  spectrum  for  (E^)  given  in  this  report  for  the  Landau 
damped  region  we  finally  have  for  (12) 


3 

St 


/ 


dv  o(v)  f(v) 


(13) 


where  <(6E)2)  is  the  total  energy  density  of  Langmuir  waves,  and  nQ  is  the 
ambient  electron  density. 

The  enhanced  number  density  of  exciting  electrons  is  essentially  determined 
by  (13).  Thus  if  we  let  n*(6300)  denote  the  enhanced  density  of  electrons 
exciting  63OQ&  airglow  then  we  can  write 


3 

St 


J dv  cr(v)f  =  <<a(v)>) 


3 

St 


n*(630o) 


(14) 


where  ((a(v)))  is  defined  as  the  average  cross-section  per  exciting  electron. 
But,  from  the  right  side  of  (7),  together  with  (6),  we  see  that  the  main 
contribution  to  3f(v)/3t  comes  from  v  wuu/k.  Hence  we  approximately  have 


from  (l4) 


3 

St 


J dv  a(v)  f(v)  «  a(^) 


3n*^o300) 

3t 


Substituting  (15>  in  (13)  yields 


3n»(6300) 

St 


<(6E)g> 

4rfnkr 


(15) 


(16) 


hut,  for  energies  less  than  4  ev,  n(v)  is  approximately  given  by 
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J  * 


(17) 


a(v)  a  vc 


so  that  (16)  becomes 


9n*(6300) 

5t 


<(SE)S> 


4m  kT 
o 


(18) 


This  describes  the  rate  at  which  exciting  electrons  are  produced.  The  saturation 
energy  density  <(5E)2)/(4mokT)  is  given  in  this  report  as  a  function  of  P  and 
(kve/io),  where  k  is  the  matching  k.  For  a  typical  modification  over  Boulder 
P  «  3  at  the  Landau  damped  height  for  which  k/^  «5  so  that  <(6E)2)/(4mokT)  »  10"^, 
and  v .  =  v  so  that  vT  «  1000.  Hence 

C  li 

*  5  x  10k  cn-3  Bee-1  (19) 


FLUX  OF  EXCITING  ELECTRONS 


The  flux  of  exciting  electrons  accelerated  down  the  geomagnetic  field, 
in  the  Landau  damped  region,  is  obtained  by  integrating  dn*/dt  over  the  height 
range  of  the  Landau  damped  region  for  which  P»3.  This  height  range  is  about 


0.3  k 


m 


Hence 


FLUX 


6.3  k 

m 

f 

m 


Bn* (6 300) 
dt 


■  0.3  x  10^  cm  x  ggffi P°). 


(20) 


Consequently,  f^r  typical  modification  conditions  the  flux  of  exciting  electrons 
is  with  (20)  and  (19), 
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(21) 


FLUX  «  10^  electrons  -cw"^  -sec"^ 

To  compare  this  flux  with  the  experimental  photon  flux  various  loss 
mechanism  of  exciting  electron  must  be  taken  into  account  such  as  collision 
with  nitrogen.  When  this  is  done  we  find  that  the  flux  predicted  in  (21 ) 
can  explain  the  observed  50  Rayleigh  enhancements  of  the  63OQA  airglow  over 
Boulder. 

In  deriving  (l 3 )  and  (l8)  -  (21)  we  have  only  had  to  calculate  an  average 
of  the  modified  distribution  function  f(v).  The  task  that  now  remains  is  to 
solve  (l)  for  f(v)  at  all  v. 
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APPENDIX  E: 


LINEAR  THEORY  OF  DIFFUSION  MODES 


Martin  V.  Goldman,  Donald  F.  DuBois, 
Richard  Berger,  and  Edward  Williams 


University  of  Colorado 


We  are  concerned  .here  with  the  linear  properties  of  first  order  collisional- 
induced  density  diffusion  across  magnetic  field  lines  In  a  homogeneous  plasma. 
Section  I  gives  the  conditions  for  the  simplest  theory  which  ignores  angular 
dependence  and  thermal  fluctuations.  Section  II  examines  the  diffusion  mode 
for  potential  relevance  to  ionospheric  modification  experiments.  Sections  III 
and  IV  relax  some  of  the  approximations  in  I. 

The  elements  and  assunptions  which  go  into  the  simplest  linear  fluid  theory 
of  the  diffusion  mode  are  listed  below. 

1.  Electric  fields  produced  by  the  diffusion  are  assumed  electrostatic 
(only  Poisson's  equation  is  necessary): 

2.  The  pressure  is  given  by  an  isothermal  equation  of  state,  p  «  N0. 

The  condition  for  this  is  |  ne<u  «  K^k2,  where  Kex  -  5neR2Vel  iS  the  trans¬ 
verse  electron  thermal  conductivity,  R  is  the  electron  g/roradius,  and  u 
and  k  are  the  frequency  and  wavenumber  of  the  diffusion  mode.  This  condition 
is  only  marginally  satisfied  in  the  simple  theory  presented  here,  although 

the  theory  is  easily  generalized  in  Section  IV  to  include  the  thermal  equations. 

3.  Viscosity  terms  are  neglected.  This  neglect  is  justifiable  provided 
the  wavelength  of  the  diffusion  mode  is  sufficiently  long,  or  kR^« 

where  R±  is  the  ion  gyroradius.  A  more  general  theory,  in  which  viscosities 
are  retained,  has  been  presented  by  Dawson  and  Okuda1  who  are  concerned  with 
diffusion  phenomena  associated  with  ion-ion  collisions  rather  than  electron- 
ion  collisions. 

1|.  Electron-ion  momentum  transfer  terms  are  retained  in  the  electron 
momentum  equations.  These  terms  and  electron  cyclotron  motion  terms  dominate 
the  inertial  terms  (|u|  «v  ,  «  u  ). 
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5.  The  Ion  charge  density  response,  to  the  electrostatic,  field  kE^ 
is  the  scatic  Debye-screened  response. 


JLBU 
k  h 


where  k^  is  the  ion  Debye  wavenumber. 

6.  In  the  linear  theory  we  may  assure  quasineutrelity:  ng  ~  n^.  Off 
resonance,  the  linear  dielectric  function  for  the  diffusion  mode  is  given  by 


e  *‘Xe  +  -55—  »  where  Xe  is  the  electrostatic  susceptibility  defined  by 
k 

n,pe  * 

and  to  be  evaluated  below. 

7.  Diffusion  along  the  magnetic  field  lines  occurs  on  a  much  faster 
time  scale  than  the  diffusion  across  the  field  lines  we  are  considering  here 
(satisfied  when  u>ce  »  vel). 

The  most  important  property  of  the  diffusion  mode  for  our  purposes  is  its 
extreme  aspect  sensitivity.  This  diffusion  mode  exists  only  for  wave  vectors 
k  which  are  extremely  close  to  perpendicular  to  the  magnetic  field  lines . 

In  Section  in  of  this  paper  we  show  that  the  mode  changes  from  diffusion-like 
to  wavelike  over  an  angle  to  the  perpendicular  to  given  by 

a  -  Vel  k 
“pe  "m 

Far  the  parameters  over  Boulder  this  is  of  the  order  of  10"^  radians. 

8.  For  the  present,  assume  ^  in  the  z-direction  and  first  order  spatial 
variation  of  all  quantities  is  in  the  y-direction.  With  these  approximations 
the  first  order  electron  fluid  equations  are  given  by 

<«  £+no37V° 


0  -  -  0jr  en0  (El?  +  |  x  B,)  -  „el  vn, 


i  m 
o  e 


£2 


Combining  these  gives,  to  within  the  above  approximations ,  for  the  fir.*;  order 
electron  density,  n. 


(ID 


an  u  r2  a2n  Veieno  3  _  „  n 

St  -  Yei^  7“ ; 2  “  — 2“  ST  h  °* 

3  y  mewce 


Fourier  transforming  equation  (11),  we  find  fbr  the  electron  susceptibility 
defined  in  (6), 


(12) 


X^oO 


lv 


R?k2 


ill 

ei  e 


k2  to  +  iv  .R2k2 
el  £ 


2  2 

where  R  is  the  electron  gyro-radius,  defined  by  R  =  0  /m  u>  .  r/fe  shall 
o  c  e  c  cg 

use  X  in  the  next  paper.  Setting  equal  to  zero  the  dielectric  function  e 
e 

(defined  above  in  (6))  gives  the  dispersion  relation 
(13)  w  -  -  iDk2, 

where 

(»)  ds  velR^(i  +  ei/ee). 

Equation  (13)  is  characteristic  of  a  diffusion  process  in  which  D  is  the  diffusion 
coefficient.  To  see  the  process  more  clearly  we  can  bypass  equation  (12) 
through  (1*0  and  put  directly  into  equation  (11)  the  static  ion  response  n^  to 
E  given  in  (5)»  and  the  quasineutrality  condition  ng  s  n^.  This  gives  the 
"diffusion  equation, 

(15)  |r  -  =  0. 

3t}  3y2 


Figure  1  sketches  out  how  an  initial  density  disturbance  uniform  axung  the 
z-axis  but  confined  to  y  =  0,  diffuses  out  in  y  as  a  function  of  time. 
Figure  2  shows  the  related  particle  motions. 
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II.  Potential  Relevance  to  Central-Line  Phenomenology 

The  diffusion  mode  is  of  interest  in  connection  with  experimental  results 

because  it  has  an  aspect  sensitivity  confining  this  mode  to  less  than  about  0.1°. 

At  present  we  have  no  plasma  theory  to  account  for  the  observed  results.  However, 

it  is  interesting  to  note  the  following: 

Hie  157  Mhz  central-line  data  results  probe  a  several-Hz-wide 
field  aligned  disturbance  of  wave  number 


(16)  k  *  .066  cnfl  (X  *  95  cm) 

Inserting  this  into  equation  (13)  gives 

(17)  |w|  -  15  Hz 

where  we  have  used  the  following  parameters: 
400  Hz 


ei 


0.2  eV 


(18) 


<0  V, 

e  i 

B  *  0.5  gauss 


w  «  8.8  x  10  rad/sec 
ce 


R  =  2.1  cm 
e 

Thus,  aspect  sensitivity  and  observed  frequencies  suggest  the  possible 
participation  of  the  diffusion  mode  in  the  central  line  scattering,  although'-' 
there  is  at  present  no  theory  which  predicts  the  strong  intensities  at  the 
probed  wavelengths  or  the  dependence  at  these  wavelengths  on  the  modifier 
Intensity.  In  the  next  paper  we  describe  how  the  modifier  can  drive  unstable 
diffusion  modes  at  wavelengths  between  50  and  500  meters  with  low  thresholds 
(close  to  that  of  the  decay  instability). 
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(22) 


HI.  Behavior  of  Diffusion  Mode  for  k.  »  Bo  *  0. 

In  order  to  determine  the  behavior  of  the  diffusion  mode  for  wave  vectors 
k  at  an  angle  to  Bo,  we  include  spatial  variation  of  first  order  quantities  in 
both  the  y-  and  z-  direction.  For  the  present,  we  do  not  assume  that  the  ions 
are  static  or  that  |u>|  «  vg^.  With  these  exceptions  we  retain  the  ass unpt ions 

in  Section  I.  The  first  order  fluid  equations  Lire  given  by: 

3d, 


(19) 

(20) 


♦  n0  1  •  2,  -  0, 

3V,  e.  V  0 


where 


n  a  3 
V  «  y 


3y  +  z  3z 


and  J  is  a  species  index,  i.e.,  J  **  i  for  ions  and  J  *  e  for  electrons.  These 
equations,  when  coupled  with  the  linearized  Poisson's  equation, 


(21)  ■  -  4tt  /Le.n  , 

j  J  J 

are  sufficient  to  obtain  the  dielectric  function  for  non- perpendicular  wavenuntoers . 
Equations  (19)  -  (21)  are  F curler  analyzed  in  space  and  time.  The  three  components 
of  the  F curler  analyzed  equation  (20)  are  solved  for  Vj  in  terms  of  the  scalar 
potential  <t>  and  the  perturbed  density  n^.  Then  equation  (19)  can  be  used  to 
express  the  density  in. terms  of  the  potential.  When  this  result  is  folded  into 


Poisson's  equation,  we  find  the  linear  dielectric  function  E(k,w),  given  by 


e(k,  «)  » 


sln^cj 

"  -  -  k%  +  ^ei) 
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The  normal  modes,  found  by  setting  e(k,  w)  =  0,  are  given  by 

2 


(25) 


“lveik  f 

“±  '  2  I2 

^Dl^ce  ' 


2“k(1  +  f>*“ce  3*,14H> 


l 


where 


f  S  1  - 


^a2  w2  ‘*2 


-El 


% 


v2  k2 
el  K 


(1  +  §*  \ 

K  01  ’ 


In  ootainlng  equation  (25),  we  have  made  use  of  the  Inequalities 

k^Tl/ne«ce«  l  and  k^  »  k2.  For  perpendicular  propagation,  a  *  0,  R  -KL, 
and  the  normal  modes  are  a>+  *  -  1  and  U_  -  -  1  k2  D  so  that  we  recover 
the  results  of  equations  (13)  and  (lH).  However,  for  finite  o,  R2  rapidly 
becomes  negative  and  the  roots  u>±  become  complex  which  corresponds  to  heavily 

danped  oscillations.  The  angle  a,  at  which  this  crossover  occurs  is  found  by 
setting  R2  =  0  or 


(26) 


kv 


_ei _  (1  +  ft  >-*s 

Ti  }  * 


pe 


Therefore,  as  mentioned  in  Section  I,  the  cross-field  diffusion  node  is  limited 
to  angles  less  than  10“^ 


i 


0 


radians. 


IV.  EFFECT  OF  THERMAL  FLUCTUATIOMS  ON  DIFFUSION  MODE 


As  we  have  seen  In  Section  II  the  diffusion  mode  exists  only  for  wave- 
vectors  nearly  perpendicular  to  the  external  field  B^.  Therefore,  in  examining 
the  thermal  effects  we  shall  assume  as  in  Section  I  that  k  •  f-  0.  Moreover, 
we  retain  Eqs.  (1)  and  (4)-(7)  but  not  Eqs.  (2)  and  (3).  We  consider  only 
the  first  order  electron  fluid  equations.  These  are  the  equation  of 
continuity. 


*r+no4-V°> 


the  equation  of  motion  in  the  y-direction, 

(28) 


,  30.  0  , 

«  _  _ *  o _ on  ,  Tf 

0  s;  ¥  ei  y 


0 


3  m  v  .  . 

J  e  el  3y 


32  V 

-ri 


+  —  +  w  v  , 

itL  <*y  ce  x* 
e 

the  equation  of  motion  in  the  x-directlon. 


(29)  0  -  -v  .  V  +  | 


ei 


30, 


ei  x  2  m  ui 
e  ce 


-5±  -  U  V  , 
3y  ce  y  ’ 


and  the  heat  balance  equation  for  electrons, 

30.  0  3V  0  v 

(30)  W  ^  ♦  |  e0  -  -  -s- 


ce 


3y 


n  0. 


32  0, 


mewce 


"2  Vei 


3y 


2m 

- -  V  A  » 

m^  el  1* 


where  n  ■  2/3(^ -66) . 

In  obtaining  Eq.(30),  we  assumed  equal  ion  and  electron  background 
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tenperatures  and  neglected  the  Ion  temperature  fluctuation.  The  viscous  terms 

in  Eq.  (29)  were  neglected  because  they  are  of  order  k2R  2  «  1.  For  the  same 

e 

reason  all  viscous  terns  except  one  are  neglected  in  Eq.  (28).  The  viscous 
heating  terms  in  Eq.  (30)  are  neglected  because  they  arc  quadratic  in  first 
order  quantities.  The  first  two  terms  on  the  right  hand  side  of  Eq.  (30)  are 
the  divergence  of  the  electron  heat  flux  and  the  last  term  is  the  collisional 
heating  of  electrons  by  ions  (for  equal  ion  and  electron  temperatures  and  no 

p  P 

ion  thermal  fluctuations).  This  last  assunption  is  Justified  provided  k  R.<:»1. 
In  this  regard,  we  note  that  Eq.  (3)  is  no  longer  necessary  as  viscous  terms 
have  been  included  where  they  are  inport  ant.  Equations  (2l)  and  (27)-(30)  are 
Fourier  transformed,  and  the  resulting  system  of  equations  solved  for  the  di¬ 
electric  function  e(k;w)  where 


(31) 
and 

(32) 


e(k,'»>)  =  1  +  k^/k2  +  Xe(k,w) 


(Xe) 


-1 


V 


u 


iA  2v 
e  ei 


+  1  + 


oj/6 


2m 


6w  +  ivel  +  (n-3/2)  ik2R„2v 


e  ei 


If  we  include  only  the  first  two  terms  in  curly  brackets  in  Ew.  (32),  we  recover 

the  resu-ts  of  Section  I.  The  additional  term  in  Eq.  (32)  is  due  to  thermal 

fluctuations;  and,  for  frequencies  u  -  ik2Revel,  it  is  clearly  not  negligible. 

The  viscous  terms  however  can  be  consistently  neglected  provided  A  2  «  1. 

e 

The  dispersion  equation,  e(k,u>)  *  0,  has  two  roots  u>  +  defined  by 

(33)  «±  -  n*  lk2Re2vel 

n±  -  _  1.9  -  Y  ±  [y2  -  3y  +  .tf* 

Y  -  .62  (~)A2R  2 
e 


♦ 
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Pot  the  ionosphere  u  -  and  k  -  0.1  k_  so  that  y  «  1.  Thus,  vdth 

UC  UC 

the  Inclusion  of  thermal  effects,  the  dispersion  relation  has  two  diffusion  modes 
where  there  was  one  (Eq.13)  without  thermal  effects.  However,  for  0e  =  0^  both 
these  modes  (Eq.  33)  differ  in  value  from  the  mode  uncorrected  by  thermal  effects 
(Eq.  13)  only  by  10J£.  Since  each  mode  still  has  zero  real  frequency,  we  don't 
expect  the  Inclusion  of  thermal  effects  to  alter  the  conclusions  of  Appendices 
E  or  F. 
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FIGURE  1 


FIRST  ORDER  DIFFUSION  IN  A  HOMOGENEOUS  PLASMA 


Ell 


FIGURE  2 


FOOTNOTES 
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APPENDIX  F 


STIMULATED  DIFFUSION  SCATTERING 


M.V.  Goldman  and  D.F.  DuBois 
University  of  Colorado 


Ia  this  note  we  discuss  s  new  electromagnetic  parametric  Instability 
In  which  an  O-wave  heater  acts  as  a  pump  to  drive  unstable  the  diffusion 
mode  described  in  the  previous  note  together  with  a  scattered  O-wave.  The 
threshold  and  growth  rates  for  this  Instability  (to  be  called  stimulated 
diffusion  scattering,  or  SDS)  are  computed  under  the  simplifying  assumption 
that  the  pump  linewldth  is  zero.  We  also  consider  possible  relevance  to 
spread-f . 

Figure  (la)  shows  the  simple  geometry  for  which  the  calculation  has 

first  been  performed,  The  incident  and  scattered  wave  vectors  are  orthogonal 

to  and  the  wave  vector  of  the  diffusion  mode  therefore  also  orthogonal 

(the  lowest  threshold  occurs  for  k..,_  .  *  2k.  as  shown).  Figure  (lb) 

aizzusion  me 

shows  a  more  realistic  stimulated  scattering  geometry  for  the  ionosphere  in 
which  k^^^  is  at  an  angle  to  the  geomagnetic  field,  and  is  such  that 

is  still  perpendicular  to  Bq.  We  do  not  expect  too  much  quantitative 
difference  between  (la)  and  (lb) ,  although  we  here  only  look  at  (la) . 

We  assume  the  validity  of  the  fluid  equations  with  pump,  A^  cos(d)ot  -  ky) 
and  scattered  radiation  A^.  These  are  combined  with  (ie  Poisson  and  wave 
equations.  The  vector  potentials  and  A^  belong  to  the  incident  and  reflected 
O-waves,  respectively.  The  use  of  a  monochromatic  pump  is  possibly  a 
serious  oversimplification  here.  In  the  fluid  equation,  the  velocity  variable 
is  related  to  canonical  momentum  by 


«  -  1  +  (*,  cos  <V  -  koy)  +  Aj)  . 


This  velocity  variable  u,  given  by  the  canonical  momentum  divided  by  the  mass 
has  been  used  by  many  and  perhaps  most  recently  by  Bodner  and  Eddleman*  in 
their  treatment  of  stimulated  Brillouin  scattering  (SBS),  which  is  similar 
to  our  treatment  of  SDS.  We  assume  that  the  zero  order  fluid  velocity  is  zero 

(u  -  0). 
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The  dominant  terms  In  the  fluid  equations  for  the  first  order 
electron  density  and  canonical  momentum  are  given  below. 


+  nQ  V  •  u^  *  0 


3t  o 


■  -  —  (u  x  B  )  +  ?n  +  v  .  L  -  —  A.1 
3t  me  — 1  -o  mn  —  1  el  fl  me  — 1 J 


-  -  [e  +  —  V 
m  1—1  me  — 


A  •  A  cos  (co  t  -  k  y)  I 
-o  -o  o  o'  J 


V  •  E L  -  4u(pe  +  p±) 


k  cos  (“o'  - ky)  +  °o±i] 

c  3t  c  m  L  -» 


The  momentum  equation  (2)  Includes  momentum  transfer  to  the  Ions  through 

electron-ion  collisions  and  the  force  associated  with  the  VA^  •  Aq 

non-linearity.  (From  a  Hamiltonian  point  of  view  this  force  arises  from  the 
2 

pA  coupling  which  adds  to  the  p$  coupling.)  Equations  (1)  and  (2)  tell 
os  that  the  electron  charge  density  responds  to  the  Indicated  nonlinearity 
In  exactly  the  way  that  it  linearly  responds  to  an  electrostatic  field  E^; 
that  is,  through  a  linear  electrostatic  susceptibility  (cf  our  first 
note).  When  this  information  is  put  into  Poisson's  eq.  (3),  we  obtain 


e(k.w>  E^k.cu)  -  !e  *1  <“  *  “o> 


in  which  the  pump  beats  against  the  scattered  wave  to  act  as  a  source  for 
driving  the  diffusion  mode.  Maxwell's  equation  for  the  scattered  wave,  driven 
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by  the  diamagnetic  current  Is 


F(w  -  h>  )2  -  w  2  -  c2  |k  -  k  l2  +  el  p  1  A.(w  -w  )  - 
[  O  p  -o'  jj)  -  (1)  J  1  o 
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In  this  current  the  pump  beats  against  the  first  order  charge  density  which 
is  again  expressed  via  the  linear  susceptibility  Xfi  and  leads  to  equation  (6). 

Equations  (5)  and  (6)  comprise  the  coupled  equations  for  the  electro¬ 
static  field  E^  associated  with  the  diffusion  mode  and  the  vector  potential 
associated  with  the  scattered  0-wave.  The  regenerative  interaction 
of  these  two  modes  through  the  pump  leads  to  the  Instability.  The 

effects  of  E  are  contained  in  the  conventional  driven  velocity  v  .  The 
— o  O 

linear  susceptibility  Xg  which  is  essentially  the  coupling  coefficient  is 
derived  in  Appendix  E.  The  results  of  the  linear  theory  are  briefly  summarized, 
and  the  diffusion  mode  dispersion  relation  Is  written  as 

e  ■  1  +  — 5—  +  I  .  Uncoupled  diff.  mode: 
k 1  e 


E  -  0  «  u  -  i-  .  -1»  A  !  (1  +  A 

D  e 


(7) 


where  Is  the  diffusion  time. 

The  scattered  radiation  consists  generally  of  a  downshifted  Stokes 
frequency 


to  -  to 

s  o 


(1) 


c2(k 
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and  an  upshifted  anit-Stokes  frequency 


0)  +  W  a 
o 


The  wave  vector  k  -  k  is  associated  with  the  Stokes  frequency,  and  k  +  k 
~ °  —i o  — 

with  the  anti-Stokes  frequency.  Since  the  frequency  U)  is  assumed  «  u 

o 

the  pump,  and  Stokes  frequencies  are  approximately  equal  and,  hence,  we  must 
have 


|k  -  k|  *  k  . 

1  — o  1  o 

This  is  possible  only  for  k  *  2kQ  or  k  «  kQ.  The  former  is  backward  scattering 
and  the  latter  is  forward  scattering.  If  the  anti-Stokes  mode  is  to 
participate,  w  cannot  be  far  from  a)  ,  and  then  we  must  also  have 

A  O 

|k  +  k|  "  k  for  anti  Stokes. 

— o  —  o 


This  is  obviously  only  possible  for  forward  scattering  which  we  do  not 
treat  here.  For  backward  scattering  the  anti-Stokes  mode  would  be  too  far 
off  its  linear  frequency  to  participate,  and  we  have  the  strictly  3-mode 
interaction  described  by  eq.  (3),  (5)  and  (6). 

The  threshold  obtained  from  solving  the  dispersion  relation  implied  in 
equations  (5)  and  (6)  is  given  by 


+  1)  (1  + 
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For  the  Ionosphere,  with  *  0^,  the  threshold  value  of  /vg  (directed 
velocity  squared  to  thermal  velocity  squared)  is  12  ^2*  Ve^/a>o.  This  Is 
about  a  factor  of  two  higher  than  the  threshold  for  stimulated  Brlllouin 
scattering  or  the  parametric  decay  instability.  Both  SDS  and  SBS  have  an 
advantage  over  the  decay  instability  in  the  sense  that  neither  SDS  nor 
SBS  involves  Languir  waves  and  Landau  damping.  This  leads  to  a  much  wider 
range  of  altitudes  for  which  these  instabilities  can  in  principle  occur. 

We  ct mpute  the  threshold  below  for  a  typical  set  of  parameters: 


(  v  2/v  2)  .  -  1.3  x  10  ^  for 

D  e  threshold 


V  .  ■  300  Hz,  co  ■  2nf  ■  A  .  10  radians, 
el  o  o 


On  the  other  hand,  a  heater  with  local  vacuum  intensity  of  50  vm/wT 
Incident  upon  a  plane-layered  plasma  with  scale  length  100  km  gives  rise  to 
a  maximum  local  value  of 


v  2/v  2  -  2  x  10”3, 

1st  max 

which  is  15  times  the  threshold  for  this  instability. 

The  frequency-matching  condition  for  SDS  is  given  by 
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|k  -  k  |  =  k  ->  k  =  2k- 
1 —  — o  o  —  ~0 


and  leads  to  the  conclusion  that  the  instability  is  strongest  when  the 

wave  number  of  the  low  frequency  mode  is  twice  that  of  the  incident  O-wave. 

There  is  some  latitude  in  the  frequency-matching  condition  insofar  as  it 

affects  the  threshold,  but  since  the  difference  between  u)  and  a.  is 

o  T 

always  smaller  than  either  frequency,  the  momentum-matching  condition 
k  -  2k^  remains  approximately  satisfied. 

The  dispersion  relation  for  the  scattered  wave  which  corresponds 
to-mlnimum  threshold  (or  maximum  growth  rate)  turns  out  to  be  exactly  the 
uncoupled  dispersion  relation: 

“  -  (AjCk  -  kj)  i  u>p2,+  c2|k  -  kj2 

Thus,  o)q  -  0)  is  a  real  frequency  given  by  1/VT  .  The  coupling  has  there¬ 
fore  converted  the  diffusion  mode  into  a  wave— likemode  of  (low)  frequency 
qf  order  of  the  former  diffusion  rate.  This  situation  is  exhibited  in 
Figure  2  which  shows  the  uncoupled  normal  modes  (dots)  and  the  coupled 
(new)  normal  modes  (crosses) . 

Two  expressions  are  given  for  the  growth  rates  above  threshold  associated 
with  SDS: 
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The  first  is  valid  only  for  growth  rates  less  than  Tp  >  and  the  second  on  y 
for  growth  rates  larger  than  vel.  Me  have  not  ret  completed  the  numerical 
work  necessary  to  give  the  growth  rates  in  the  intermediate  range.  To 
compute  the  growth  rate  we  first  need  the  diffusion  rate 


—  -  2u  ,k2R  2. 
XD  ei  e 


The  wave  number  k  is  **ice  the  local  pump  wave  number.  We  look  first  near 
the  first  Airy  maximum,  where  the  electric  field  vector  of  the  pump  is  largest. 
There  the  local  pump  wave  number  is  related  to  the  free  space  wave  number  by 


local 


-c  2  \ul) 


where,  L  is  the  density  scale  length.  For  L  - 
radians,  eq.  (13)  gives 


100  km  and  ui  -  4  x  10 
o 


k  ■  .066  —  . 

0 local  C 


The  diffusion  rate  is  then  given  by 
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(15) 


—  -  8.3  x  10~*  Hz,  or  T  =  20  minutes 
XD  ° 

This  extraordinarily  slow  rate  implies  that  a  more  accurate  theory  is  required 

-3 

taking  into  account  the  finite  pump  linewidth  when  it  exceeds  10  Hz. 

Inserting  (15)  into  equation  (10)  we  find  for  the  growth  rate, 

Imbi  ■  5.2 

XD 

which  violates  the  conditions  of  validity  for  equation  (10),  but  is  probably 

roughly  accurate  (eqn.  (11)  gives  a  far  worse  violation  of  validity 

conditions).  This  leads  to  a  growth  time  of  about  4  minutes,  which  is 

rather  long  in  comparison  with  common  spread-f  onset  times. 

-1  2 

Since  the  growth  rate  goes  as  and  hence  as  k  local»  we  nave 

examined  the  growth  rate  at  somewhat  lower  altitudes,  where  the  ratio  of 

the  local  pump  wave  number  to  the  free-space  wave  number  is  larger  than  .066. 

However,  the  Airy  function  swelling  factor  for  the  pump  is  smaller  at  these 

altitudes  and  the  gain  in  wave  number  appears  to  be  generally  offset  in 

intensity,  so  that  the  magnitude  of  the  growth  rate  is  not  appreciably 

changed  from  that  at  higher  altitudes.  (As  an  example,  at  the  100th  Airy 

node  k2  is  22  times  larger  than  at  the  first  node,  but  the  swelling 

’  local 

factor  is  down  by  a  factor  of  10  and  the  growth  rate  is  numerically  unchanged.) 

Work  is  in  progress  on  several  other  features  of  this  instability. 

In  addition  to  folding  in  a  finite  pump  bandwidth,  the  competition  of  this 
instability  with  other  parametric  instabilities,  such  as  SBS  and  the  decay 
instability,  must  be  studied.  Stimulated  Brillouin  scattering  may  be 
especially  interesting,  because  the  growth  rates  are  larger,  although  the 
unstable  ion-acoustic  waves  are  not.  field  aligned. 
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NONLINEAR  SATURATION  OF  PARAMETRICALLY  EXCITED 
FIEID  ALIGNED  IRREGULARITIES 


Jeraae  Weinatock  and  Bandel  Bezzerides 


Goldman  and  Dubois  have  shown  that  a  low  frequency  unstable  wave  propagat¬ 
ing  across  the  magnetic  field  la  excited  during  modifications  (see  these 
proceedings )  •  Hero  we  calculate  the  nonlinear  saturated  amplitude  of  this 
wave •  The  parametric  excitation  of  this  wave  can  be  Illustrated  as  follows? 


low  frequency  diffusion  mode 


0-Mode  (Scattered) 


Low  Frequency  Mode  at  Threshold 


«  -  k2  D  -  1  k2  D  , 
l  xa  i  ia 


(D 


where  Is  the  wavenumber  and  Is  typically  of  the  order  of  tens  of  meters  $ 
D  n  is  the  amblpolar  diffusion  coefficient  normal  io  the  magnetic  field. 


v  r8 
e  e 


1.  c?_ 


R) 


with  ve,  f>e,  ve  the  electron  collision  frequency,  gryofrequency  and  thermal 
speed. 
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The  linear  growth  rate  /  of  the  mode  is  roughly  given  by 


-(k)  „  [p  f(lt)  .  i 


k2  D  ,  2coa29  *  1 


where  P  is  the  ratio  of  pump  intensity  to  threshold  intensity,  0  is  an  angle 
in  the  plane  perpendicular  to  the  magnetic  field  between  the  wave  vector  and 
the  vertical  plane,  f(k)  is  unity  for  the  matching  k  and  is  very  narrowly 
peaked 

The  nonlinear  wave  equation  for  this  mode  is  given  by 


sf 


(/(h)  -  /L(k>) 


(3) 


where  1^  ie  the  spectral  intensity  at  wave  vector  k,  y  is  the  nonlinear 
damping  decrement  due  to  waves,  and  S  is  the  noise  source. 

/V)  lo  determined  by  calculating  the  porturbation  of  oloctron  and  ion 
orbits  due  by  the  unstable  waves.  This  result  is  found  to  be  [see  Eq's.  (39) 
and  (52)  of  We  instock  and  Williams,  Fhys.  Fluids  l4,  1^72  (l97l)]. 


where  c  is  the  speed  of  light,  B  is  the  magnetic  field  strength,  z  is  a  unit 
vector  along  the  magnetic  field,  and  is  the  real  part  ?f  the  frequency  of 
the  unstable  mode.  In  the  present  case 

u)Rw'/ff,(k')  (5) 

To  approximately  solve  the  wave  Eq.  {3)  for  the  saturated  value  of  1^ 
we  set  dl^/dt  «  0,  and  neglect  S  as  being  relatively  small  to  obtain 

V®‘(k)  «  /(h)  (6) 
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Substituting  (2)  into  (6)  ve  bare 


/L(  )  „  L  Scoafe^l  fM  _  1  fc2 
L  cos20  J  A 


Equations  (4)  and  (7)  can  be  solved  for  the  total  wave  energy  density 

V 

f  ^  "r  *  ^b®  solution  is 

J(2tr)3  m  .  _  _2y.  2  _  \2  J2. 
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vbere  (<k^))  ■  is  the  mean  ware  vector  component  in  the  norml  plane  defined  by 


«k j»  ■  k2  r 

X  X  j  I 


p  p 

sin  0  cos  0 
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as  is  the  electron  plasma  frequency.  The  value  of  ((k  ))  •#  is  found  to  be 

pe  * 

roughly  given  by 

«kj;»  wio’1  k^  (1C 

The  mean  square  density  fluctuation  (6n)  at  saturation  can  be  obtained 
from  (8)  by  making  use  of  the  relation 

1  mkT\^  <u 

“  *D 


so  that 


<|E?I>  -  0*™0a  *5  <h£l> 

*D 

Substituting  (18)  Into  (8),  with  <E^>  -  1(h)  *  Volum,  w»  have  finally 
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Equation  (13)  is  the,  mean  square  density  fluctuation.  For  typical  values  of 


ftj 


w  10 


P  »  12 


<(k^>>  “  10 


ve  find,  from  (13),  that 


(SJ 


?  -4 

«  1.1  x  10 


We  note,  from  (13 ),  that  the  level  of  density  fluctuations  depends  upon  the 
angular  spectrum,  hut  does  not  depend  upon  the  absolute  value  of  the  wavelength. 
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